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water can be stored during times of over-supply for consumption during the 
dry periods. Most of the artificial storage is provided in the form of surface 
reservoirs. These reservoirs change the regime of the rivers drastically, 
particularly with respect to their sediment characteristics: 


1. Any large storage reservoir will permanently store practically the 
entire sediment load of the stream. 

2. Part of the sediment storage occurs in the storage volume proper, part 
in the channel and valley bottom upstream by backwater effects. 

3. The sediment load leaving the reservoir is usually nil or is restricted 
to very much reduced releases that have the purpose of flushing some of the 
sediment through the storage volume. 

4. The hydrograph of flows available in the stream channel downstream 
from the reservoir is entirely changed and with it the annual sediment load 
carrying capacity. 


Problems in and Above the Reservoirs.—Sediment deposition inthe reser- 
voir causes loss of storage volume; deposition above the reservoir causes 
back-water effects, swamping and flooding in the valley. These damages have 
hitherto been accepted as unavoidable and toalarge degree as uncontrollable. 
Many examples of actual reservoirs and lakes have shown that the deposition 
of the sediment may follow many different patterns. Often, changes of patterns 
are clearly caused by particular influences, such as snags or vegetation, that 
suggest possible engineering control of the deposition process. If it is realized 
that both the storage volume in the reservoirs and the area of the extended 
valley bottoms above the reservoirs represent large and often irreplaceable 
values, it appears to be extremely wasteful to accept these losses, as they 
may occur, without question. A planned program of valley sedimentation, on 
the other hand, may not only reduce the loss of storage space in the reservoir 
but also create large valuable valley areas that are sufficiently high to pre- 
vent flood damage if used agriculturally or industrially. Such a program may 
also prevent the development of swamps with their health hazards, 

What type of research is needed in connection with the development of 
such areas for sediment deposition? The first stepis probably the systematic 
observation and description of the various deposition processes as they occur 
in existing prototype conditions, covering all materials from the bed load to 
the finest, for an extended range of river conditions. However, even such a 
very extended survey is only able to describe various chance combinations of 
the many variables that may or may not be important in their effect on the 
final result. Only asystematic study of the problem under the controlled condi- 
tions of a laboratory will permit the separationof the effects of the individual 
parameters, 

Such laboratory studies, in turn, become effective only after the similarity 
conditions for deposition processes are known. Today, most model studies 
with movable bed deal with an approximately uniform bed material and try to. 
find the geometric patterns of scour, transport, and deposition for that material 
under given flow conditions. The study of reservoir deposition would need the 
inclusion of all the fine components of the load in order to study the more 
complicated process of total load deposition. 

After the development of methods operating with a large range of grain 
sizes under deposition conditions, and after scale effects for such processes 
are known, one can study specific ‘questions: 
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Under which flow conditions is sediment deposited in a dense form; under 
which in a loose form? Is strength, or resistance against shear, in a sediment 
deposit related to density? 

Under which flow conditions are sediment deposits segregated as to grain 
sizes and under which are mixtures deposited? 

What are the relationships between range of grain sizes, density, and 
strength of various deposition patterns? 

What are the agricultural values of the various types of deposits? 

Can a systematic variation of the flow conditions achieve deposits that a 
continuous flow is not able to produce? 


These are some of the most pertinent questions in connection with reser- 
voir and valley deposits. Such questions may arise in connection with the 
engineering control of large-scale deposition processes and maybe answered 
by systematic studies. It is certain that many more such problems will arise 
as soon as either deposition studies or deposition engineering in the prototype 
are undertaken. 

The entire field of engineering control of deposition processes in and 
above reservoirs is practically unknown today. It is unquestionably possible 
to develop rules, methods and other tools for that purpose in the same way 
that similar know-how has been developed for the deposition of sediment in 
near-equilibrium streams for the purpose of channel stabilization. 

There are mainly two reasons for the lack of such development. First, 
most of the reservoirs of significant size are located upstream of the area in 
which, at the time of construction, valley development was important. There 
is usually not enough damage expected in the reservoir proper for a genera- 
tion or two to-warrant the expenditure of funds and effort to prevent it, That 
means that the incentive is missing in cases in which new areas are developed. 
Secondly, where reservoirs are built into an already developed area, the local 
landowners above the reservoir do notbenefitfrom it, but usually incur a loss 
of usuable land in various ways, and always feel that they are not sufficiently 
compensated for such losses. The operators of the reservoir, private or public, 
are then in the precarious position of having hostile neighbors in the entire 
area in which such controlled deposition may be practiced. 

It is sometimes difficult for an individual to prove legally that he has been 
damaged by a very slowly progressing deposition process. In contrast, if 
engineering measures are undertaken to control deposition in the valley bot- 
tom, many local owners who may actually benefit from the operations, will go 
to court and claim damages. Thisis arather costly undertaking for all parties 
involved. It can probably be avoided by organizing all landowners into a district 
that can then as a unit negotiate with the operator of the reservoir, after a 
concrete plan has been established for the operationof the deposition process. 
Unfortunately, sufficient general knowledge is not available today on these 
processes on which such a plan canbe based. The economy of any such plan is 
usually even more in question. Only specific research in this field can provide 
the necessary background for such developments. 

Sediment Problems Downstream from the Reservoirs.—The construction of 
a reservoir will have its greatest benefits in the river valley downstream 
through flood control and increased all-year water supply. The river valley 
downstream must be expected to develop economically at a much increased 
‘rate whether it remains agricultural or changes to urban and industrial de- 

velopment. The river that brings the fresh water from the reservoir will also 
t 
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be called on to drain the valley bottom and, depending on its size, may even 
provide inexpensive transportation, Because the river channel is already 
provided, the laymen are apt to accept its shape and course as given and 
stable. When the majority of the high flows are held back by the reservoir, 
most or all of the flood danger appears to be eliminated. 

To the engineer who is familiar with the management of alluvial rivers, 
the picture looks less rosy. He recognizes that the construction and operation 
of the reservoir severely reduces the sediment supply and radically changes 
the hydrograph of the river, These changes are so drastic that in most cases 
there is very little similarity between the new and the pre-reservoir river. 
Only the channel is the same. It is a sediment-lined channel that had usually 
reached at least near-equilibrium shape for the old river condition, It would 
be an extreme coincidence if the channel were still in equilibrium at the new 
flow and sediment conditions. The sediment supply is reduced. But also the 
channel capacity to transport sediment is somewhat reduced by the reduction 
of flood flows and the increase of medium flow durations. Usually the sediment 
supply is reduced more than the carrying capacity, causing the channel to 
erode. Often this development may be reversed in lower reaches after some 
sediment carrying tributaries have joined the main stream with the reduced 
sediment capacity. The immediate response of the river will be scour in the 
upper reaches. It will have acquired a fullload when it arrives at the junction 
with any sediment carrying tributary. Itmay beginto aggrade at and below that 
junction, If the river channel is left free to develop at will, the upper reaches 
will, after a given time, adjust themselves to the new conditions. They will 
stabilize at a lower sediment rate, andthe sediment load to the lower reaches 
will be reduced. It must be expected that during that process of adjustments 
the sediment supply to all parts of the river system will steadily change and 
call for continued readjustments. 

After a channel section has finally stabilized its profile, for instance by 
coarsening its bed sediment, it is, unfortunately, not yet completely stable. 
The banks, that usually consist of a finer material than the bed, also begin 
to fail. They also seem to be stabilized only as long as finer material is 
deposited from the stream at about the same rate at which it scours, With 
the finer sizes missing in the load, the banks retreat steadily and the river 
channel widens. The flows begins to meander between banks and attack them 
with increased force where they impinge on them. It is not sufficient to sta- 
bilize the bed alone; the same must be done with the banks. : 

The valley bottom below a reservoir has been described as the area in 
which the maximum economic development must be expected to take place. 
The valuable ground will be in great demand and developments of all kind 
will soon begin to encroach on the river, A gradually enlarging river channel 
that is usually already too wide for the reduced flood flows is, therefore, a 
very undesirable solution. 

Under these changed flow conditions, the engineer will ask again a number 
of questions when be tries to design a more desirable channel below a major 
reservoir. | 


What methods may be used to stabilize the profile of various streams 
downstream from reservoirs? 


What methods may be used to stabilize the banks of various streams pelowl 
reservoirs ? 


What cross section will be stable for such a changed stream? 
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What changes will occur in the composition of the stream bed due to the 
changed flow conditions? 

Is it possible to stabilize channels for a range of sediment supply rates 
for each flow? 

Is it more advantageous to stabilize the river channel before, immediately 
after, or any time after closure of the dam? 

What types of structures are most effective for the stabilization of such 
channels and which materials should be used in their design? 

Should stabilization begin at the dam and progress downstream or should 
it begin simultaneously over a considerable length of stream? 

Should the changed stream after stabilization follow a meander pattern 
Similar to the original stream or may it be straightened? 


Many more such questions can and will be asked, particularly if special 
interests such as navigation, water quality control, and others are involved. 

Now one must translate these engineering questions into research problems. 
This has already been attempted at the University of California, Berkeley, 
Calif. sediment laboratory for a number of years, Answers have been obtained 
from the resulting studies that, in many instances, point the way to the solu- 
tions of practical problems. First, the general conclusion was drawn from the 
previously stated arguments that any designed and stabilized channel must 
have artificially stabilized banks. The width between these banks is not auto- 
matically attained by the river, but mustbe chosen by the designer. This is in 
strong contrast to the alluvial river in its natural equilibrium state, in which 
the width reaches an equilibrium value depending on parameters of the flow, 
the sediment load, the vegetation and other factors. 

The basic problem of the changed river as it arises in its most general 
form can then be stated as follows: Given the hydrograph and the sediment 
supply by amount and composition for all times, what slope, bed composition, 
depths and flow velocities will result as a long-range equilibrium for a given 
chosen channel width? Or the slope may be given and the channel width 
required. This basic question in its bold generality has not been solved, not 
even approximately. But a first step to its solution has been made at the 
University of California. 

An equilibrium flow with a very heterogeneous load was disturbed by both 
an increase and then by a decrease of the sediment load. It was found that an 
increased load will behave quite rationally as may be expected from trans- 
portation laws for uniform sediment. The decreased load, on the other hand, 
resulted in a strong coarsening of the bed with progressively reduced scour. 
It was found that for each flow there exists a range of sediment sizes that 
move more slowly than the bulk and that all sizes above a limit practically 
do not move at all. These “slow-moving” and “non-moving” sizes appear to 
be the tools that nature uses to retard and even to prevent excessive down- 
ward scouring of channels. This process of bottom shingling and the resulting 
very limited bed scour has been observed in various rivers below large 
reservoirs. More work is needed for its quantitative prediction, particularly 
on the effect of flow variability. ; : 

In many rivers such non-moving sediment sizes are not contained in the 
river bed in sufficient amounts to limit down-scour to an acceptable depth. 
How can a similar protective layer of loose stone be introduced into the gabe 
sed artificially to limit down-scour? How does the necessary amount of ae 
Jepend on the rock size? These are all most important problems that see = 
solved experimentally and for which a complete answer does not exist. 
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In river reaches in which there is still an ample sediment supply in the 
bed sizes, but a lack of fines, the banks may need stabilization whereas the 
bed is Still in equilibrium. The question arises in such rivers how deep the 
foundation of the banks must go. In any event must the toe of the bank be be- 
low the deepest instantaneous scour holes ofthe bed that may develop near the 
bank, This calls for a systematic and quantitative study of the shape of bed 
irregularities for the various flow and bed conditions. This extremely com- 
plex problem has recently found some interest in various quarters, but no 
theory is developed as yet by which quantitative river predictions can be made. 

These are some of the most important developments that may be expected 
to occur in our rivers and an attempt has been made to predict the most im- 
portant research problems that will arise with them. It is most interesting to 
remember that many of our notorious sedment carriers were clear and deep, 
quietly flowing, and tree lined streams at the time agriculture was brought to 
this land. They were changed into what they are by reckless methods of soil 
exploitation. It may now be necessary to change many such streams back to 
their original state because their waters have become sufficiently precious 
to conserve them. 


CONCLUSION 


The large storage reservoir is seen as the most important causefor new 
sediment problems is our rivers. Sediment transportation research should 
concentrate on the study of river conditions upstream and downstream from 
such reservoirs, 
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COST ALLOCATION OF WATER PROJECTS IN CALIFORNIA 


By Amalio Gomez,! F. ASCE 


SYNOPSIS 


This paper describes briefly the separable costs- remaining benefits method 
of allocating the cost of multiple-purpose water projects; summarizes the 
allocations made by the Corps of Engineers for several reservoir projects in 
the Central Valley of California; compares the results obtained; and makes 
some observations concerning the application of the method, based on actual 


experience. 


INTRODUCTION 


Cost allocation is as old as civilization itself. It probably began with the 
problem of sharing the cost of making war or of running the government. It 
may be said that cost allocation is a new name for describing a process quite 
similar to that of assessing taxes. It is now used primarily to describe the 
process of distributing the cost of amultiple-purpose water project among the 
functions that the project serves. It must be remembered that no fully satis- 
factory procedure has ever been devised for assessing taxes equitably in all 
cases. Similarly, no fully satisfactory procedure has ever been devised that 
will automatically yield equitable results in allocating the cost of a multiple- 
purpose water project, in all cases. The problem, therefore, concerns the 
computation of approximate values that reasonable men consider fair and just. 
Sound and experienced judgment is still one of the basic ingredients of any 


cost allocation. 


To extend the closing date one month, 
CE, This paper is part 
an So- 


_ Note.—Discussion open until August 1, 1961. 
a written request must be filed with the Executive Secretary, AS 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the Americ 


ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961. : 
1 Che., pieaniay and Reports Branch, Engrg. Div., U. S. Army Engrg. Dist., Sacra- 


‘mento, Calif, 
; 
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Experience in allocating the cost of multiple-purpose water projects goes 
back some 30 yr. Many methods have been devised and used with more or 
less satisfactory results. Some of the most common methods used include:2 


Proportionate benefit method 
Proportionate use method 

Priority of use method 

Alternative justifiable expenditure method 
Separate projects method 

Equal apportionment method 

Incremental method 

Direct cost method 

Separable cost-remaining benefits method 


METHOD 


The method commonly in use by Federal agencies is the separable costs- 
remaining benefits method. The underlying principle of the method is that all 


TABLE 1.—COST ALLOCATION 


Flood 
control, 
in 
$1,000 
» (2) 


Item 


1. Benefits 
2. Alternative cost 


3, Benefits limited by alternative 
cost (lesser of items 1 and 2) 


. Separable costs 
Remaining benefits (items 3-4) 
Allocated residual cost 

Total allocation (items 4+ 6) 


ee ee 


project functions should share equitably in the savings to be realized by the 
multiple-purpose development. In the application of the method, project bene- 
fits for each function are first limited by the cheapest alternative cost of ob- 
taining the same benefits. The method then assigns to each function all 
“separable cost,” that is, the costs that were incurred by virtue of adding each 
function to the project. The remaining costs are called “joint costs,” or 
“residual costs,” and are distributed among the various project functions in 
proportion to the benefits remaining after deducting each separable cost from 
each project benefit. The method is best described in a report to the Inter- 


2 «Report on the Allocation of Costs of Federal Water Resources Development Pro- 
jects,” House Committee Print No. 23, 82nd Cong., 2nd Session, U.S. Govt, Printing Of- 
fice, Washington 25, D. C. 
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Agency Committee on Water Resources.3 This report was prepared by the 
Sub-Committee on Evaluation Standards, that developed the separable costs- 
remaining benefits method. A fair understanding of the method can be gained 
from the following quotation and Table 1 taken from the report: 


“The recommended method of cost allocation is illustrated below 
for a multi-purpose project for which the total project costs amount 
to $1,765,000. These include investment costs and operation, main- 
tenance, and replacement costs, all reduced to a common time basis, 
and are expressed either as an average annual amount or a present 
worth amount.” 


APPLICATION 


Present practice of the Corps of Engineers is to make preliminary cost 
allocations in connection with preauthorization reports. The preliminary cost 
allocation is brought up-to-date when funds for initiation of construction are 
being requested, and again as soon aspractical after construction is initiated. 


TABLE 2,—CONSTRUCTION STATUS 


Reservoir Project 


(1) 


Pine Flat Kings River Completed 
Isabella Kern River Completed 
Folsom American River Completed 
Success Tule River Under construction 
Terminus Kaweah River Under construction 
New Hogan Calaveras River Under construction 


Black Butte Stony Creek Under construction 


The final cost allocation is prepared just before the project is placed in opera- 
tion. The practice followed during the last few years varied somewhat with 
each project, but, generally, it was to make a preliminary cost allocation 
before initiation of construction and fix the percentage distribution of project 
cost among its various functions at that time. Later, these percentages were 
applied to the actual construction costs. 


PROJECT FUNCTIONS 


Multiple-purpose reservoir projects built by the Corps of Engineers in the 
Central Valley of California serve two or more of the following fais 
(a) navigation; (b) flood control; (c) irrigation; (a) municipal and eae 
water supply; (e) hydroelectric power; (f) recreation; and (g) fish and wild e. 
Cost allocations for Corps of Engineers projects are made by the Chief of 
Engineers, generally with the concurrence of the Bureau of Reclamation and 
other affected Federal, State, and local agencies, With some important ex- 


_ 3 «proposed Practices for Economic Analysis of River Basin Projects,” U. S. Govt. 
Printing Office, Washington 25, D. C., May, 1958, 
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ceptions, costs allocated to navigation and flood control are borne by the 
Federal Government. Costs allocated to irrigation are paid for by the water 
users, generally without interest, over a 40-yr period, under a contract with 
the Secretary of the Interior, made pursuant to the provisions of reclamation 


TABLE 3.—SUMMARY OF ECONOMIC 


Pine Flat, Kings River | Isabella, Kern River 


Reservoir capacity 
Flood-control storage 
Irrigation storage 

Inactive storage 

Average annual runoff 

Average annual new water 

Ave, ann, new and rereg, water 


First cost 
Interest and amortization 
Maintenance cost 
Total annual cost 


Flood-control benefits 
Irrigation benefits 
Power benefits 


Total annual benefits 


Flood control 
Irrigation 
Power 

Total 


Flood control 
Irrigation 
Power 

Total 


; Includes $590,000 for $59,000 acre-feet of municipal and industrial water supply. 
b After deducting $2,140,000 creditable to generating, switching, and transmission 
facilities built by the U.S. Bureau of Reclamation, the cost of which is not included in 
the cost to be allocated. 
C Includes $8,810,000 allocated to municipal and industrial water supply. 
Exclusive of cost of power facilities built by the Bureau of Reclamation, 
© Includes $26,000 allocated to municipal and industrial water supply. 


law, including the 160-acre limitation and other restrictions. Costs allocated 
to municipal and industrial water supply are repaid with interest, under a 
contract approved by the Secretary of the Army. Costs allocated to hydro- 


electric power are repaid with interest. The power is sold by the Secretary 
of the Interior. 


Le a ae a ee 
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Until recently, recreation and fish and wildlife have been considered as 
important, but secondary byproducts of the primary project functions, How- 
ever, recreation and fish and wildlife are becoming more important evcry 
day in water resources planning. The Fish and Wildlife Coordination Act of 


DATA AND ALLOCATED COSTS 


Folsom, Ameri- 
can River 

(4) 

1,000 acre-feet 


Success, Tule 
River 


(5) 


Terminus, Ka- 
weah River 


(6) 


New Hogan, Ca- 
laveras River 


Black Butte, 
Stony Creek 


1958 (P.L. 85-624) provided that in certain pureed eles and eh ai 

projects, fish and wildlife will be considered for possible inc ue oe R a 
_ ject function on the same basis as any other project function, an a piri 
allocated thereto may be considered wholly or in part arenas Pus 
the projects discussed here, no project costs have been allocate 
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tion and fish and wildlife so far, but it may be that in some cases when the 
final cost allocations are made, some costs will be charged against these 
functions. 


COST ALLOCATIONS MADE 


Cost allocations have been made by the Corps of Engineers for seven 
multiple-purpose reservoir projects in the Central Valley of California. 
The names of these projects and the construction status are listed in Table 2. 

The Pine Flat cost allocation was not based on any specific method of cost 
allocation, but gave weight to the results obtained by the use of several 


TABLE 4,—COMPARISON OF COST ALLOCATION TO IRRIGATION, 


Project 
een Pine Isa- Ter- | New | Black 
Flat bella | Folsom| Success | minus} Hogan | Butte 
(2) (3) (4) (5) (6) (7 


Allocated first cost per 
acre-foot of average annual 


runoff, in dollars 40.52) 17.06 


Allocated first cost per 
acre-foot of new and re- 


regulated water, in dollars 88.36 |128.07 


Allocated first cost per 
acre-foot of reservoir 
storage usable for 


irrigation, in dollars 20.52] 48.67 


Annual cost to water users 

on basis of repaying capital 

cost in 40 equal payments 

without interest but includ- 

ing M. O., and R, charges: 
Per acre-foot of new 
water, in dollars 


Per acre-foot of new and 
reregulated water, in 
dollars 


methods. All other allocations have been made by the separable costs- 
remaining benefits method. Table 3 contains some descriptive data on the use 
of storage space as well as a summary of first and annual costs for each 
project. It also contains a summary of allocated first costs and maintenance, 
operation, and replacement costs. A further summary for each specific pro- _ 
ject with pertinent comments is given later. 

It is extremely difficult to pass judgment onthe fairness of a given alloca- _ 
tion. For this reason, results obtained at one project are often compared with 
results obtained at other projects. Such comparisons can serve a useful pur- 
pose. However, they must be properly qualified because there are too many 
unrelated factors that affect the results. Among these factors are: (a) physical 
characteristics of the dam and reservoir site; (b) characteristics of the 
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stream such as runoff and extent to which already developed; (c) service area 
in which water is used; (d) developments in flood plain; and (e) the price level 
under which the project was constructed. With these qualifications in mind, 
the comparisons involving the cost allocated to irrigation given in Table 4, 
may be of some interest. 

The preceding comparisons have been made on the basis of cost allocated 
to irrigation, because such costs are repayable and there is a tendency to 
judge the fairness of the allocation on the basis of repayable costs. 


PINE FLAT DAM AND RESERVOIR 


The cost allocation for the Pine Flat Project4 was made in 1947. It was 
based on an estimated first cost of $35,760,000, of which $33,500,000 was 
for the dam and reservoir, and $2,260,000 was for penstocks for future 
power generation at the site. The allocation was as follows: 


Floddteontrol i227 .> . $19,250,000 
TeripatiOne sh. cae e ste. ot tas ae 14,250,000 
1S Sieee 5 C clee Be eo a 2,260,000 

SLOtAlS Seer hee. ak $35,760 ,000 


The War Department Civil Appropriation Act, 1947, appropriatedfunds for 
the initiation of construction of this project subject to the proviso that: 


“, . . None of the appropriations for the Kings Riverand Tulare Lake 
project, California, shall be used ... until the Secretary of War... 
with the concurrence of the Secretary of the Interior, shall have made 
a determination as to what the allocation shall be: provided further... 
that the agreement of concurrence shallbe made not later than 9 months 
from the date of the enactment of this Act.” 


Because of the increase in price level and some design changes, when the 

dam was finished in 1954, the total cost was about $39,385,000, of which about 
$1,000,000 was for penstocks for future power. Because the allocation had 
been made at the specific request of Congress, and it had been furnished to 
and accepted by Congress, it could not be revised or done over again, on the 
basis of the actual construction cost, without a new directive from the Con- 
gress. This procedure has not been repeated with any of the subsequent pro- 
jects, 
The project is being operated for irrigation on the basis of interim con- 
tracts between the water users and the Bureau of Reclamation. A permanent 
contract is expected in the near future. The project has been eminently suc- 
cessful as illustrated by the following accomplishments during the 6-yr period 
since its completion in 1954: 

Cash payments received from water users. ....- - $4,627,000 

Cash payments received for storage of power water. - $704, 000 

Estimated flood damages prevented. ....-- + + + $20,000,000 


ISABELLA DAM AND RESERVOIR 


The cost allocation for Isabella Dam and Reservoir was made in 1955, 
when construction of the project was essentially completed. In this allocation, 


i a Printed in House Document 136, 80th Congress, 1st session U. S. Gov’t. Printing 


Office, Washington 25, D. C. 
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annual cost and benefit data were capitalized on the basis of 2.5% Federal 
interest rate and a 50-yr amortization period, andthe allocation was made on 
a “present worth” basis, Table 5 contains a summary of the allocation. 

Even though no power generation at the site is involved in this project, the 
regulation provided by the project does benefit existing downstream power 
plants, and a cost allocation commensurate with the power benefits was made. 

The project has been operated for irrigation since April, 1954 under terms 
of an interim contract between the water users and the Corps of Engineers. 


TABLE 5,—SUMMARY OF THE COST ALLOCATION. 


Benefits 
Alternative cost 49,089 
Benefits limited by alternative 

costs 33,550 
Separable costs 1,503 
Remaining benefits 32,047 
Allocated joint cost 22,942 
Total allocation 24,445 


Allocation of first cost 


A permanent contract is being negotiated between the water users and the 
United States Bureau of Reclamation (USBR). Payments received so far for 
the interim irrigation service total $1,134,000, and flood damages prevented 
to date amount to $2,500,000. 


FOLSOM DAM AND RESERVOIR 


Folsom Dam was built by the Corps of Engineers. Folsom power plant and 
the Nimbus dam and power plant were built by the USBR. The entire project 
is being operated by the USBR as apart of the Central Valley Project. A pre- 
liminary cost allocation for Folsom Dam and Reservoir was made by the 
Corps of Engineers in 1953. The results of this allocation have been used by 
the USBR in various repayment analyses of the Central Valley Project. The 
preliminary cost allocation is given in Table 6, Note that the power benefits 
and the cost allocated to power are relatively smiall. The reason is that the 
power benefits were assigned to the power plants and to the dam and reser- 
voir in such a manner as to provide for an equal ratio of benefits to cost 
for the power plants and for the dam and reservoir. Only the power benefits — 
credited to the dam and reservoir are shown in Table 6. R 

An interesting sidelight in connection with this project is that no sooner. 
had the dam been completed when the greatest flood of record, that of — 
December, 1955 occurred, and the Sacramento metropolitan area was spared — 
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a great catastrophe. It has been estimated that during the 1955 flood, the dam 
prevented about $20,000,000 in flood damages, that ismore than the total cost 
allocated to flood control. 


SUCCESS DAM AND RESERVOIR 


A preliminary cost allocation for the Success Project was made in 1956. 
The percentage distribution of the project cost betweenflood control and irri- 


TABLE 6,—PRELIMINARY COST ALLOCATION, FOLSOM DAM AND RESERVOIR 


Flood 
control, 
in $1,000 
(2) 


Item 


Annual benefits 
Alternative annual costs 
Annual benefits limited by 
alternative costs 
Separable annual costs 
Remaining annual benefits 
Allocated joint annual costs 
Total annual costs 
M. O. & R. costs 
Interest and amortization 
First costs 


gation was fixed at that time. Such percentages will be applied to the actual 
construction cost when the project is completed. A summary of the cost 


allocation is given in Table 7(a). 
A contract is now being negotiated between the USBR and the water users 


for the amount to be repaid. The project is now under construction and is 
Scheduled for completion in 1962. 


TERMINUS DAM AND RESERVOIR 


A preliminary cost allocation for Terminus Dam and Reservoir was made 
in 1958. The percentage distribution of the project cost between irrigation and 
flood control was fixed at that time. These percentages will be applied later 
to the actual construction cost. A summary of the actual allocation is given 


in Table 7(b). 
A contract is being negotiated between the water users and the USBR for 


repayment of the part of the project cost allocated to irrigation. The project 
is under construction and is scheduled for completion in 1963. 


: NEW HOGAN DAM AND RESERVOIR 

: A preliminary allocation of cost for the New Hogan Dam and Reservoir 
vas made in 1959. This allocation fixed the percentage distribution of the pro- 
4 


: 


i 


j 
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ject cost between flood control and irrigation, but the amounts will not be 
finally determined until the construction costs are known. A summary of the 
allocation is given in Table 7(c). 

It is interesting to observe that construction funds were appropriated for 
fiscal year 1960, with the stipulation that construction not be initiated until a 
firm irrigation repayment contract had been executed. Sucha contract between 
the State of California and the USBR was executed in March, 1960, and actual 
construction was started soon thereafter, This contract was entered into with 


TABLE 7,—COST 


(a) Success Dam and Reservoir (b) Terminus Dam 


Flood Irrigation, Total, Flood Irrigation 
control, in in control, in $1,000 
in $1,000 $1,000 $1,000 | in $1,000 

55 
442 


Annual benefits 

Alternative annual costs 

Annual benefits limited 
by alternative costs 

Separable annual costs 

Remaining annual benefits 

Allocated joint annual 
costs 

Total M.O. and R, costs 

Interest and amortiza- 
tion 

First costs 

Percentage of total costs 


147 


the prime objective of permitting initiation of construction at an early date. 
The contract allows the USBR seven years tonegotiate directly with the water 
users for the repayment of the amount allocated to irrigation. If the Bureau 
does not arrange for complete repayment in that period, the State of California 
will make the remainder of the payments due under the contract. 


BLACK BUTTE DAM AND RESERVOIR 


The preliminary cost allocation for the Black Butte Dam and Reservoir 
was made in 1959. It fixed the percentage distribution of the project cost be- 
tween flood control and water conservation. Water conservation includes irri- 
gation and possibly municipal and industrial water supply. The fixed percent- 
ages will be applied later to the actual construction costs. A summary of the 
cost allocation is given in Table 7(d). 

As in the case of the New Hogan Project, funds appropriated for fiscal 
year 1960 could not be utilized until a firm repayment contract was executed, 
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Such a contract was executed in March, 1960 by the USBR and the State of 
California. The contract is identical to the one just described for the New 
Hogan Project. 

FUTURE PROJECTS 


In addition to the projects described heretofore, a cost allocation has just 
been completed for the Oroville Project of the State of California. The purpose 


ALLOCATIONS 


(c) New Hogan Dam and Reservoir 


and Reservoir 


(d) Black Butte Dam and Reservoir 


| Total, Flood | Irrigation, Total, Flood Water Total, 
: in control, in in control, | conservation, in 
| $1,000 in $1,000] $1,000 $1,000 in $1,000 in $1,000 | $1,000 


of such allocation was to determine the justifiable amount of Federal contribu- 
tion in the interest of flood control, sucha contribution having been authorized 
by Congress subject to the amount being fixed by the Secretary of the Army 
and approved by the President. A preliminary cost allocation has also been 
made for the proposed Hidden Reservoir on Fresno River. The reports con- 
taining these allocations are now being considered by the Chief of Engineers. 
If approved, the results of the allocations will soon be made public. 

Preliminary cost allocations are also being made in connection with pro- 
spective multiple-purpose storage projects on the Stanislaus, Tuolumne, Chow- 
chilla, Merced, Mokelumne, and Yuba Rivers. 

These possible future projects are noted simply to illustrate that the sub- 
ject of cost allocation in the Central Valley of California is important and 


bears on many current and future projects. 


APPRAISAL OF THE METHOD 


f A superficial examination of the separable costs- remaining benefits method 
may lead to the conclusion that there is finally available a specific and definite 


, 
i 
| 
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method that will produce the same results no matter who is making the cost 
allocation, Far from it. The method assumes that a great deal of data are 
available that were gathered during the project formulation stage. Survey 
report investigations, that include project formulation and economic feasibility 
studies, are carried on generally with limited funds and often short cuts are 
needed. Therefore, it is seldom that all the data required for cost allocation 
purposes are available. Even the data that were originally available have be- 
come partly obsolete due to the time elapsed since they were collected and 
due to the change in economic conditions. The data most frequently deficient or 
unavailable are the cost of the cheapest alternative projects and the separable 
costs. The development of these data at reasonable cost requires the exercise 
of mature judgment. The most difficult part involves the determination of the 
magnitude and location of the physical works required to obtain the specified 
services in each case, with reasonable assurance that the combinations 
selected are the most economical. Normally, this is a time-consuming trial 
and error job. Hence, short cuts by experienced men are almost a practical 
necessity. 

However, in spite of these drawbacks, the author believes that, as a general 
rule, the separable costs-remaining benefits method is by far superior to any 
other method with which he is familiar. It also has the advantage that it is 
beginning to be better understood and generally accepted. The allocations made 
by the Corps of Engineers in the Central Valley of California, by using this 
method, have been generally accepted as being fair and equitable by the water 
users and by the organizations involved. 
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VECTOR ASPECTS OF DYNAMIC SIMILARITY 


By R. C. Kolf,! M. ASCE, and W. L. Reitmeyer? 


SYNOPSIS 


A study of the laws of dynamic similarity emphasizing the refinements 
imposed by the vector nature of the fluid flow variables is presented, As 
examples of the method proposed, the problems of turbomachinery and geo- 
metrically distorted models that have needed separate treatment in the past, 
because of their similarity distortions, are shown as logical extensions of the 
same theory. 


INTRODUCTION 


The bases for the prediction of the characteristics of flow in a fluid me- 
chanics design problem are usually grouped in three general categories: (1) 
That of “engineering experience,” (2) the laboratory method of studying each 
problem by means of scale models, and (3) the process of theoretical analysis. 
In most of the problems that face hydraulic engineers, the shape of the physical 
boundaries that confine the flow are complex and they are often in motion. 
Consequently, the prediction of the flow characteristics by purely analytical 
means is many times impossible because of the lack of presently available 
theoretical methods that are adaptable to such situations. The second method 
of prediction from the study of hydraulic models, utilizing whenever possible 
the best features of the other methods, has proved very useful in such situ- 


ations, 
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of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961. 
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Notation.—The letter symbols adopted for use in this paper are defined 
where they first appear, in the illustrations or in the text, and are arranged 
alphabetically, for convenience of reference, in the Appendix. 


THEORY OF DYNAMIC SIMILARITY 


The laws of dynamic similarity are the criteria used in the design and 
analysis of hydraulic models. Depending on the type of model to be built and 
the characteristic to be studied, these criteria take on various degrees of 
sophistication. When only geometric similarity is of interest, the prediction 
equations can be written in terms of a single scale ratio, Ly, in which the 
subscripts m and p refer to model and prototype respectively. 


m 

ee ee puree iene s (1) 

r 

Mp 

A 

ae MA pe Dips ony Walaa 5 na mee (2) 
P 
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In Eqs. 1, 2, and 3 V is the volume, L refers to the characteristic length, 
and A denotes the cross-sectional area. 

A second class of variables involves time measurements as well as 
length measurements. For kinematic similarity a second scale ratio, namely 
the ratio of characteristic times, ty, must be established. 


m r 
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in which V denotes the velocity vector and a is the acceleration. 

A more complete listing of these prediction ratios is available in many 
standard reference works,? The significant characteristic to be noted at this 
point is that Eqs. 4 and 5 relate vector quantities, The division of one vector 
by another, as indicated by these vector ratios, has no known physical sig- 
nificance; hence the relationship so expressed is the ratio of the scalar or 
absolute magnitudes of the vectors. If this quantity is to have any pertinence 
in terms of similarity, it is necessary that each velocity or acceleration 
vector in the model have the same direction as the corresponding vector in 
the prototype. Then, in addition to Eqs. 4 and 5 being fulfilled, the flow pat-. 
terns must be identical. This condition is ordinarily fulfilled, or at least 
approximated, by insisting on geometric similarity as a prerequisite for 
kinematic similarity because, for each boundary geometry, the ideal flow 


3 ASCE Manual of Engineering Practice No. 25, Hydraulic Models, 1942. 
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wet is uniquely defined. Distortions of the flow net occur because of shear 
ind must be accounted for by the further conditions of dynamic similarity. 

Dynamic similarity requires that the ratios of allhomologous forces in the 
wo systems be the same. Limiting this study to cases in which only the 
rressure force, Fp the viscous force, Fy, the gravity force, Fg, and the 
nertia force, Fj, are significant yields 


Fim _ Fo)m _ elm _ Fim 


Fi), Lee : (g),, : (Fy) SE os «nine rnke eater (6) 


p 


ir, because the resultant force on a fluid element is called the inertia force 


a Pe Pies UO oat Fs ee Een Te (7) 
n 
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gain, insisting that the ratios indicated by Eqs. 6 and 8 are merely the ra- 
ios of the absolute magnitudes of the vectors involved, and that the homologous 
ectors on model and prototype must have the same directions in order to be 
imilar, it can be seen that Eq. 6 is a limiting solution of Eq. 8. Other solu- 
ions of Eq. 8 are conceivable, but to define similarity of the resultant (in- 
rtia) force without specifying similarity of the component forces would re- 
ult in a distortion of the geometry and the kinematics of the system in a 
1anner that has thus far been undefined. There is no present utility of such 
system. 

The symbolic relationships for the forces can be developed as follows: 


Fag eI Ee a RBI oe (9) 
F, =Ap Ac Ap bane i eT sg EK: (10) 
o tte abies air Ba betwee, eee, Bet (11) 
Fy May Avice Vit Tse; eR Sawn 
nd 
DF = Fj = Ma-pv v ia? We es (12) 


n which W is the weight, y refers to the specific weight, p indicates the 
ressure, denotes the dynamic viscosity, V is the velocity vector and p re- 
ers to the mass density. Eqs. 9, 10, 11, and 12 in the final form shown, all 
wolve a proportionality constant. For instance, in the case of the gravity 
rce, F, a y L’. In order that the force ratios in Eq. 6, have the same value 
egardless of the particular force studied, these proportionality constants 
hould be the same for the model as for the prototype. This will be the case 


: 
t 
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if the two are geometrically and kinematically similar. A hierarchy of simi- 
larity conditions is therefore as follows: 


1. Geometric similarity may be obtained exclusive of the other classes. 

2. In order to obtain complete kinematic similarity, geometric similarity 
must also be present, but in the ideal flow sense, dynamic similarity is not 
necessary. 

3. For complete dynamic similarity, both geometric and kinematic simi- 
larity are ncessary. 


The various dynamic similarity design criteria are fashioned by rear- 
ranging the terms in Eq. 6, 


(Fi) Fi) 2 
®. my “ a OF i, Fe. Se (13) 
cee * (Fi), el ae Rk (14) 
one ou arr Rigi s errs cree 

eA (Fi) (Fy) oor aah 
®). 2 ®), i A eC de (15) 


in which F is the Froude number, R refers to Reynold’s number and E indi- 
cates Euler’s number. ; 

Other combinations are, of course, possible, but they would not be inde- 
pendent of those shown in Eqs. 13, 14 and 15 and are, therefore, unnecessary. 
If a geometrically similar model is fashioned according to one of these equa- 
tions in such a way that it does not violate the other two, it will generally be 
dynamically similar to the prototype both in the sense of the absolute magni- 
tudes of the variables and in the sense of the direction of the vector quanti- 
ties. This is facilitated in many problems by the absence of one of these types 
of forces as a motivating factor, that obviates the necessity of conforming to 
all of the three equations (Eqs. 13, 14, and 15). Also, because the inertia 
force has been included in this study as the composite representation of all 
the other forces, one of the forces (and, therefore, one of the three equations) 
may be considered as being dependent onthe others. Thus, a particular model 


design will generally necessitate observation of only one or possibly two of 
these ratios, 


SPECIAL CASES 


The vector conditions implied in the development of these laws of simi- 
larity limit the sole use of the Froude, Reynolds, and Euler numbers to prob- 
lems in which the direction of the velocity and acceleration vectors is de- 
termined by the boundary geometry. In many problems this is not the case, 
Thus, designs that include varying degrees of curvilinearity, with or without 
the added complication of geometric distortion, have oftenbeen handled in the 
past as separate problems, and not as logical extensions of the same theory. 
The limits of the possible extension of the general laws have not been estab- 
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lished and, therefore, their effectiveness as tools of the hydraulic engineer 
has not been fully realized. 

An example of such a problem that has been successfully handled by ex- 
tending the laws of dynamic similarity is the vortex chamber. Essentially, 
this device consists of a circular or spiral shaped chamber with tangential 
and axial connections. It may be an enclosed (pressure force) system4 as in 
the Thoma counterflow brake or an open system (gravity force) as in the 
Portland combined sewer diversion.® If only the Euler number and Reynolds 
number conditions must be satisfied, a particular value of the one parameter 
will be accompanied by the same value of the other in each system, The 
following equations then apply: 


ia ee ieee 
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ot Geet D he ctie gh: Ail) Geniton aaa (18) 


in which Q is the volume rate of flow, C denotes the coefficient of discharge 
and H refers to the head. In the vortex chamber problem, this analysis was 
insufficient because, for a particular boundary geometry, there exist an in- 
finite number of possibilities for the direction of the velocity vectors de- 
pending on the strength of the vortex developed. This then is a case in which 
geometric similarity does not necessarily fix flow net similarity. In order 
to fix this similarity of vector directions Eq. 6 was amended to 
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This adds another necessary group to F, R, and E. Using the symbolic ex- 
pression for centrifugal force (that is a component of the inertia force), it 


becomes 
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and 
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in which m is the mass, u denotes the tangential component of the velocity 
vector and K is the Kolf vortex number, This expression relates the normal 
component of the inertia force vector (and, therefore, also the acceleration 
vector) to the tangential component. Along with the geometric similarity 
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conditions and the other dynamic similarity criteria F,R, and —E, complete 
similarity is achieved. In this problem, Eqs. 16 and 18 are modified to 


2 2 
V5- R= t (Rk) hd Sapien pt Rees 8. (22) 

and 8 Pp 
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Because viscous shear is significant only at lowvalues of R (that normally 
do_not occur in water flow problems), laboratory studies of the vortex cham- 
ber yield a curve establishing the dependence of the orifice coefficient, C, 
on K.6 


HYDRAULIC MACHINERY 


Another study that falls into this category is that of hydraulic machinery. 
Valid similarity laws have been developed for specific use in this field. How- 
ever, these criteria (such as specific speed, unit speed, unit discharge) have 
not been represented as amplications of the previously described methods, 
but rather as substitutions for them, thus restricting the view of the dynamic 
similarity concept. The following presentation, on the other hand, follows 
directly from the basic theory presentation. 

The pressure force and the viscous force should be considered in the case 
of a centrifugal pump, with the inertia force being again included in the study 
as the resultant of the two, If complete dynamic similarity is assumed, it 
should be possible to show the Euler number as a function of the Reynolds 
number, but this simplification of the problem would not prove useful. Be- 
cause of the variability of the angular velocity of the impeller, kinematic 
similarity in the vector sense cannot be fixed by geometric similarity alone. 


The “speed ratio,” y that is the ratio of a characteristic velocity, V, to the 


tangential component of a characteristic velocity vector within the impeller, 
u, has been used to fix this kinematic similarity. This ratio is recognized to 
be the reciprocal of K from the vortex chamber consideration. It is, there- 
fore, proportional to the square root of the ratio of inertia force to centrifugal 
force. Using the relationship u « N D and V « Q D-2 


GA. Wee Ce 
Un by D3 
in which N is the angular velocity and D denotes the diameter. Most of the 
problems of engineering interest involve high Reynolds numbers, The pres- 
sure force is then the predominant force, and in the same way that Eq. 16 
was modified to the form of Eq. 22 for the vortex chamber, the Euler number 
must be modified by Eq. 24 or an equivalent expression, é 


~ 
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Fig. 1 shows a plot of Eq. 26 for a vertical end-suction pump operated as a 
laboratory model, with an impeller diameter of 20-5 /16 in,, and a geometri- 
Cally similar prototype with an impeller diameter of 50 in. The speed of the 
model was 1267.5 rpm and the prototype speed was 514 rpm. The character- 
istic velocity, V, in Eq. 26 was arbitrarily chosen as the average velocity in 


the exit pipe. 


VERTICAL END 
SUCTION PUMP 


75% EFF. 


Q 


ND3 x102 


BIG. 1 


To the scale shown it was not possible to differentiate the data for the two 
pumps, the “spread” being less than the thickness of the line. Because of Aes 
dynamic similarity defined by the curve, each point represents a gts: 
efficiency for this design. Curves of this nature have also been used w: e 
Same success to represent the data for a single pump operated at various 


speeds. 

From Eq. 20 
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It is of interest to note that 
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It will be recognized that the Eqs. 29 through 34 describe the “laws of 
homologous operation” of turbo machines. 


GEOMETRIC DISTORTION OF OPEN CHANNEL MODELS 


The problems involved in applying the laws of similarity to the vortex 
chamber and to turbomachinery might be classified as resulting from a 
“kinematic distortion,” That is, true vector similarity is not assured by geo- 
metric similarity. The situation encountered when geometric distortion is 
necessary, although falling in a separate category, has many of the same 
characteristics. 

Geometric distortion is often a necessary evil in river hydraulics. The 
Reynolds and Froude criteria must be observed inthis case. Although viscous 
shear is not a predominant force except at low values of Reynolds number, 
this restriction places a lower limit on the size of the model in order that the 
flow will be well into the turbulent range, and the modifications of the velocity 
profiles due to shear will be like those of the prototype. This lower limit on 
the size is, in many instances, unreasonable for economic reasons. When this 
is the case, geometrically distorted models are used in the hope that, although 
complete similarity is impossible, Similarity with respect to an isolated — 
phenomenon may be achieved. If the phenomenon to be studied involves a 
“movable bed,” the conditions of sediment transportation must also be con- 
sidered. 

Hans Einstein, F. ASCE and Ning Chien” have given an orderly method for ~ 
the design of such models, Several severe discrepancies have been noted, — 
however, using these presently accepted techniques. This may be illustrated 
by considering the Froude number in which the characteristic length is con- ; 
sidered to be the depth. The reason for this is that accelerations of the fluid — 
are accompanied by changes in the elevation of the water surface. The justi- ‘= 
fication is evidently an assumed two-dimensionality of the flow with the velo- 
city variations in vertical, longitudinal planes more important to the problem 


ee SSSSSSSSSSSSSSSSSe 
7 “Similarity of Distorted River Models,” by H. A. Einstein, N. Chien, Transactions, — 
ASCE, Vol. 121, 1956. =a 
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than the transverse variations, This basic assumption can be shown to be a 
common source of error. 

All ratios of homologous accelerations should have the same value if kine- 
matic similarity is to be achieved. 


ay = Sr = 1 ie, 910. 6: 16\ 6.10 '& © we 61 © 6b w: 10) 01m 6 (35) 


Let us suppose that the grain size of the moveable bed for the model is 
considered along with the velocity of the stream necessary to satisfy the 
Froude conditions based on depth, so that approximately similar sediment 
transport conditions are produced. The most common type of geometric dis- 
tortion is to use a depth ratio, d,;, that is greater than the length ratio, Ly. 
Good engineering practice requires that the model and scales chosen be veri- 
fied by reproducing known prototype developments with similar flows. Mis- 
cellaneous effects that have been observed in such model studies include8 
(1) excessive scour at bends and the beginning and end of obstructions such 
as deflecting dikes, (2) improper location of scour, (3) excessive bar building, 
and (4) formation of eddies in which none occurred in the prototype. All of 
these difficulties occur in regions of significantly curvilinear flow. The trans- 
verse velocity distribution in such regions (scour) is related to the normal 
force experienced by the fluid. The transverse currents (bar building) can 
also be linked to this phenomenon. If the modei has been designed to give 
similarity of inertial conditions in the straight portions of the channel, the 
distorted normal acceleration will result in dissimilar conditions in the curved 
sections. Because the acceleration ratio is fixed, Eq. 35, another criteria 
that should be observed is 

2 
Vy 
a, = R, =n SARE Cas ch eae pe oe (36) 


Supposing that the depth ratio, d,, is greater than the length ratio, Ly, by a 
factor of 6, that is not uncommon, and that the Froude condition has been ob- 


served, then 
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It should be noted that, because Ry = Ly, both of these equations cannot be 
satisfied simultaneously. Eq. 36 will have been violated by the same factor 
6, resulting in severe distortion of the transverse conditions and greatly in- 
creasing the effects of the normal acceleration. Several possibilities of 
alleviating this condition exist, but there is no published record of their at- 
tempted use. One method would be to model the radius of curvature of severe 
bends to the depth ratio instead of the length ratio. This would amount to 
partially eliminating the distortion in the zones that normally cause difficulty. 
This type of distortion would introduce a psychological problem to the observer 
because it would not “look” like the prototype condition. It would not eliminate 


the need of the verification step. 


8 «Observed Effects of Geometric Distortion in Hydraulic Models,” by K. D. Nichols, 
Transactions, ASCE, Vol. 104, 1939. 
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In rivers that normallyflow at relatively large depths (small Froude num- 
bers) it seems that the Froude number as computed from the depth loses its 
significance. In such channels, small accelerations of the fluid do not greatly 
affect the depth. Here then, it would seem better to evaluate the parameter 
using the length ratio, thus satisfying Eq. 36 in preference to Eq. 37. 


CONCLUSIONS 


1. The limitations imposed by the vector nature ofthe fluid flow variables 
necessitates reformulation of the statement of the laws of dynamic similarity 
in order to extend their utility into new realms of application. In some cases 
this extension may be accomplished by using the “Kolf vortex number” in 
addition to the classical similarity criteria. 

2. A large class of problems in which variable amounts of “vorticity” are 
imposed by moving portions of the flow boundaries (such as pumps, turbines, 
and so on) may be included in the general study of similarity by adding this 
parameter (ratio of centrifugal to inertia force) that assumes like directions 
for homologous vectors. 

3. Many of the discrepancies noted in geometrically distorted models 
arise from a violation of the normal acceleration requirements in zones of 
curvilinear flow. Recognition of this fact indicates possible model adjustments 
for compensation. 


ACKNOWLEDGMENTS 


The authors are grateful to the Marquette University Department of Theo- 
retical and Applied Mechanics that, through the interests of W. G. Murphy, 
Chairman, gave full support to this study. 

Financial support was given by the National Science Foundation, Engineer- 
ing Sciences Program, 

Data on geometrically similar pump operation wasfurnishedby Fairbanks, 
Morse, and Co, 


ADDITIONAL READING REFERENCE 


1, Engineering Hydraulics, by H. Rouse, John Wiley and Sons, New York 
1950. 


’ 


‘ 
7 
; 
y 
3 
% 
3 
7 
% 
3 
3 


APPENDIX.—NOTATION 


The following symbols adopted for use in this paper and for the guidance 
of discussers, conform essentially with “American Standard Letter Symbol 
for Hydraulics” (ASA Y10.2-1958) prepared by a committee of the Americ 


HY 2 


DYNAMIC SIMILARITY 
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Standards Association with Society representation, and approved by the As- 
sociation in 1958: 


Peeraye ce 4 


a 
A 
Cc 


D 
E 


d 
F. 
Fy 
Fi 


F 
p 

F 
Vv 

F 


g 
H 
K 
L 
m 
N 
p 
P 
Q 
R 
R 
s 
t 
u 
V 
Vv 
WwW 
tf 
m 
p 


= acceleration 

= cross-sectional area 

= coefficient of discharge 
= depth 

= diameter 

= Euler number 

= centrifugal force 

= gravity force 

= inertia force 

= pressure force 

= viscous force 

= Froude number 

= acceleration due to gravity 
= head 

= Kolf vortex number 

= characteristic length 

= mass 


= angular velocity 


pressure 
= power 

= volume rate of flow 

= radius of curvature 

= Reynolds number 

= distance along streamline 

= time 

= tangential component of velocity vector 
= velocity vector 

= volume 

= weight 

= specific weight 

= dynamic viscosity 


= mass density 
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THEORY OF WAVE AGITATION IN A HARBOR 


By Bernard Le Méhautél 


SYNOPSIS 


A first order theory for the value of the wave agitation in a simple basin 
caused by an incident periodical gravity wave is presented. This study is a 
natural result of experimental results and previous theory. 

A harbor may be considered as a combination of discontinuities separated 
by a definite distance. The ultimate one, often in the form of a beach, is a 
total obstruction that may dissipate a part of the incident energy. The value 
of the agitation in such a basin is computed when the motion is predominantly 
two dimensional, as afunction of the amplitude and period of the incident wave. 
A practical example applies the results of the theory directly. 

The influence of the shape of the entrance, beach, wave traps, and so on 
on the value of the agitation are analyzed. The obtained formula for the magni- 
tude of agitation emphasizes that the ratio of the length of the basin to the 


Note.—Discussion open until August 1, 1961. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE, This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 


ciety of Civil Engineers, Vol, 87, No. HY 2, March, 1961. ; , 
1 Ingenieur-Docteur, Special Lecturer in Graduate Studies and Research Associite 


at Queen’s Univ., Kingston, Ontario, Canada. 
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wave length is important because of multiple reflections, Resonance condi- 
tions are computed. 


Results obtained by theory are compared with those obtained experimentally. 


INTRODUCTION 


This theory of agitation in a harbor could be considered as a natural out- 
come of the experimental results presented2 in 1954 and3 eras and of the 
author’s theory of periodical gravity waves on a discontinuity .4 


FIG, 1.—LAYOUT OF A BASIN 


A harbor is considered hereinas a combination of discontinuities separated 
by distances that define the lengths of the basins. The terminal discontinuity 
is, of course, a total obstruction. This total obstruction, frequently formed 0 
a beach, dissipates a noticeable part of the incident wave energy. 

With this assumption, the value of the agitation in a harbor may be com- 
puted from a general theory, that considers the behavior of waves acting on a 
combination of a number of discontinuities. This broad theory may be applied 


2 “Two Dimensional Seiche in a Basin Subjected to Incident Waves,” by Bernard Le 
Méhauteé, Proceedings, Coastal Engineering, 5th Conference, 1954, s 
“Mouvements de Resonance & Deux Dimensions dans une Enceinte sous l’action 
d’Ondes Incidentes,” by F, Biesel and B, Le Mehaute, La Houille Blanche, July, August, 
1956, No, 3, pp. 348-374, 


4 Periodical Wave on a Discontinuity,” by B. Le Méhauté, Journal, ASCE, Vol. 86 
No. 9, September, 1960, 
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to a number of such phenomena as waves passing over a shelf, agitation be- 
tween two barriers and wave resonators, agitation in a succession of basins, 
and so on. 

However, because of its importance in engineering practice, and because 
of the experimental data now at hand, only the case of a relatively simple 
basin, as shown by Fig. 1, is considered herein. It is, of course, a particular 
case in the more general theory of waves on a combination of discontinuities. 

As in the author’s previous work,4 great use is made of the complex 
number method calculus. 

It will be noticed that this method of complex number calculus provides a 
Simple and powerful tool for the analysis of gravity wave motion. Solutions 
will be obtained by this method with a minimum of computation for cases 
which would usually be solved by finding a complicated potential function, @. 
The method of the potential function, ¢, is always tedious if there are disconti- 
huities in the boundary conditions. 

An example is presented which is consistent with the theory previously 
presented. 

The theory is limited in that only termsof the first order of approximation 
are considered. That is, it deals only with the linear solutions in which, for 
example, the convective inertia effect, that is proportional to the square of 
the velocity, is neglected. This particular assumption seems valid for the 
majority of cases. Also, the theory deals only with the main vertical and 
horizontal components of sinusoidal motion. This stipulation tends to make 
the theory invalid as the discontinuity is approached, but it seems sensibly 
valid at distances of the order of two or three times the depth from the dis- 
continuity. Three dimensional components, if existing, must have a limited 
effect. They are theoretically negligible when the width of the basin is less 
than half the length of the wave inthe same domain. However, even for greater 
width, experiments indicate that this theory provides a fairly good approxi- 
mation for the value of agitation, as long as no transverse resonance occurs. 
Transverse resonance is likely to occur when the width of the basin is very 
close to a multiple of half wave lengths.4 

A linear periodical gravity wave may be definedby a rotating vector whose 
modulus is equal to half the wave height and whose argument is equal to the 
phase of the free surface at agiventime and at a given point. In the same way 
a clapotis is defined by two vectors rotating in opposite directions. 

Notation.—The letter symbols adopted for use in this paper are defined 
where they first appear, in the illustrations or in the text, and are arranged 
alphabetically, for convenience of reference, in the Appendix. 


COMPLEX VALUE OF THE AGITATION IN A BASIN 


Consider a basin aligned in the direction of the incident wave as shown by 


Fig. 2. 


The incident wave, coming from the open sea arrives at the mouth of this 


; n 
basin with a half amplitude Aj and a phase Ay atagiventime. _ 
_ This wave is considered to be defined by the complex number A, of modu- 


lus | Ay | = Ay, and argument Ne Similarly, the wave reflected by the over- 


tT 


. J 
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all layout (entrance and basin) is defined by a complex number By of modulus 


a NN 
| B,| = By, and argument By. If there is no loss of energy at the disconti- 
nuity and inside the harbor 


De Aepchesiad SB py lace eM « leaping ae (1) 


The wave Aj is partially transmitted and partially reflected. A 
The transmitted wave is a; A, of half amplitude oa, A, and phase (@y 


+ Ay) and the reflected wave is Py Ay of semi-amplitude #1; Aj, and phase 


(By + Ay). The terms a1 and By refer to the coefficients of transmission and 
reflection, respectively. They are defined by complex numbers of moduli aj 


IN 
and 6; and arguments a and £1, respectively. 
These moduli (a4 and #1) are equal to the ratios of the transmitted and re- 
flected wave heights to the incident wave height, respectively. Similarly, 


WRN 
‘ 
oI 


| 


FIG, 2.—NOTATION OF PRIMARY WAVES 


their arguments ay and a are equal to the change of phases, caused by the 
discontinuity, of the transmitted and reflected waves, respectively. 

Now, consider a wave defined by a complex number Dy arriving from the 
basin at the discontinuity. A part Bo Dy is reflected within the basin, and 
another part ag Dy is transmitted to the open sea. The terms a9 and Bo are 
two complex numbers defining the transmission coefficient towards the open 
sea and the reflection coefficient towards the basin. 


The four complex Saat a A, a9, ) or eight parameters moduli 4, 
84, 29, By and arguments Or A, Qo, Bo completely define the flow character- 
istics of the entrance of the basin. 

Now, consider the four translatory primary waves Ay, By in the open sea 


and Cy, Dy inside the basin at the entrance to the harbor. Because all thea 
waves have the same period 


Cy = ay Ay + By Dy ah uae 


and fo 
By = aq Dy + By Ay PRD. eae te es Teh ee ee 
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The wave C1, moving inside the basin away from the entrance, arrives at the 


end of the basin, with a different phase A This phase is a function of the rela- 
tive length of the basin with respect to the wave length in the basin 


Ne 271 
r= 


inwhich | is the length of the basin and Lo the length of the wave in the basin. 


The function D is negative because there is a negative phase difference be- 
tween the wave at the entrance and the same wave arriving at the back. The 
wave height remains constant if friction is neglected. However, if friction is 
included, as produced for example by a wave trap or by lateral beaches, then 
the amplitude of the wave is reduced in the ratio r(z <1). Asa first degree 
of approximation, r may be written 


in which K is a friction coefficient. 
Hence, the incident wave at the back of the basin becomes 


ai Raritan as aay Ree ts) 
in which r is a complex number of modulus r, and argument %: 
TN 
Tae 


in whichi=v-1. a 
In the case of a single basin, the obstruction is complete, and the wave Co 
is reflected without change of phase. If Do is this reflected wave, then 


1 ays 
Cet 3) ce peer ee rs eka aN 


If reflection is assumed to occur without any loss of energy, then 


Co = Do Cp ohh Mok SA Eee TSs) 


and ~  o = 
Gye Dg lacie) ++ no eben 


However, the loss of energy, such as, that due to the breaking of a wave 
incident on a terminal beach, must often be considered. Hence, it is more 
accurate to write a a5 
Dy, = Py COi oie serisientgs csi tie sie (9a) 


a Dg = PC Mee Pate eo ae ee OB) 


in which p is a coefficient of reflection ( eS 1) that depends on the steepness of 
the incident wave, the slope, the roughness, perviousness, and the length of 
the beach. The coefficient p may be computed by use of theory 5 or be obtained 
experimentally as in the chart given by Messrs. Greslov and Mahe.” Asa 


onstitue par un 


ici s ion d’une houle’ sur un obstacle c 
5 “Etude du coefficient de reflexion nie 


plan incliné,” by Greslou and Mahé, Proceedings, Coastal Engineering, 


1954, 
6 «Le pouvoir reflechissant des ouvrages maritimes exposes a l’action de la houle,” 


by M. Miche, Annales des Ponts et Chaussées, Mai, Juin, 1951. 
: azcE—aE—E—oeoecre’rerv—o? 
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first approximation, p may be considered as a function of the ratio of the 
length of the beach to the width of the basin, assuming that all the wave energy 
arriving at a beach is dissipated in the breaking zone. 

However for a long wave, such as a seiche type wave for which the period 
is of the order of 100 sec the reflection on a beach is almost total even when 
it has a very gentle slope so that we have Eq. 8b 


Up: = Dgtves 20'S th. See bas (8b) 


The wave Do moving towards the entrance of the harbor becomes D1, such 
that 


Eqs. 1, 2, 5, 9a, and 10 form a system of five equalities that permit the 
computation of the five unknowns Cj, Co, Do, Dy and B, as functions of the 


incident waves Ay and the coefficients Q4, Qo, By (Bo i)is each a complex 


number of modulus | i! =1 and argument = 0, that is the real number 1. 
(It should be noted that iis, in this case, different from the imaginary num- 
ber i= V- 1 used in Eq. 6.) 

The following formulas are obtained from the previous equations: 


a ay = 
Cy = eg iae Ay Sle ls ete(n Wales ole eae (11) 
i-p r2 Bo 
kt eer 
es ie eee eS: (12) 
= pt ~ Bs 
= pra am 
Do = = ane wie EE a Tk (13) 
i tp Po 
p ae Q 
a 1 a 
D1 = = oe Ay id re ee ee (14) 
and Sat RR) 
pF? a a \ _ 
By = By + Pee Ay + Ye: Ghar ar ethep wate eee (15) 
i-pr Bg 
At the entrance the total value of the agitation in the basin is given by the © 
equation 
2 a(i +pr ai 
2(¢, +D,)= = Go) = Ay ra ee ee, we ee (16) 
d = PT Bg 
and at the back by the equation 
2T ay (1 - p) ; 
Con + Do):® acer yc Ay a(t Wiel evel svter eoret > (17) 
ae Oe ae =))) 
The value in the open sea, at the mouth, is 
at = an ee pare a ao pS 
2(Ay + By) = 2(7 + By +——_ Ay Boe) ole otters: (18) 
i =p FP - Bo 
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The difference of phase between (Ay a Bj) and (Cy + Dy), that is argument 
Ay+By 
=—— > Sives the difference of phase of the agitation bétween the two sides 
Cy + Dy | 
of the entrance and enables an estimate to be made of the current at the en- 
trance. 
Verification.—A simple verification of these formulas consists in comput- 
ing the absolute value of By when the loss of energy in the overall layout is 


zero, that is when p = 1 and r = 1. In this case, By should be equal to Ay. 
This condition is obtained if 


pay Fe pels 
PT 2 
By + 58 ee oe A sigs se Wed et sis at es care (19) 
LSE oy ERE S/S) 
or 
bs Syn = aD at Pd 
(i = By) (i =" Og Bo)= PY Gy a so ee sass (20) 


Introducing Eqs. 19 and 20 and utilizing the relationships previously pre- 
sented 


A a + By =i Ser ne ats ee es ee (21a) 
and Ps 

ag —_ = 

ror By RPL Se) eaieMiat abet cs Vanes ate aes Tome (21b) 


it is found that Eq. 20 is verified; that is Ay = By. 

Moreover the value of the agitation in the basin at the back of the harbor, 
when it is assumed that all the incident energy is destroyed by a terminal 
beach such that p = 0, is found to be 


21 Cot+Dy!=21 Co! =29 a Ay ....--. ( 22a) 


At the entrance ung: ‘y 

21 Cy +D, 1 =21 Cy 1 =2 ay Ay......--. ( 22) 
that is also consistent, because these values represent the transmitted wave 
from the open sea into the basin. 


ABSOLUTE VALUE OF THE AGITATION 


The value of the agitation at the rear of the harbor is given by the absolute 
value 2 | Co + Do |; at the entrance in the basin by the absolute value 21 Cy 
+D | ; and at the entrance inthe open sea bythe absolute value 2 | Ay + By is 
The amplification factor is given by dividing these values by the amplitude of 


the incident wave, Aj. 
At the back of the harbor, the factor is easily found from Eq. 17 to be 


Ag 2° 'CG # Do! 2r (1+ p) ay ae 
Ay” " 3 ae SES 
AY Ay [1 +(p r2 Bo)” - 2pr2 Bo cos (27+ »)| | 
% 
From previous theory4 SNE 

B= (F< AZ al) 2--++eee nese . (24) 
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ad eae EC Ss 
COS By = Dia cc (25) 
2(1- A Za;) 
in which L. 1 
2 ‘2 
Cerne 'a Ww wink aoe See ae (26) 
Ly 


The terms lo and ly define the widths of the domains, 2, the open sea, and 1, 
the basin, respectively and Ly and Ly are the wave lengths in these two do- 
mains. The term L is a function of the period T and the depth h: 


gt 2a7h 


In addition (m = =) 2 mo ho 


sinh 2 my hy 
Hence, it is possible to express the value of the agitation as a function of Z, 
A, and qj. 

The term ay may be expressed as a function of the shape of the disconti- 
nuity by the empirical formula 


1/2 /1,\1/4 
b i 2 
pn) Tn Se Ca De Aes eel See (29) | 
“t (i) (2) Ly 


cesta PE 


in which b is the width of the entrance of the basin. 

Finally, Eqs. 3, 23, 24, 25, 26, 28 and 29 permit the computation of the val- 
ue of the agitation at the rear of the basin for a given incident wave of ampli- 
tude 2 Aj and period T. 


PRACTICAL FORMULAS 


The value of the agitation at the back of the basin is given by the theoreti- © 
cal formulas 
2r(1+p) ay Ay 


ar eet kar Obie San he a (23) 
re [1 +(pr? Ba) - 2p r2 a cos (2+ 4)/* 2 ’ 


in which r is the damping coefficient of the amplitude of the wave travelling ‘ 
from the entrance to the back of the basin; p is the coefficient of reflexion of © 
the back of the harbor (p is small when this back is a beach); ay is the coef- 
ficient of transmission for the wave arriving from the open sea into the basin; 
and is computed from Eq. 29; &9 is the coefficient of reflexion of the entrance 
for the wave travelling inside the Lag <3 the direction of the open sea. 89 is” 
given as a function of a, by Eq. 24; Tis the difference in phase in the basin - 
between a wave at the entrance and a wave at the back of the harbor. Dis 


ws 


eet. ee 
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given by Eq. 3; By is the difference in phase between the wave reflected by 
the entrance inside the basin inthe direction of the basin and the wave travel- 
ling in the direction of the open sea also inside the basin; 8b is zero in the 
case of a total obstruction and equal to 7 in case of a total opening in a wide 
and deep open sea; and is given by Eq. 25. 

In these formulas, use is made of coefficients: the shoaling coefficient 
A(m = 2 7/L) as defined by Eq. 28 and the change of domain coefficient Z 
caused by any change from the open sea to the basin, and defined by Eq. 26. 

Example.-—-Consider a basin as shown by Fig. 2 in which the length of the 
basin 1 = 1450 ft; the width of the basin lg = 500 ft; the length of the mouth b 
= 187 ft; and the depth of the basin ho = 32 ft. A terminal beach exists with a 
very gentle slope and a length of 340 ft. The open sea (or another basin) is 
66 ft deep (hy) and 625 ft wide (11). The incident wave has a period of 10 sec 
and a height of 2 Aj = 15 ft. 

Preliminary Computation.—A 10 sec period wave in the basin in a depth 
hg = 32 ft, has a wave-length Lg = 300 ft. A 10 sec period wave in the open 
sea in a depth hy = 66 ft has a length: Ly, = 400 ft. Hence 

Lo 


oe = 0.75 


L L 
Because lo = 375 ft > 3. and ly = 750 ft > $ the theory is not strictly ap- 


plicable but should provide a fairly good approximation to the value of the agi- 
tation. A transverse resonance is unlikely because lg and 1, are not close to 


L2 li 
any multiple of > and > respectively. 


Computation of the shoaling factor A by means of Eq. 28 proceeds as fol- 
lows: 


4 _ (47m) (32) _ 
NY eS ho = 399 = 1-34 
sinh 2 mg hg = 1.78 

2, Tito neha (66) x: 
ABS 5 a re SAN 


and 
sinh 2 my hy = 3.90 


Substituting into Eq. 28 yields 


= 1.14 
Substituting 


into Eq. 26 yield the change of domain coefficient 
Z = 0.75 x 0.80 = 0.60 
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Deriving the values 
é 1/4 


l 
(71) = (1.25) 2/4 = 1.06 
lg 


py/? (181)? « o.61 
Ip 500 ‘ 
2 . 2 " 
Lo 1+ 1.14X0.75 


As Mg Pe 


and 
1.08 


and substituting into Eq. 29 yields the transmission coefficient 
a4 = 1.06 x 0.61 x 1.08 = 0.70 


From Eq. 24 the reflection coefficient is found to be 


91/2 
By = [1 - 1.14 x 0.60 x (0.70) | = 0.82 
Substituting into Eq. 25 
A _ 2- (0.70)? (1,14 + 0.60) _ 


ana 2x 0.82 a8 
the phase of the entrance is found to be 
By = 45° 


From Eq. 3 the phase lag of the basin is found to be 
af. 9 21 10x d+ 120° 
Lo 


Bonee. cos (24+ 8» ) = cos (120° + 45°) = - 0.97 

It is assumed that the friction effect inside the basin is given by putting r 
= 0.95. 

The terminal beach reflection coefficient p is then computed. When the 
beach length d is small relative to the width of the basin lo, the coefficient of 
reflection of the back of this basin can be computed by assuming that the wave 
energy dissipated is proportional to the beach length so that: 


lo 
When 2 is of the same order as lg, it is more accurate to compute the reflec- 
tion coefficient as a function of the beach length, its slope and the steepness 


2C 
of the incident wave re In this example, it is assumed that p = 0.1. 


Now the value of the amplification factor may be computed by the use of © 
the theoretical Eq. 23: 5 


oy 0.95 x (1 + 0.1) x 0.70 

a [1 + (0.10 x 0,90 x 0.82)” - 2(0,10 x 0.90 x 0.82) (- 0,97) | 1/ 

and the value of the amplitude of the agitation at the back of the harbor is 
Ag = 0.62 x 2 Ay = 0.62 x 15 = 9.3 ft 


=1.24 
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. . . . L 
It is likely that this value is slightly over-estimated. Because lo > = 

5 Saree 
and 1; > > it is apparent that some secondary three dimensional motions 


are Superimposed on the primary two dimensional motion. 
The Importance of the Ratio of the Length of the Basin to the Wave Length.— 
It should be noted that a slight variation in the term (2%) in the expression 


2p r2 Bg cos (2 “ + 2) usually has a large effect on the value of the agitation. 

If 1 is of the same order as Lg,a slight variation of T, does not change the 
value of the cosine given by the previous formula nor, consequently, the agi- 
tation very much. In practice, this is the case for long waves such as tides 
in a bay or seiches in a basin. But as far as ordinary wind generated waves 
and swell are concerned, | is generally much greater than Lg, and, a very 
small change of the period of the incident wave T, causes a change in the value 


of a that introduces a great variation in the value of cos (2 S + 22) and conse- 
quently in the value of the agitation. 

Because waves in the field are always irregular it is necessary to consider 
the “degree of agitation” that could be found by computing the mean value of 


the agitation for a complete cycle of (27+ By ), or more simply of 2% when 
the wave length, that is the wave period, varies. 

It is interesting to note that the “degree of agitation” does not depend on 
Bo, whereas the value of the agitation for a given period is very much influ- 


A 
enced by the value of Ao. 


RESONANCE 


Maximum agitation occurs when: 


| T-pr? Bp =1-pr2 py aA, GU Bes (31) 
a A ee fai inh ac oP a ul ge (32) 
that gives 
DL ear hilar EE gig ER (33a) 
n> oT 3 ertert rt rie oe 
q 5 28 ee (33b) 
2nm7+ Bo 


If the water is shallow such that L = Tvgh _ , then the periods of resonance 


are given by 
Wyeth Reavy (ROR int gies emg.) 


Vg ho regs 


The corresponding maximum amplitude of agitation is equal to 


it 


k 2r(1+p) ay 
*Gmax 1 - pr” Bp 


the values of p and r both tend to re- 
tuitively. 


Ay - eee es case USD) 


i 
the effect of friction is apparent where 
free the value of the amplitude of resonance, as would be expected in 
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Minimum agitation occurs when 


| repre ee = teppei BR (36) 
That is, when 

27+ By = 1(+ 2nT) See Per peat pte (37) 
en Lig ges lin Pat noe 

an yy eee PR RS SEN i ae ae ee 8 

as Siar anna aa, enna ar Rape (38a) 
and 

Loaztl er (38b) 

2na7+7+ Bo 


The minimum amplitude of agitation is equal to 


2r(1+p) a 
A TR 
Smin 1+pr Bo 


EXPERIMENTAL VERIFICATION 


The theory presented herein has been verified by a number of experiments 
conducted in the Sogreah Laboratory at Grenoble, France. Details of this ex- 
perimental work have been presented previously. 2,3 For each case experi- 
mentally analyzed, the following results were obtained: 


1. Thecurves ofthe coefficient of amplification = are given as functions 


of the relative length of the basin for various shapes of the entrance and vari- 
ous characteristics of domain 1, the open sea, and domain 2, the basin. That 
is they are given as functions of the opening of the basin and the depths and 
widths of these two domains. The theoretical curves are computed assuming 
no loss of energy in the basin, that is 


No terminal beach; p = 1 
No damping of the wave: r= 1 


These theoretical curves are compared with the experimental results in Figs. 
3, 4, 6, and 7. 

2. The maximum value of the agitation that occurs, of course, at reso- 
nance has been analyzed and presented in Fig. 5 for the case of a basin closed 
by a simple obstruction. The computed values presented assume 

a) no friction 
and_ bb) a constant value of the friction terms p and r. 


These two theoretical curves are compared with the experiment. ; 
Successive computations have been made for the following harbor entrances: — 


1. Obstruction in a flume. (The domains 1 and 2 have the same charac- > 
teristic.) (Figs. 3, 4 and 5) 

2. Change of width, constant depth, and full opening. (Fig. 6) 

3. Change of depth,. constant width, and full opening. (Fig. 7) 


It has also been verified that this theory could be applied to a number of 
other combinations with the same degree of accuracy. 
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FIG, 7.—AGITATION IN A BASIN WITH A CHANGE OF DEPTH 
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The experimental and theoretical data presented in Figs. 3, 4, 5, and 6 are 
considered to exhibit fairly good agreement. 

In the case of a basin closed by an obstruction, the friction has, as would 
seem logical, an important effect particularly on the maximum value of the 
agitation that occurs at resonance. Here it is necessary to take into account 
the friction coefficients r and p, however their values are usually difficult to 


determine accurately. The theory, particularly concerning the values of - 
2 


causing the resonance, is well verified. 

For the case of a simple change of width, the friction effect may be neg- 
lected, (p = 1, r = 1) except when a terminal beach exists: p<1. The com- 
puted values of the maximum (resonance) agitation is particularly well veri- 
fied by the experiments, but there is a slight discrepancy between the theo- 


A 
retical values of 89 causing the resonance and those obtained by experiment. 
A 
This means that the value of cos £9 calculated by introducing ony ts (14/12) 1/4 
A 
is incorrect as far as the value of By is concerned. As previously seen, the 


degree of agitation caused by an irregular wave does not depend on Os Hence, 
the computation of the degree of agitation caused by an irregular wave is not 
affected by this discrepancy. 

In the case of a simple change of depth, the occurrence of resonance as a 
function of 1/Lo is well verified. However, the magnitude of this maximum 
obtained experimentally seems greater than that computed from the theory. 


SUMMARY 


The results presented and cited indicate that the theory is generally well 
verified by experiments and canbe successfully applied to compute the period 
and the amplitude of resonance, particularly if friction effects are small. 
Above all, the degree of selectivity of resonance given by the theory is always 
in accordance with the experimental results. This is shown by the relative 
thickness of the peaks of the agitation curves. 

The following is found both by theory and experiment: 


1. That the occurrence of resonance increases with a change of depth and 
a change of width, and decreases with an obstruction such that it tends to zero 
when this obstruction tends to be total; 

2. That the value of the agitation at resonance increases as ly /lp, Ly /Lo 
or hy/hg and lo/b increase. It will be noticed that the value of agitation tends 
to infinity with lo /b, that is in the case of total obstruction, when the friction 
is neglected, but at the same time the occurrence of resonance tends to zero. 
However, when friction effects are taken into account, the value.ot the agita- 
tion tends to zero when lg/b tends to infinity. ee 

3. That resonance occurs when 


Lo 1\ Le 


for all cases; 


aa 
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in the case of a complete obstruction; 


in the case of a change of depth without obstruction; 
4. That the minimum agitation occurs when: 


L 
(v2) Bcr<ien 


5. That the agitation is proportional to the amplitude of the incident wave, 
even inthe case of viscous friction as long as the quadratic friction is consid- 
ered negligible. 


It is relatively easy to obtain approximately the reflection coefficient ofa 
beach in nature, that is, it is easy to find the value of p. It is certainly more 
difficult to determine the damping of a wave. Damping is particularly impor- 
tant in laboratory scale tests in which the motion in the boundary layer is 
likely to be laminar. In nature, as a first approximation, this effect could 
probably be neglected, that is r = 1, except when a wave trap or the presence 
of lateral beaches could produce important frictional effects. 


CONCLU SION 


The motion of a wave on a discontinuity has been analyzed previously.4 
In this paper, considering the entrance and the back of a harbor as disconti- 
nuities, a theory has been formulated that permits the computation of the 
value of the agitation in a basin of relatively simple shape. Such an idealized 
basin is rarely met in practice, but a number of cases may be considered as 
very similar. This theory is valid to explain not only some phenomena of 
agitation in harbors but seiche motion and even tidal motion in bays and on 
the continental shelf. 

The main interest of the theory is to show quantitatively in an idealized 
case and qualitatively in some more complex cases the importance of a num- 
ber of parameters such as the shape ofthe entrance, the slope of the terminal 
beach, and also the importance of multiple reflections and the number of wave 
lengths contained in a basin. It has been shown that for the case of resonance, 
the agitation in a basin could become greater inside the harbor than in the 
open sea, 

Only the main topics have been emphasized inthis paper, but it is possible 
by the set of formulas presented herein to study the previously unknown in- 
fluence of a very great number of other parameters. 

The power and relative simplicity of the method of computation by come 
plex numbers in the analysis of a case that, even though idealized, is fairl 
complicated, should be emphasized. 

This theory is limited to the computation of the agitation in a harbor sub- 
jected to a wave of constant amplitude and constant period. In the field, waves 
are irregular. Hence, to be completed this study must be followed by anothe 
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in which the probability of obtaining a degree of agitation is to be computed 
as a function of the incident wave spectrum. 
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APPENDIX.—NOTATION 


The following symbols, adopted for use in this paper, are listed here for 
easy reference and for use by discussers: 


Ay = complex number of modulus, 

Ay = (semi-amplitude of the incident wave) and argument, 

a = (phase of the incident wave at the entrance to the harbor). 
A = shoaling coefficient (see Eq. 28). 

Ag = agitation at the back of the basin, 


Ag = agitation with resonance. 
max 


Ag nin = minimum agitation, opposite of resonance. 


B, = complex number of modulus. 

By = (semi-amplitude of the reflected wave seawards) and argument, 

By = (phase of the reflected wave at the entrance to the harbor), 

b = length of the entrance. 

C, = complex number of modulus. 

C. = (semi-amplitude of the wave travelling towards the basin at the en- 
trance) and argument-phase. 

/\ 

Cy = at the entrance. 

Co = complex number of modulus: Cg and argument: Co : half the wave 

height and phase of the incident wave at the back of the basin, re- 

spectively. 

Dy = complex number of modulus. 

Dy = (semi-amplitude of the wave arriving from the basin to the entrance) 

: and argument-phase. 

= at the entrance. 


NN 
= complex number of modulus: Dg and argument Do: half the wave 
height and phase of the reflected wave at the back of the basin, re- 


spectively. 


> Ppl a> 


> AP oP! 
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= depth of the open sea. 
depth of the basin. 


A 
= complex number of modulus: i=1, and argument: i= 0. 


coefficient of friction. 


= wave length at the open sea and: my = a ; 

= wave length in the basin and: m2 *= ae 

= width of the open sea. 

= width of the basin. 

= length of the basin. 

= a positive integer: n = 0, ; ee fe 

= coefficient of reflection of the back of the basin- _Kl 
= coefficient of decay of the wave-height in the basin r = i : 


complex number of modulus: r and argument: S 


difference of phase between the waves at the two extremities of the 
basin. 

= change of domain coefficient (see Eq. 26). 

= complex number of modulus: aj and argument: 4. 


= coefficient of transmission of the wave Ay: ratio of the transmittec 
wave height to the incident wave height: Ay. 


= change of phase of the transmitted wave at the entrance towards the 
basin. 


= complex number of modulus: @» and argument: @o. 


= coefficient of transmission of the wave Dy: ratio of the transmitte: 
wave height to the incident wave-height Dj. 


= change of phase of the transmitted wave at the entrance seaward 
= complex number of modulus fj and argument By. 


= coefficient of reflection of the wave: Ay: ratio of the reflected wav 
height seawards to the incident wave height: Ay. 


= change of phase of the reflected wave height at the entrance goa 
wards. 


= complex number of modulus: fs» and argument: Bo. 


= coefficient of reflection of the wave: Dy: ratio of the reflected wav 
height towards the basin to the incident wave height: D,. * 


= change of phase of the reflected wave towards the basin at the er 
trance. 


4 
= length of the beach at the back of the basin. © 
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ESTUARIAL SEDIMENT TRANSPORT PATTERNS 


By H. A. Einstein,! F. ASCE, and R. B. Krone2 


SYNOPSIS 


The sediment deposits most characteristic of estuaries are muds which con- 
tain a large percentage of clays. These clays givethe mud its particular qual- 
ities. The transport of these muds is largely determined by flocculation char- 
acteristic of suspended clays. The flocculation process is described and an- 
alyzed. The process of deposition and scour of the bed of these sediments is 
also described. All of these processes are finally combined in a description 
of the process of shoaling in Mare Island Strait, an artifically deepened side 
arm of San Francisco Bay. 


INTRODUCTION 


Sediment movements in a large estuary or bay system appear to an inter- 
ested observer to have several curious aspects. Sediment materials, trans- 
ported by land drainage, enter an estuary in varying amounts with the inflowing 
waters. Except during exceptionally high river inflows, little sediment may 
leave the estuary tothe ocean. It is apparent that river-borne sediment is re- 
tained by the estuary. During periods of low sediment inflow to an estuary, 
however, shoaling inchannels and harbor areas may occur at much larger rates 


a ee 
 Note.—Discussion open until August 1, 1961. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE, This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961. 
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than are observed at periods of high inflow rate. The amount of sediment re- 
moved for maintenance of navigable channels often exceeds the total sediment 
inflow to the estuary. That the estuarial sediment movement consists of more 
than the deposition of the currently river-borne sediments is shown also by 
varying degrees and locations of turbidity, often visible at the water surface, 
that occur at times unrelated to the varying river inflow. 

The estuarial environment obscures even the simplest information on sedi- 
ment movement. Tidal fluctuations at the mouth of an estuary or bay are pro- 
pagated throughout the estuary with velocities and amplitudes that depend on 
local bottom and shoreline configurations. These tidal changes, often modified 
by resonance together withriver inflow, result in roughly periodic, continuously 
changing circulation patterns withinthe estuary. Water current velocities range 
from zero at slack water to high velocities, particularly at narrows, during the 
intermediate periods. 

Salinity changes are also characteristic of an estuarialenvironment. During 
high runnoff periods, fresh river waters may extend to the ocean. During nor- 
mal or low flows, however, estuarial waters are commonly saline, with salt 
contents near that of sea water in the estuary and for some distance up inflow- 
ing rivers. The salinity often varies with tidal currents, particularly near the 
river mouth. Regular wind patterns are often characteristic of the coastal a- 
reas of whichestuaries are a part. In additionto modifications of water depths 
and currents by wind, even moderate winds generate wave disturbances that 
seriously alter the character of the sediment transporting medium in shallow 
water. Extensive shallows formed by accumulated sediment and cut by sparse 
dendritic drainage channels are characteristic of bays that might by their na- 
ture be expected to have an important role in the estuarial sediment regime. 

The San Francisco Bay system in California, for example, has all of the 
characteristics mentioned. Drainage from the entire Central Valley area is 
discharged at the northeastern end of the bay, and numerous small streams 
draining the neighboring area enter the bay at widely distributed points. The 
bay system covers an area of 485 sq miles, and 70% of this area is less than 
15 ft deep at mean lower low water. Large areas less than 6fft deep occur. 
The bay is about 55 miles long, extending both north and south of a narrow 
mouth. Resonant tidal variations, excited predominantly by diurnal tide, con- 
tribute to water interchange between the north and south bays, and river dis- 
charges modify the northern current patterns. 

Many similar bays exist, and qualitative aspects of the sediment regime in 
this bay should also be applicable to them. Recent studies of the sediment be- 
havior in the portions of the bay near the entrance of the predominant land 
drainage have revealed the sediment transport processes that comprise the 
sediment regime. This paper presents a brief description of the factors that 
largely determine this regime and the participation of these factors in over-all 
sediment movement patterns. 


THE SEDIMENTS 


Samples of sediment from the bottom of points distributed around the bay 
system indicate that, except at the mouth of the predominant inflow channels, 
the sediments are composed 50% or more of clay minerals with the largest 
particle size ranging upto about 100 mv. The mineral composition of the fines 
supplied by the watershed includes many clay minerals of which mostare clas- 
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Sified as illite and mortmorillonite, with smaller amounts of the kaolinite group 
added. The bay sediment shows the same composition. 

These minerals respond rapidly to changes in their hydraulic and chemical 
environment because of their small size, their large specific surface areas, 
and their crystal structures. The minerals are commonly thin plates, flexible 
sheets, or tubes, andrange in effective diameters from a few micronsto small 
fractions of a micron. The arrangements of atoms in these minerals are such 
-hat the large surfaces of the minerals have a negative charge, and when sur- 
rounded by water containing salts, cations accumulate near the surfaces ina 
cloud. In waters containing small amounts of salt, these clouds are large and 
interact electrically, preventing the particles from coming close together, Such 
suspensions are called “dispersed” and are typical of suspensions in many 
streams. Since small particles settle only slowly, dispersed particles tend to 
remain permanently in suspension in flowing water, and represent there part 
of the “wash load.” 

When these particles are suspended in water containing a gram of sea salt 
ver liter, or more (sea water contains approximately 34 g per 1), the effect of 
he cation cloud is repressed by the abundance of anions in solution, and par- 
icles do not repel one another at the larger distances. They can approach one 
another freely, and an attractive force, that operates only at short range, will 
oredominate and mutually attach the particles. They then stick together, and 
f successive particle collisions occur, a multi-particle aggregation, or floc, 
vill form. The shapes of clay particles cause the structure of a floc tobe open 
ind largely to contain water as does a sponge. The reduced density and the in- 
sreased size change the hydraulic character of suspended sediment. 


FLOCCULATION 


The two requirements for the formation of flocs from dispersed suspensions 
re the reduction of the repulsive forces and the movement of one particle to 
nother. The first requirement, the salinity of at least 1 g per 1, occurs in 
nost of an estuarial system to the extent that nearly every collision results in 
dherence. Collision of particles can occur in several ways: (1) by Brownian 
notion, or random thermal movements of small suspended particles; (2) by 
nternal shear motion of the water; and (3) by differential particle settling ve- 
acities. Since flocculation determines the hydraulic transport character of 
uspended clayey sediment in an estuarial environment, it is useful to consider 
he relative importance of these different ways of bringing particles together. 

A theory of rapid flocculation was first given by Smoluckowski in 1917-18. 
‘his theory, and extensions by later workers, is given by J. Th. G. Overbeek. 3 
mly simplified formulas are used here to illustrate the relative importance of 
he flocculation mechanisms under conditions which cause every interparticle 
ollision to result in a stable multiplet. 

_ The probability of collision from Brownian motion, I, can be approximated by 


poe bal te Se re Re ee 
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3 «Kinetics of Flocculation” by J. Th. G. Overbeek, Colloid Science 1, Elsevier Pub- 
Shing Co,, Amsterdam, 1952, pp. 278-3-1. “ 
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in which k represents Biltzman’s constant, 1.38 x 10716 ergs; T denotes ab- 
solute temperature, in K; vo is the number of primary particles per cu cm; 
and 7 represents the viscosity, in poises. Suspended sediment concentrations 
in channels of San Francisco Bay commonly range between a few and 100 parts 
per million. Taking 0.3 mv as a representative spherical diameter and 15 C 
as a typical temperature, I = 0.038 per sec. As primary particles join other 
primary particles and flocs, the total number of multiplets is reduced, increas- 
ing the distance between them and reducing the frequency of collisions. 

Each time a collision occurs to all particles, the number of particles is re- 
duced to one half of the previous number, but the particle weight is doubled. 
The rate per unit time at which the particles grow in mass during any part of 
the growth period is proportional to their size, but inversely proportional to 
the time between collisions, whichis itself inversely proportional tothe number 
of particles or proportional tothe weight of the particles. This means that the 
rate of particle growth is constant, or for the prior given 0.3-mv, original size 
particles, the growth for any size floc in 26 sec is equal tothe mass of an orig- 
inal floc. 

The time to build particles of 100-mv average size andof a specific gravity 
of 1.10 from particles of 0.3-mv diameter and 2.65 specific weight is 


100\3 (1.10-1) _ goo pk 
26 (5) 26521) ~ 58 x 10° sec = 670 days 


while the growth from an average size of 0.3-mv to a floc size of 10 mv takes 
only 16 hr. This process, which does not require any outside driving force, 
must be expected to create the smaller floc in a short period, wherein the large 
flocs will probably appear only if one of the other effects enters the picture. 
Only a much faster build-up of flocs could actually maintain an average size 
much larger than 10 mv if the natural breakage of these rather fragile struc- 
tures is taken into consideration. 

Collisions of suspended particles or flocs resulting from internal shear of 
flowing water can be calculated by finding the probability of one particle pass- 
ing through a sphere of radius R within which a lasting bond between particles 
is formed. The simplest spherical radius would be the sum of the radii of the 
colliding particles. For a water shearing rate, du/dz, and a concentration of 
particles V, the probability of collision J is 


A comparison with the probability of collision from Brownian motion yields 
the ratio 


PACS Ho iter ce ere see ecco ees (3) 


The relative importance of internal shear strongly grows with the floc or par- 
ticle size. For a shearing rate of 1 per sec and a particle or floc diameter of 
2 mv, the ratio would equal one. A 10-mv floc would already give J/I = 120. 
It can be concluded that after flocculation by Brownian motion has begun, in- 
ternal shear rapidly accelerates flocculation. 
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An internal shearing rate of 1 per sec probably does not occur in open flow 
except near the boundaries and in strong turbulence. The bulk of sediment sus- 
pended in the flow must form flocs of moderate size, therefore, without sig- 
nificant increase in flocculation rate by shear. Under wave conditions, unstable 
flow, flow around pilings or other objects, the flocculation rate, enhanced by 
local turbulence, can “snowball” and produce flocs of sizes that rapidly settle. 

When particles having a wide range of sizes are suspended together, larger 
ones settle faster than smaller ones, and increased frequency of collisions oc- 
cur. H. Miller4 studied this mechanism and showed that a significant effect 
is found when the radius, r, of the settling particle is 


r= 


E ery 
ngp 


in which g is acceleration of gravity and p is the density of the particles. The 
term in brackets corresponds to a 5.0-mv diameter solid particle or a 6.3-mv 
floc. The particles caught must be larger than 


which corresponds to a 2.1-mv solid particle, or toa 2.6-mv floc. It can be 
seen that this mechanism is most important for collection of small flocs by 
larger ones. This mechanism intensifies the formation of large flocs after ac- 
celerated flocculation by shear, and further enhances deposition near pilings 
or flow-disturbing works. 


DEPOSITION 


While the individual clay particles are so small that Brownian motion alone 
is almost able tokeep them in permanent suspension, the flocs of such particles 
are sufficiently great and heavy to make them settle out in still water. Large 
flocs have a settling velocity measured in feet per hour. Even the slack water 
periods in tidal areas are sufficiently long to permit considerable amounts of 
flocculated clay to accumulate on the bed. 

Laboratory experiments show that at a density of about 10 g per 1, the clay 
flocs form a continuous structure. This density is obtained by settling of in- 
dividual flocs. Settling near this density occurs at a much reduced velocity 
compared to that at low concentrations, because of the increased resistance of 
the waters escaping from betweenthe many flocs. One speaks here about * hind- 
ered settling.” At a density of 10 g per 1, the flocs make bodily contact and 
form a continuous cohesive mass of the same density and strength as the in- 
dividual large flocs. This mass can support a measurable shear stress with- 
out permanent deformation, at least for restricted duration. Upon prolonged 
stress, it will very slowly creep. Such creep may be either a shear motion or 
a compaction, or both, depending on the type of stress applied. Compaction is 
particularly important because it causes the strength of the deposit to increase 
aS more and more bonds are made. 


_ 44, Muller, Kolloidchem, Beihefte, Vol. 27, No. 223, 1928. 
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FIG, 1.—RADIOACTIVITY DISTRIBUTION 


Deposition occurs not only in still water, but also in flowing water. Up toa 
flow velocity of 0.5 fps, deposition is the same as in still water, that is, the 
entire flocs are deposited at and bonded to the bed. As the velocity increases, 
some of the loosely attached particles are occasionally torn off the bottom while ~ 
smaller flocs or parts of larger flocs can still deposit on the bed. A denser 
and tougher layer of clay develops at the bed surface, but at a much slower 
rate. These stronger bonds between clay minerals canbe visualized to develop 
where they attach their flat sides to one another, in contrast to edgewide con-— 
tact in the loose flocs. The more flat contacts are made, the stronger becomes 
the structure. A wide range of such structures occurs in nature, each being 
characteristic for a particular velocity of flow during deposition. This picture 
is further complicated by secondary compaction of the deposits due to gravity. 


SCOUR 
The scour of these clay deposits has not been studied with as much detail 
as has their deposition. It is unknown whether scour is caused predominantly 
by the shear stress at the bed or by pressure fluctuations. It has been ob=_ 
served, however, that scour does not occur by chipping large complexes from 
the bed, but rather by the continuous wearing off of small particles. The speed 
of scour is reduced with increasing concentration of the flowing water, indi 
cating atype of equilibrium condition to be approached between scour and ee 
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CONTOURS IN HUNDREDS OF 
COUNTS PER MINUTE 


TOM OF MARE ISLAND STRAIT 


AT BOTTOM OF MARE ISLAND STRAIT 


sition. No attempt has been made, so far, to define the conditions of such an 
quilibrium. 


FIELD OBSERVATIONS 


In returning to the problems of the previosuly quoted San Francisco Bay, 

ome particular observation may be described that resulted from a prototype 
tudy which had the purpose of providing a basis for the verification of the 
Yorps of Engineers’ hydraulic model of the Bay at Sausalito, Calif. One of the 
nost important problems is the very high rate of mud deposit in Mare Island 
trait, the lower endof the Napa River channel, just above its junction with the 
nain Bay system. This junction is at Carquinez Straits, a narrows of the Bay 
n which the fast tidal currents prevent deposition of mud. The deposited sed- 
ment was identified as identical with that of the bay proper, not with that of 
he Napa River. 
: The most perplexing observation was made in connection with the time at 
yhich the sediment was deposited. It was found that deposition occurred pre- 
iominantly during the summer months when the bay water contains the lowest 
ediment concentration and when the rivers carry practically clear water. It 
as, thus, necessary to find the immediate source of the material somewhere 
n the bay, even if the material originated from the major river system. 


i 
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In order to establish the exact process by which the sediment enters Mare 
Island Strait, a number of specific observations were made. First,a suspended 
charge of bay sediment was released at the beginning of flood tide at the mouth 
of Mare Island Strait, after it was identified by radioactive gold label. This 
trace sediment was found to have been almost fully deposited inthe straits dur- 
ing the first tidal cycle. It was found, furthermore, not to have moved signifi- 
cantly during the next few days during which the activity was still sufficient to 
be identified. Within the Strait, distribution of the radioactivity was the same 
as that of the normally deposited sediment (Fig. 1). Characteristically, the 
highest rates of deposition occur on the two sides behind piles andother struc- 
tures which provide local turbulence that causes fast growth of large flocs. 
This is shown in Fig. 1 by the 500-count-per-min contour along the west shore 
above the finger piers. The experiment also showed that the majority of the 
shoaling material enters in suspension through the mouthof Mare Island Strait, 
which enters the area from Carquinez Strait. 

Next, it must be asked, from where: may suspended sediment originate dur- 
ing the summer months. Immediately west (this is, downstream) from the 
junction of Mare Island Strait, the main channelof Carquinez Strait widens into 
San Pablo Bay. This large body of water has wide, shallow, tidal flats in which 
the regular summer winds generate strong waves. These break near the en- 
trance to Carquinez Strait and suspend large amounts of sediment. The tur- 
bidity at the surface, resulting from this action, can regularly be observed 
from the air. 

Finally, another test with radioactivetracer was madeto study the removal 
of sediment from the tidal flats during periods of wave action. Bay sediment, 
labeled with radioactive gold, was placed in suspension above the bottom during 
a quiet period. The tracer was placed over two shallow areas of San Pablo 
Bay, and part of the tracer was deposited with the natural sediment. These 
deposits were established at different depths. Later, a moderate wave con- 
dition occurred and it was found to have disturbed the shallow deposit some- 
what, but to have left the deeper deposit undisturbed. The disturbance appeared 
in the form of elongations of the radioactivity pattern not in the direction of 
the wave action, but in the direction of slow local currents. The same effect 
had been observed previously in the laboratory. Wave action will suspend the 
sediment in the water, but does not move it a significant distance. Transport 
of sediment from the mud flats is predominantly by suspension in moving water. 


CONCLUSIONS : 


During studies of the motion of flocculated sediment in estuaries, both in 
the laboratory and in the San Francisco Bay system, it was found that the clay 
sediments flocculate slowly as they advance from the rivers into water with 
salinity of 1 g per 1 and more. A faster development of the largest flocs which 
settle rapidly is caused by a general turbulence which exists in the extended 
tidal flats. This sediment supply occurs during winter floods. “ 

During almost daily summer winds, some of these winter deposits are stirred 
up by wave action, and high concentrations of flocculated particles occur. Thes 
are moved by local currents into the quiet areas of the artifically deepene 
channels, particularly where pilings and other harbor structures cause a 
formation of settleable flocs. Here the sediment deposits and necessitates “7 
} 
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1sive maintenance dredging programs for the preservation of navigable ship- 
lanes, 
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VIBRATION PROBLEMS IN HYDRAULIC STRUCTURES 


By Frank B. Campbell,! F. ASCE 


SYNOPSIS 


The paper is principally concerned with the vibration of control gates and 
regulating valves commonly used in hydraulic structures, such as dams. 
The principal findings of 12 yr of vibration investigations by the Army Engi- 
neers civil works activities are reported. 

The classical theory of vibration is reviewed to select those portions ap- 
plicable to the subject. Evidence of Coulomb damping is demonstrated and 
the results of field tests are analyzed todetermine the constant friction force. 
Both laboratory and field tests are reported. 

Emphasis is placed on exicting forces expected to be found in hydraulic 
structures. The Von Karman vortex trail, self-excitation involving reflecting 
pressure waves and other hydraulic pulsating phenomena are examined. Cer- 
tain practical conclusions are drawn. 


INTRODUCTION 


The study of vibration theory has had an interesting history. The physical 
phenomena involved with musical instruments was of interest to the scientist 
of the eighteenth century. J. L. Lagrange solved the physical problem. 
- Note.—Discussion open until August l, . To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol, 87, No. HY 2, March, 1961. 

I Chf., Hydr. Analysis Branch, Hydraulics Div., U.S. Arm, 


ment Sta., Vicksburg, Miss. 
_ 2 Recherches sur la Nature et la Propagation du Son,” by J. L. oaeetyaa Couvres 
759). ‘ 


de Lagrange, Vol. 1, Paris, 1867, p. 39 (Miscellanea Taurinesia, t. 1, 1 
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of the vibrating string in 1759. A great advance tothe physics of vibration was 
contributed in the nineteenth century by Lord Rayleigh.3 Many refinements to 
Rayleigh’s work have since been accomplished. 

It is the purpose of this paper to emphasize to civil engineers, the impor- 
tance of certain types of mechanical vibrations inhydraulic structures, These 
particular vibration problems are not adequately covered, evenin mechanical 
engineering literature. 

Scope.—Special emphasis is placed on the problems of exciting forces 
which arise from certain hydraulic phenomena. 

The vibration problems arising in connection with the operation of gates 
and valves commonly employed to control and regulate the flow of water 
through or over dams will be stressed. The word operation is used because 
much research needs to be done before a hydraulic structure can actually be 
designed with confidence in the knowledge that the mechanical elements in- 
volved will not vibrate under certain conditions of operation. 

The writer is associated with prototype investigations to determine causes 
of vibrating gates or gate elements. Therefore, this paper will constitute, in 
a sense, a Summary report of the U. S. Army Engineers experiences of un- 
desirable vibration, which have been brought to the attention of the Waterways 
Experiment Station. The problems of vibration in hydraulic machinery will 
not be treated, nor will the problems of water hammer and surge tank design. 


APPLICABILITY OF CLASSICAL THEORY 


It is important to examine the classical theory of vibration to determine 
which portions of the science are useful to the solution of the specific prob- 
lems which might be encountered in hydraulic structures. The theory may 
be found in various publications such as those of C. R. Freberg, E. N. Kemler, 4 
J. P. Den Hartog,° S. Timoshenko, D. H. Young,6 M. ASCE, L. S. Jacobsen, 
and R. S. Ayre,? M. ASCE. 

A definition of symbols is shown in Fig. 1. The case of undamped free 
vibrations is used for the purpose of defining basic physical quantities. The 
device shown is a mass, defined as W/g, with an elastic suspension having 
a spring constant K, expressed in 16 per in. In this case, assume a cube of 
steel supported by a tension spring. There will be an initial static elongation 
5s- ) 

If the cube is forceably pulled downward and released, a vertical vibration 
will ensue. The basic differential equation is 7 

WwW . 


ek PER = Ora. wk pee . eee 


The symbol x signifies the second derivative of displacement with respect to | 


time. This derivative is recognized as acceleration. The first term is, 
AE a ae eel Del Scat hata ot neers ata Oe el beta er eee AS 

3 «The Theory of Sound,” by J. W. S. Rayleigh, First Amer. Edition, New York, Do- 
ver, 1945. 


4 “Elements of Mechanical Vibrations,” by C. R. Freberg and E, N, Kemler, Second 
Edition, John Wiley and Sons, Inc., New York, 1949, 

5 “Mechanical Vibrations, "» by J. P. Den Hartog, Fourth Edition, McGraw-Hill Book 
Co., Inc., New York, 1956, 


6 “Vibrations Problem in Engineering,” by S. Timoshenko and D, H. Young, Third 
Edition, Van Nostrand, New York, 1955, 


7 “Engineering Vibrations; 2 by. L, S, Jacobsen and R.S, Ayre, First Edition, McGraw- 
Hill Book Co., Inc., New York, 1958, 
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therefore, the product of mass and acceleration, which represents the inertial 

force of the vibrating body. The second term has the product of the spring 

constant and the displacement, which is referred to as the restoring force. 
The useful solution to the basic equation is 


One period of the cosine curve is 2 7 or wp t. The subscript n is used to 
denote the natural frequency of the elastic system. The symbol w is commonly 
ealled the angular frequency although various writers use different symbols. 


(a) Coulomb Damping 
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It may be readily seen that the natural frequency is inversely proportional 
to the square root of the initial static elongation, 5g. For the elastic system 


described, this value is expressed as 
W (6) 
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The phase angle circle is included to demonstrate the initial displacement 
position. Certain writers, notably Jacobsen and Ayre,? have developed 
graphical phase plane solutions which appear to be useful in certain types of 
damped vibrations. 

Damped Free Vibrations.—When a restraining force is imposed on a vi- 
brating system, the type of damping must be known before analysis. The basic 
differential equations for the three common types of damping are: 


Constant friction (no single complete solution) 


© X+KX + Ce (sink) = 0% 30a Gentteatas (8) 
Viscous damping 
a 2 Ck Pee oe ered eee (9) 
with the solution: 
i olution c.g 
Kimcerd se, (A cos @n’thvars soa etees (10) 


Turbulent damping (no direct simple solution) 
= 0) 
7x +Ga(X) KKM = 0 A OA Oe (11) 


Their applicability to the problem of vibration in hydraulic structures is 
important to this paper. It should be noted that each of the differential equa- 
tions contains inertial force involving mass and acceleration plus restoring 
force involving the spring constant and displacement. To each of the equations 
a friction or damping force is added. 

The constant friction force has the sign of the velocity. The force always 
opposes the motion and, thus, changes sign at each half cycle. Therefore, no 
single complete solution can be written (Eq. 8). Each part of the solution must 
be considered between the largest positive and negative displacements of 
each half cycle. 

The second term of the viscous damping equation (Eq. 9) contains the 
damping coefficient which is multiplied by the first power of the velocity. 
The general equation has a relatively simple and complete solution. The first 
factor on the right contains an exponential function which involves the damp- 
ing constant C, and mass, The second factor contains a cosine function of 
the natural frequency. A large part of the classical theory of vibration is based 
on viscous damping. Unfortunately, the viscous damping solution does not 
seem to be generally applicable to problems encountered in high head hydrau- 
lic structures. ( 

Turbulent damping (Eq. 11) isthe third case to be considered. The turbulent 
damping force is a function of the square of the velocity. This is the common 
type of force with which hydraulic engineers ordinarily deal in their problems 
of pipe friction loss. Unfortunately, the equation is a secondorder, second de- 
gree differential equation andhas nodirect simple solution. W. E. Milne®,9 has 
obtained an elaborate exact solution. He has also prepared tables for the 
evaluation of the solution which have been found useful in surge tank problems 


Soe eed Vibrations,” by W. E, Milne, Univ, of Oregon Publication, Vol, 2, No, 2, 


2) 4 
9 «Table of Damped Vibrations,” by W. E. Milne, Mathematical Series, Vol. 1, No. 1, 


Univ. of Oregon, March, 1929, 
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for hydro-electric power plant design. The technique involved is beyond the 
scope of this paper. 

The problem of attenuation of displacement in damped free vibrations is con- 
‘Sidered next. Constant friction is characterized by an arithmetic decrement 
as indicated in Fig. 2(a). It is important to note that only the constant friction 
force and the spring constant are involvedin the determination of the arithme- 
tic decrement. This is an important relationship as will be shown later. 

In the case of viscous damping, the attenuation for each cycle is called the 
logarithmic decrement and can be evaluated when the viscous damping co- 
efficient, the mass, and the frequency are known. The attenuation is affected 
by an exponential law so that the displacement is damped rather rapidly in the 
‘beginning, as demonstrated on Fig. 2(b). After the damping factor becomes 
negligible, consider that a steady state exists. Reference willbe made to steady 
‘state vibrations, subsequently. 

Forced Vibration.—The aspect of a damping force which decreases the 
energy of the vibrations has been considered. Itis now appropriate to consider 
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FIG. 3.—FORCED VIBRATIONS- FIG, 4.—FORCED VIBRATIONS- 
VISCOUS DAMPING CONSTANT FRICTION 


an exciting force which can increase the energy and hence, the amplitude of 
the vibrations, The basic equation for forced vibrations with the classical 


viscous damping is 


WX +Cy+K + KX = F cos (uy t) eet doko Ss eae eee 


The solution shown in Fig. 3 is forthe case of steady state vibration. In other 
words, the damping effect of the exponential term has become negligible. The 
exciting force F is considered to be applied with a cosine variation. The 


solution for the steady state vibration is 


(13) 
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The first factor of the solution is the ratio F/K, which represents the static 
elongation that would be produced by the maximum exciting force with the 
existing spring constant. If the cosine function in the numerator is placed 
equal to 1, the factor one would be divided by the radical, which is ordinarily 
known as the Magnification Factor. Note that the Magnification Factor involves 
ratios of forcing frequency to natural frequency for the various damping 
ratios Cy/Cy¢. The denominator Cy, of this ratio is the critical damping coef- 
ficient which can be evaluated according to 


Cyc = 2 Wy en eae (14) 


The critical damping coefficient involves only the natural frequency of the 
vibrating system and its mass. 

In considering the subject of damping friction, it would be appropriate 
to quote from Jacobsen and Ayre! as follows: “Friction is therefore, a 
priori, a more complicated property to deal with than either inertia or 
restoration.” 

A graph can be constructed of the Magnification Factor as a function of the 
frequency ratios. The various damping ratios are represented by a family 
of curves. For a damping ratio of zero, the curve represents an undamped 
system, Near the condition of resonance where the ratio of frequencies is 
equal to one, the Magnification Factor and the amplitude approaches infinity. 
The system is safe against wildly fluctuating vibrations in the resonant range 
if the damping ratios are substantially greater than 0.3. 

Constant Friction.—As previously mentioned, the differential equation for 
the case of constant friction has no simple continuous solution. Nevertheless, 
Den Hartogl1 developed a remarkable solution for maximum displacement in 
the case of steady state vibration with constant friction, His equation for maxi- 
mum displacement is 


F 
D=& C0" ents) fame ces (Se 
in which 
“ 1 
A= a ocpiisivnn ip (e! « Ue Dhetiareynet aoe 
f 
eee 
and . 
Wn T Wy 
B= ta 6 goa et pM chs 46 se ee eee wees ) 
Wf +5 2 we (i? 


The solution again contains a ratio of exciting force to spring constant. The 
variables A and B under the radical are two different functions of the ratio of 
forcing frequency to natural frequency of the system. The second term under 
the radical contains a square of the ratio of the constant friction force to the 


OI, 
10 “Engineering Vibrations,” by L.S. Jacobsen and R.S, Ayre, First Edition, McGraw- 
Hill Book Co., Inc., New York, 1958, p. 196, | 
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exciting force. The magnification can be expressed as the ratio of the maxi- 
‘mum displacement to the static elongation attributable to the exciting force. 
Magnification = _ 2 i eT ay Re Ri (18) 


K 

When the second term under the radical (Eq. 15) becomes larger than the 
first, the evaluation of maximum displacement involves the Square root of a 
negative number. Eq. 15 is only of value in the determination of magnification 
above the dashed line shown in Fig. 4. The family of curves represent various 
‘ratios of the constant friction force to the maximum exciting force. For a 
ratio of forces greater than 7/4, the magnification approaches infinity. 

By a complex procedure, Den Hartog11 was also successful in evaluating 
‘Single points underneath the dashed line andhe constructed curves represent- 
‘ing high force ratios. The friction is sufficient in the region of high force 
Tatios to cause a stoppage in the vibratory motion. It is interesting to note 
that the two lines representing force ratios of 0.85 and 0.95 from the Den 
Hartog investigation have maxima at frequency ratios substantially smaller 
than the resonant frequency. It might be expected that, in the case of constant 
friction, if the force ratios are considerably greater than 0.8, there would be 
little likelihood of developing large vibratory amplitudes, Inthe case of viscous 
damping, we surmised that the situation would be fairly safe if the damping 
ratios exceeded 0.3. 

The engineer is then faced with the problem of determining the constant 
friction force and the exciting force. Some test data are available on friction 
factors, which should be applicable to hydraulic gates. On the other hand, 
experimental data on the magnitude of exciting forcesin the various hydraulic 
phenomena are practically nonexistant. The problem of excitation in hydraulic 
structures will be presented later. 

Hydraulic Gate Friction.—A substantial amount of work hasbeen conducted 
on sliding friction and rolling friction research. However, much of the recent 
experimentation has been concerned with the effect of various types of lubrica- 
tion on journal friction. It is believed that ordinary lubricants are not effective 
over a long period of time whenusedon submerged high head hydraulic gates. 

The U. S. Bureau of Reclamation!2,13 has conducted tests on both sliding 
and rolling friction, in connection with their studies of gates for hydraulic 
structures, These laboratory tests were conducted in thedry rather than with 
water surrounding the element as would be expected in an outlet works, for 
example. It appears that the coefficient for sliding friction is in the neighbor- 
hood of 0.1, whereas the coefficient of rolling friction may be approximately 
0.001 depending on the type of metals used. In the field of hoist design for 
high head gates, provision must be made for adequate capacity to break the 
seizure of metals after the gate has been closed under a high head for a 
substantial period of time. Therefore, with application ofa judiciously chosen 
factor of safety, the coefficients mentioned would probably be entirely ade- 


: : 5 ; 5 man Paden 
11 “Forced Vibrations with Combined Coulomb and Viscous Friction,” by J . P, Det 

Hartog, Transactions, ASME, Paper APM 53-9, presented at an June 1931 National Con- 
vention of Applied Mechanics at Purdue Univ, in Lafayette, ‘ 
_ 12 «Tests on Rollers,” by N. G. Noonan and W. He Eaness atl S. Bur. of Reclam., 
Tech, Memorandum No. 399, Separate No. 26, September 26, 5 
_ 13 «Report of Tests on Coefficients of Friction” by N. G. Noonan and W. H, Strange, 
U. S. Bur, of Reclam., Tech. Memorandum No, 432, January 30, 1935. 
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quate for purposes of hoist capacity design. On the other hand, it is believed 
that further investigation should be made of frictionforces in the prototype to 
determine a more accurate value for use in the problems of gate vibration. 


GATES ON ELASTIC SUSPENSIONS 


Model Studies.—Hydraulic model studies of amassonan elastic suspension 
have been found to be practical and useful in a study of exciting forces. The 
Waterways Experiment Stationl4 made studies for the Omaha District beginning 
in 1947. These studies simulated the early Garrison control gate design as 
well as the Fort Randall control gates. 

The mass of the gate in the prototype was simulated in the model, An ex- 
tension spring was selected which produced a static elongation similar to that 
expected in the prototype. Considerable care was taken to reduce friction in 
the model by the use of minature roller bearings. Both the vibratory motion 
and the forces on the gates were measured. 


(2) FLAT BOTTOM GATE +x, | O5SEC 
ete Lecesees 
t 
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FIG, 5.—FLAT BOTTOM GATE- FIG, 6.—DEVELOPMENT OF STANDARD 45° LIP 
VORTEX TRAIL RESO- 
NANCE 


One of the objectives of the tests was to determine whether a vortex trail - 
from a partly open flat bottom gate could produce vibrations, The basic theory 
of this investigation is demonstrated in Fig. 5. The initial static elongation of 
the model gate on its elastic suspension 5= 0.0126 ft. Substituting Eq. 5, the 
natural frequency of this system fp = 8 cycles per sec. 

When a flat plate is placed laterally across a stream, the Von Karman 
vortex trail is shed from each edge ofthe plate. The frequency of the shedding 
of these vortices is defined by the Strouhalnumber (St) which involves the freq 
quency, the length of the plate, and the velocity of the fluid. 


feos " 
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14 “Spillway and Outlet Works, Ft. Randall Dam, Missouri River, ” Tech, Report Now 
2-528, U.S, Army Engrs., Vicksburg, Miss., October, 1959, 
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Fig. 5 shows the flat plate to have half of its length protruding into a gate. 
chamber. This represents the partial opening of a control gate. The image 
of the flat plate above the roof of the conduit completes the normal concept of 
the vortextrail. During the early phases of the tests it was notknown whether, 
the vortex trail can be shed from half of a plate as shown in this situation. 
However, it was soon learned that this phenomena can existfor the protrusion 
of a gate leaf into a conduit. 

Various experimenters!9 have determined the Strouhal number for a flat 
plate normal to the flow to be approximately 1/7. If the vortex trail can pro- 
duce an exciting force on the bottom of the gate which wouldcome into reson- 
ance with the naturalfrequency of the elastic system, it would be expected that 
the gate would produce vibrations of large amplitudes. For a 17-ft prototype 
gate opening, L,, = 0.333 ft,and a model velocity V,, = 18.5 fps, the theory in- 
dicated that there would be essential resonance. For a normal plate § = 1/7 
= 0.143 and fg = 7.92,thus the ratio of forcing frequency to natural frequency 
is 

te 9 


Pie “ 
te Sat ian & 0.99 


It was actually determined that whenthe model gate was set to simulate a 17-ft 
prototype gate opening, violent vertical vibrations were produced. Therefore, 
it was obvious that resonance between the frequency of the vortex trails from 
a gate projection acting as an exciting force on an elastic system could actually 
be produced in the laboratory. 

Split Leaf Gates.—In connection with the subject of the vortex trail as an ex- 
citing force, mention should be made of split leaf spillway gates. With this sys- 
tem, one leaf rests on top of the other. In times of flood-flow both leaves are 
removed. 

After the top leaf has been removed and the bottom leaf is in the process of 
raising, water can flow over andunder the partially raised leaf. The full vor- 
tex trail is then effective in producing periodic pulses, alternately, on the top 
and bottom of the gate leaf. 

_ The Waterways Experiment Station studied this problem with a model of the 
Old River Control Structure gates. 

Improvement of Gate Lips.—It was learned from the studies of Ft. Randall 
gate vibration that a flat bottom gate had highdownpull. The tests at the Water- 
ways Experiment Station further showed that a flat bottom gate also was more 
susceptible to vibtations. Various lip extensions from the downstream edge 
were tested. The flat bottom gate and the lip extension schemes are shown in 
Fig. 6(a) and (b). Typical oscillograms translated in terms of vertical dis- 
placement from the accelerometer records are shown in Fig. 6(d) and (e). 
A fairly regular motion of 8 cycles per sec maybe seen in Fig. 6(d), whereas 
no such regularity is shown for the case of a lip extension (Fig. 6(e)). 

- The vortex trail was probably springing essentially from the lip extension 
rather than from the upstream edge of the flat bottom gate. The engineers of 
the Waterways Experiment Station tried the next logical step which was a gate 
with a 45° bottom surface (Fig. 6(c)). It was found that the vortex trail would 


_ 15 «Modern Developments in Fluid Dynamics,” by 8, Goldstein, First Edition, Oxford, 


Marendon Press, Vol. 2, 1938, p. 571, 
, 16 «Downpull Forces on Vertical Lift Gates,” Tech. Report No, 2-477, Report No, 1, 
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then spring from the downstream edge and minimize pulses on the bottom of 
the gate. 

No tendency for vibration was found for the 45° gate lip. Furthermore, the 
downpullor rather the reduction of pressure onthe bottom ofthe gate was much 
less for the sloping gate bottom than for the flat bottom. Therefore, the 45° 
gate lip has become standard on Army Engineer gates because reduced down- 
pull is indicated and there is less tendency towards vibration of the gate. 

Prototype Tests.— Before proceeding with a report of the prototype tests, it 
is appropriate to briefly examine the important problem of the spring constant 
for cable supported control gates. Tests were made by Waterways Experiment 
Station engineers at Knightville Dam in New England in the fall of 1955. In the 
analysis of the data, it was desirable to compute the spring constant ofthe cable 
supports. If the total metallic cross-sectional area ofthe cableandthe modulus 
of elasticity are known, the spring constant can be determined. A fairly accu- 
rate value of metallic cross-sectional area can be furnished by the manu- 
facturers of wire rope. 

The only information on modulus of elasticity of wire rope at that time was 
a value of 12,000,000 psi for unstressed cables furnished many years ago by 
one of the wire rope manufacturers. The use of this value did not produce a 
theoreticalfrequency which agreed with the measured frequency at Knightville 
Dam on tests at three different gate openings ranging from 10 1/2 ft to11 ft. 

The natural frequency measured at Knightville was about 6 1/2 cycles per 
sec. Knowing the other physical data pertaining to the vibrating system, the 
modulus of elasticity was estimated to be approximately 21,000,000 psi. Soon 
thereafter, information became available from tests on prestressed wire rope. 
These results indicated that the modulus of elasticity should be in the neighbor- 
hood of 20,000,000 psi. This value is now considered to be a fairly reliable 
modulus of elasticity for wire rope which has been stressed repeatedly in oper- 
ation. 

It was clearly evident that the mass of the water contained within the gate 
must be added to the mass of the gate metal for a proper vibration analysis. 

Fort Randall Tests.—The Waterways Experiment Station 7 conducted ase- 
ries of three tests, for the Omaha District, on the control gates for Ft. Randall 
Dam in 1954 and 1955, while the reservoir was filling. Fig. 7 pertains to the 
test made hg a head of 110 ft. The spring constant for the cable suspension 
was 1.7 x 105 psi and the gate weighed approximately 47 tons. Gate roller 
tracks and side guides are shown schematically. The side guides restrain any 
large motion of the gate in the upstream direction. Half of the guides were 
spring mounted and the other half were not. 

The oscillograph records show that for large gate openings, there were oc- 
casional vertical vibrations of approximately 4 cycles per sec. These vibra- 
tions continued for several cycles and were then damped out. After an exam- 
ination of the record, it was apparent that the frictional force involved in the 
rollers was rather high, but that occasionally a surge against the downstream 
face of the gate movedthe gate upstream so that momentarily the rollers were 
not in contact with thetrack. During these intervals the gate was permitted to 
vibrate vertically. As the rollers again came in contact, the vibration was 
damped. Further examination of the damping periods of the record indicated 
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Missouri River,” Tech, Report No, 2-435, U.S. Waterways Experiment Sta., Vicksburg, 
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that the damping was of a Coulomb character. One of the damping records is 
shown in Fig. 7(b). It may be noted that straight lines connect the peaks and 
troughs of the curve. 

This was a significant finding, for although experimental evidence was at 
hand to show that the damping was predominantly caused by a constant friction 
force, the theoretical problem was known to be complex. It is pertinent to 
quote again from Jacobsen and Ayrel8 as follows: 


“Tt is therefore clear that the presence of constant friction greatly com- 
plicates the analysis of steady state response to an alternating force, 
especially if the friction is large enough to enforce stops in the motion.” 


The arithmetic decrement was then computed for several damping periods. 
Knowing the spring constant and the arithmetic decrement, the damping force 
(Ce can be computed. 


The section of the recordshown in Fig. 7(b) is for a 20 ft gate opening and the 
damping force was computed to be 2,740 lb. An average of four observations 
(gate openings of 19 ft to 21 ft) indicated a damping force of about 2,800 lb. It 
should be mentioned that in two cases studied for the same gate opening, one 
damping force appeared to be double the damping force for another period. It 
is possible that the lowerforce could have been caused by only one roller track 
coming in contact with the guide. In the case of the higher force, it appears 
that both roller tracks came in contact with the guides simultaneously. 

It was realized that although the theory of Coulomb damping is complex, the 
simple determination of the arithmetical decrement affords a means of meas- 
uring friction. It is believedthat a study of damping records can furnish valu- © 
able information on the total damping force for any particular system under ~ 
observation. There has not been an opportunity to attempt to compare the total © 
damping force measured in the prototype with a computed force based on fric- — 
tion factors as measured in the laboratory. An estimation of the frictional — 
forces of the individual components of a rollertrain may not be a simple mat- 
ter. . 

Sluice Gates.—The hydraulically operated sluice gate is a common device — 
used for controlling the flow of water through concrete dams. Prototype tests — 
were made on the sluice gates of Pine Flat Dam, Calif. for the Sacramento 
District in 1952. Measurement of the vibrations of the gate leaf was a part of © 
these tests. The gate size was 5 ft by9 ft. Fig. 8(a) shows the schematic sec-_ 
tion of the slide gates. | 

The 6 ton sluice gate was supported by a 6 1/2 in. diameter steel stem. A 
typical oscillogram is shown in Fig. 8(b). The maximum amplitude of vibra-_ 
tion was of the order of 6/1,000,000 in. This emphasizes the fact that we are 
dealing with a relatively light weight gate on a stiff suspension. } 

The oscillogram in Fig. 8(c) represents a significant observation on the 
Pine Flat sluice gate. This record represents the vibration existing when the 
gate was being moved from a 3 ft to a 4 ft opening. For such a situation, the 


SS 
18 “Engineering Vibrations,” by L.S, Jacobsen and R.S. Ayre, First Sate Meas 
Hill Book Co., Inc., New York, 1958, p. 226. 
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amplitude exceeded 1/1000 in. or a thousand times that observed for the sta- 
tionary position. It seems reasonable to believe that this larger vibration am- 
plitude was caused by alternate seizure and relief of the metals on the slide 
bearing. This is sometimes called chatter. 

Although the hydraulic head for the test just described was 110 ft, this gate 
was subsequently operated under a headin excess of 300ft. No detectable dam- 
age to the gate was in evidence. A cursory study was made of the fatigue as- 

pectin the gatestem. The stresses in the stem proper were so low as to elim- 
inate any further consideration of fatigue. No study has been made of the pos- 
sible fatigue aspects of the supporting elements at either end of the stem. 


ELASTIC BEAMS 


Flexural Vibration.—The horizontal beams of a tractor gate are normally 
stiff. The Fort Randall gate beam is shown on Fig. 9. 

For a concept of the stiffness, it may be noted that the initial static deflec- 
tion based on the mass of the gate and contained water is only 6x 10-4 in. A 
comparison is shown between the theoretical frequency of the beam with an with- 
out contained water. 

The frequency for the dry condition may be seen to be 205 cycles per sec 
and 140 cycles per sec for the condition of inclosed water. This emphasizes 
the importance of accounting for the mass of the contained water for gates of 
this type. 


PLATES AND SHELLS 


Howell-Bunger Valves.—Some years ago it was learned that weld seams on 
4-vane Howell-Bunger valves at two different projects had failed. Fig. 10 isa 
schematic representation of the vibration characteristic of both a 4-vane and a 
6-vane valve shell. Timoshenko treats the subject of vibration of thin walled 
elastic cylinders. 

Fig. 10(a) shows schematically a 4-vane valve shell vibrating with the fun- 
damental mode and Fig. 10(b) a 6-vane valve shell vibrating with a secondary 
mode. Eq. 21 gives the natural frequency for the fundamental mode and any 
of the higher modes. 


ere gil 
eer a Mj 


It may be seen that the factor involving the radicalis concerned with the mass 
and elasticity of the shell. The factor (Mi) is designated a mode factor. 


ee Nees 


or the fundamental mode (Fig. 10(a)), i= 2, and the resulting mode factor 
9 = 2.68. For the secondary (Fig. 10(b), i= 3, and the value of the mode 
actor M3 = 7.59. Since the naturalfrequency of the shells are proportional to 


atio of the fundamental mode is 2.82. 
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This ignores the existence of vanes within the shell. 


However, it appears 


that by the use of 6 vanes and producing a vibration of the secondary mode type, 


the frequency could be expectedto be much higher than for the 4 vanes and the 
fundamental mode. When the frequency is substantially higher, it would nor- 
mally be expected that the amplitude of vibrations would be lower and that for 
each cycle the total amount of energy is less for the higher modes. 
reason, it is believed that a 6-vane valve is less liable to know failure of the 


(a) FOUR VANE VALVE (6) SIX VANE VALVE 


FIG, 10,—HOWELL-BUNGER VALVE VIBRATION 


—4 Koo “ss SECOND 


HEAD ON GATE 3 FT 


FIG, 11.—SKIN-PLATE VIBRATION— 
SUBMERGIBLE GATE 


seams which connect the vanes to the shell. 


ACCELERATION RECORD 


THEORETICAL 159 CPS 
ACCELEROMETER 150-200 CPS 
SOUND 115-160 CPS 


For this 
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To the writer’s knowledge, the problem of the Howell-Bunger valves with 
the mass and elastic characteristics of both shell and vanes has not been 
treated. It would be interesting indeedto see a thoroughanalysis of the 4-vane 
valve compared with the analysis ofthe 6-vane valve. Furthermore, it is con- 
 ceivable that a 5-vane valve would cause the vibration to go to a still higher 
mode and hence have less liability for failure of weld seams. 

Rayleigh-Ritz Method.—It is appropriate to mention the Rayleigh-Ritz method 
by which problems of this nature can be solved. Lord Rayleigh3 treated only 
the case of the fundamental mode. He devised the method of writing the basic 
equations for both potential energy and kinetic energy. By equating these two 
energies and by proper manipulation of the equation, he determined the fre- 
quencies of the various types of elastic systems. Some years later Walther 
 Ritz19 extended Lord Rayleigh’s method to treat the problems of the vibration 
characteristics with modes higher than the fundamental. Since that time, nu- 
merous contributions have been made to the theory of vibration of elastic ele- 
ments using the Rayleigh-Ritz method. 

Gate Skin Plates.—The vibration of the skin plate on a submergible gate at 
Cheatham Dam on the Cumberland River offers an example of this type of vi- 
bratory motion. The operators reported a loud noise issuing from the Cheat- 
ham gate when the overflow head was about 3 ft. The gate is shown in Fig. 11. 

The Waterways Experiment Station investigated the cause of the noise in 
coeperation with the Nashville District. Accelerometers were placed to meas- 
ure the principal freedoms of motion of the gate itself as well as those of the 
upstream and downstream skin plate. The specific source of the noise was not 
known at that time. 

It was found-that the frequency of the downstream skin plate approximated 
the frequency of the sound. The frequencies of otherfreedoms of motion were 
substantially different. It was therefore concluded that the downstream skin 
plate was vibrating and causing the sound. 

The resonant air chamber inside the gate was analyzed but found to have a 
low frequency in comparison. Similarly, the vortex trail from the strut arms 
possesseda lowfrequency. As the exciting force of the vibratory motion could 
not be readily determined, itwas recommended that the skin plate be restrained 
bya support in the middle of the panel. It was reasonedthat this would greatly 
increase the natural frequency of the skin plate and eliminate the vibration which 


emitted the noise. 


EXCITING FORCES 


Self Excited Oscillations.—The action of the bow ona violin string is the 
classic example often cited. A good example of self-excitation from the hy- 
draulic engineering viewpoint is the interaction of a reflecting pressure wave 
a conduit and the vertical vibration of a gate. This is shown schematically in 
Fig. 12. 

me small vertical displacement of the gate would cause a pressure wave to 
travel upstream in the tunnel. When the wave reaches the reservoir, it is re- 
flected and changes from a positive pressure wave to a negative pressure wave 


or vice versa. 


19 «Gesammelte Werke,” by Walther Ritz, Paris, 1911. 
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If the speedof the wave is a and the distance from gate to reservoiris L, the 
period of the pulsation against the gate bottom wouldbe (4L/a sec. The pres- 
sure pulses on the gate bottom are shown on Fig. 12(b) as a square wave for 
the sake of simplicity. When the natural period of vertical vibration of the gate 
is close to the natural period of the pressure wave pulse, a condition of self- 
excitation could exist. 

It is interesting to note that the Von Karman vortex trail is a self excited 
oscillation in itself. 

Free Water Surface Phenomena.—Two free water surface phenomena are 
shown in Fig. 13. Either could constitute an exciting force on an elastic struc- 
tural element. 

The fluctuating nappe for a small head on a sloping weir has been noted by 
various observers. This phenomena was extensively observed and analyzed by 
Bruno Leo. The phenomena was also observed by the Bureau of Reclama- 
tion21 on the drum gates at Black Canyon Dam. In this case, the vibration was 
eliminated by aerating the space under the nappe. 
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FIG, 12,—SELF-EXCITATION WITH FIG, 13.—EXCITING FORCES— 
PRESSURE WAVE FREE WATER SUR- 
FACE 


Channel Waves.—The binodal wave which has a clapotis type actionfrom 
each side wall was observed in the fish ladders at McNary Dam. A trinodal 
wave was also observed in the Bonneville fish ladders. These waves can con- 
stitute an exciting force on adjoining structures. 

Other Exciting Forces.—Various other phenomena can possibly constitute 
exciting forces in vibratory motion. The successive formation and collapse of 
vapor cavities in the phenomenon of cavitation seem to have a periodicity. 
Little experimental information is available for the multitude of possible geo- 
metrical boundary situations. Some experimenters have suggested a Strouhal_ 
type number which defines the frequency of shedding of cavities. 


20 “Self Excited Vibrations at Overflow Weirs (Selbst 
gestenerte Schwengungen an | 
Uber Stromten Wehren),” by Bruno Leo, U. S, Waterways Experiment Translation No. 
boss, Vicksburg, Miss., June, 1955, 
21 “Report on Vibration Studies Made at Black Canyon Dam, by R, E, Glover, C. W. 


Thomas, and T, F, Hammett, U. S, Bur, of Reclam,, Hydr. Lab. Report N 
Colo., July 24, 1939, , ah a aaa 
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The fixed cavity which forms on the trailing side of an obstruction in high 
velocity flow has an intermittent pressure pulse. However, little research has 
been conducted on this problem. 
The toe of the hydraulic jump is known to pulse in an up and downstream 
| direction. The possible effect of this is unknown. 


CONCLUSIONS 


1. The modulus of elasticity of wire rope suspensions can be determined 
by measurement of vertical vibrations. 
: 2. The type of friction whichaffects the damping of hydraulic gates has been 
_ shown to be constant friction or Coulomb damping. 
3. The magnitude of the constant frictionforce can be evaluated by ananal- 
ysis of the oscillograms. 
4. The Von Karman vortex trail has been observed to be the excitingforce 
which can cause vibration of gates. 
5. A gate leaf with a sloping bottom sheds the vortes trail from the down- 
stream edge and therefore minimizes vibratory motion. 
6. Extensive researchis neededon the magnitude of the exciting forces be- 
fore gates can be designed against vibration. 
7. Considerable research is needed on the character of vibration of both fixed 
and traveling cavities and other hydraulic pulsating phenomena. 
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STREAM-GAGING NETWORK IN THE UNITED STATES®@;P 


By John E. McCall,! F. ASCE 


SYNOPSIS 


The national stream-gaging program is analyzed with regardto size, cost, 
comparison with other countries, and between states, Its deficiencies are pin- 
pointed and conclusions reached as to changes needed. A modified program 
of areal stream gaging is described that wouldovercome some of the existing 
deficiencies. Progress made toward adapting such a program is reported and 
remaining problems are outlined. 


INTRODUCTION 


At least four thorough studies2,3,4,5 have been made since 1949 to evaluate 
the adequacy of hydrologic data in this country. The strengths and weaknesses 
of existing data programs were analyzed at great length and a number of rec- 
ommendations were made for needed improvements. Three of these studies 


Note.—Discussion open until August 1, 1961, To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961. 
2 publication authorized by the Director, United States Geological Survey, Dept. of the 
Interior, Washington, D.C. 
Presented at the March 1960, ASCE Convention in New Orleans, La. 
1 Dist. Engr., U. S. Geological Survey, Trenton, N, J. : 
« Adequacy of Hydrologic Data to Meet Federal Needs,” Federal Inter-Agency River 
Bagin Committee, 1949, 
“A Water Policy for the American People,” President’s Water Resources Policy 
Comm,, U. S. Govt, Printing Office, Washington, D. C., 1950. 
“Water Policy,” Presidential Advisory Committee on Water Resources Policy, U. 
35 . Printing Office, Washington, D. C., 1955. 
: ope icnecrisd Appraisal Ee escen Data,” Task Group on Hydrologic Data of 
he Committee on Hydrology of the Hydr. Div., Proceedings, ASCE, Vol. 85, No. HY 7, 
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also attempted to estimate the required size of program or number of stations 
needed to provide adequate hydrologic data of the several principal types. The 
estimated number of stream-gaging stations needed, as determined in these 
three surveys, is shown on Fig. 1. Also plotted for comparison is the actual 
number of gaging stations, showing the growth rate since 1903. 

It can be seen that all three of the surveys concluded that the number of 
gaging stations should be expanded rapidly. The surveys differed only slightly 
as to the total number of stations recommended by given target dates. Never- 
theless, the growth rate of total stream-gaging stations slowed down about 
1952 in the face of these carefully considered recommendations. To under- 
stand this paradox one must look at the existing program from several different 
perspectives and then follow the evolutionof new concepts aimed at getting the 
maximum amount of hydrologic information for each dollar invested in stream 
gaging. 


THE PROGRAM FROM SEVERAL PERSPECTIVES 


Size and Cost of Program.—Operation costs for a stream-gaging station 
in the United States range from $600 per year to as much as $5,000 a year 
for afew remote stations, with an average of about $950 per station per year 
in 1960. Installation costs range from $1,000 to $15,000, and average about 
$3,000. The average life of a structure is on the order of 30 yr, making de- 
preciation, without interest, about $100 per year. 

The cost to the nation averages $950 for operation plus $100 for amorti- 
zation of capital investment for each of about 7,150 stations, or a total of $7.5 
million per year spent by the United States Geological Survey, Dept. of the 
Interior (USGS) alone. Of this, about one-third is contributed by the states 
on a 50-50 matching basis and nearly one-fourth is transferred to the USGS 
from other Federal agencies, Perhaps another million is spent annually by 
other Federal agencies within their own organization for stations closely re- 
lated to internal operations, 

The Mississippi River Commission annually publishes about 60 records of 
river stage and discharge and another 100 records of river stage only.® The 
International Boundary and Water Commission, United States and Mexico, 
publishes approximately fifty streamflow records collected in the Rio Grande 
basin." The USGS also publishes some records collected by other agencies. 

The total national investment, including state, county, andcity expenditures 
made in cooperation with the USGS, is about $9 million each year for basic 
streamflow data. This amount represents less than 0.5% of the amount spent 
in the country annually for dams, levees, irrigation works, water supply 
sources, and other works dependent on streamflow datafor safe and economi- 
cal design,8 Expressed another way, it costs about five cents per person 
per year for operation of streamflow stations at all levels of government, This: 
expenditure is modest considering the extent to which man is dependent on 
water for luxuries, necessities, and indeed life itself. Yet continuing infla- 


OE 
6 “Annual Report of Stages and Discharges of MississippiRiver and its Tributaries,” 
Mississippi River Comm., Vicksburg, Miss. i 

7 @ annual Report of International Boundary and Water Commission, United States and 
Mexico,” El Paso, Texas, ’ 


8 «Water Facts for the Nation’s Future,” by W. B. Langbein and W. G. Hoyt, Ronall 
Press, N, Y., 1959. 
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tionary trends and mounting Federal budgets make it imperative to economize 
even more if that is feasible. 

Comparison with Other Countries.—How does the United States compare 
with other countries? Fig. 2 shows the number of gaging stations per thou- 
sand square miles in relation to total area for various countries, as listed by 
Walter B. Langbein.? It is evident that area alone is not a reliable guide to 
the number or density of gaging stations needed. The range is infinite. Libya 
has no gaging stations in an area of 633,000 sq miles, while Israel has 85 
stations in an area of 8,000 sq miles. It might be inferred from this plot that 
the smaller countries have, of necessity, developed their natural resources 
more intensively than the larger countries, and that gaging-station density 
reflects this in a general way. Note that the United States is considerably 
higher in gaging station density than any other country of comparable size. 

Gaging-station density in relation to population density is illustrated in 
Fig. 3. Including population as a parameter along with area makes the re- 
lationship take a vague form that can begin to be interpreted. Note that un- 
developed countries plot low and usually toward the left. The correlation is 
still too scattered to be very useful. However, on the basis of population, the 
United States does not seem unexpectedly high or low in gaging-station density. 
As the intensity of development continues to rise and as the “population ex- 
plosion” will presumably continue, the position of the United States might be 
expected to continue to rise and drift further to the right. The change in posi- 
tion since 1930 is shown by the dashed line. Of course, most other countries 
are tending to move in the same general direction at a rate dependent on the 
speed of development of their water resources. 

Comparison of States Within the United States.—Perhaps better evidence 
that the number of gaging stations is related to population density can be 
seen in Fig. 4, as reported by Langbein.? This graph shows the relative posi- 
tion of each state in continental United States as of 1955, after empirical ad- 
justment for irrigation and waterpower. The adjustments were determined by 
trial, their effectiveness being measured by reduction in scatter of the plotted 
points. Langbein’s analysis indicated that, insofar as stream-gaging needs 
are concerned, each acre of irrigated farmland was equivalent to 4 people. 
Likewise, each hydro-power plant of 1,000 kw or more capacity was found to 
be equivalent to 10,000 people. With these adjustments the points tend to de- 
fine a rather consistent relationship between number of gaging stations and 
population density. That is to say, that in the United States the stream-gaging 
program in all states is affected by the same national policies and that the 
standard of living or stage of economic development does not vary appreciably 
between states. New Jersey with 10.5 gaging stations per thousand square 
miles and Nevada with 0.5 station per thousand square miles represent nearly 
the two extremes in both population density (2nd and 48th) and industrialization 
in conterminous United States. All other states line up reasonably well on 
this plot. 

Effect of Administrative Policies.—It is interesting to note that Alaska, 
which has replaced Nevada as our least developed and populous state since 
Langbein’s analysis (1957) would plot off the sheet to the lower left but about 
in line with the other states, Alaskaisthe only state that has never cooperated 


9 «Numbers of Stream-Gaging Stations in Various Countries with Analysis of their 
Distribution in the United States,” by W. B, Langbein, Bulletin de 1’ Association Interna- 
tional d’Hydrologie Scientifique, No, 8, December, 1957. 
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on a 50-50 financial basis, and all gaging stations in Alaska prior to 1959 
were financed completely by the Federal government. Hawaii would plot off 
the top of the sheet, having twice the gaging-station density of New Jersey, 
the previous leader in this respect, but only one-tenth the population density 
of New Jersey. While Hawaii is hardly typical, being made up of eight major 
and many smaller islands, it seems likely that the administrative policy of 
the Territorial government was a major factor in reaching such a high level 
in gaging-station density. The huge pineapple plantations probably exerted a 
major influence in this because of their irrigation and drainage needs. Hawaii, 
in contrast to Alaska, has for several decades furnished half of the stream- 
gaging dollars from local taxation. 

In 1929 the Congress first authorized the Federal matching of state or 
territorial funds for stream gaging. This principle of 50-50 matching has 
been followed ever since and the Federal-state cooperation has become the 
backbone and major source of funds for the work. Referring again to Fig. 1, 
the sharp rise in the curve since 1929 shows very plainly the effect of this 
partnership policy on growth of stream-gaging. 

Israel is another example of a country which, by administrative policy, 
decided it would have intensive stream gaging and ground water exploration. 
Taking advantage of all feasible means to “lift itself by the bootstraps” from 
a backward nation to a position of leadership in the Middle East, this tiny 
country has perhaps the most intensive and best planned hydrologic investi- 
gations’ program of any nation in the world (as of 1960). 

Preponderance on Large Streams.—The early pressures for gaging stations 
were for navigation, hydroelectric power, and flood control with the result 
that practically all older stations are on the larger streams. This had the 
advantage that approximately 80% to 90% of the runoff came under measure- 
ment rather quickly. The only surface runoff not being measured today in 
conterminous United States is in relatively small coastal streams and in tidal 
reaches of the larger rivers. In fact, some water originating in the upper 
Missouri or Ohio Rivers may be measured a dozen or more times before it. 
reaches the ocean, | 

However, in most instances this flowing water would not be measured until 7 
it had passed through a runand abrook, or perhaps a creek or two, and finally 
entered a river draining several hundred square miles. Less than one-third | 
of the gaging stations in the United States are on watersheds of 100 sq miles | 
or smaller. Less than 5% of all gaging stations are on drainage areas of 10° 
Sq miles or less despite a three-fold increase in small-area gaging-station 
coverage since 1945. Yet there are many, many more miles in aggregate 
length of small streams than there are total miles of the large streams. As 
the years pass, the smaller streams are becoming of ever greater economia 


or scientist interpolate the characteristics of the stillungaged small streams, 
which probably will always outnumber the gaged streams a hundred-fold. 

is a large order because the flow characteristics become more diverse 
variable as the drainage area becomes smaller. 
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It is interesting to ponder the question of the spacing of gaging stations, 

Assume the choice is one station for every 20, 50, 100, 200, 400, or 700 miles 
bf stream, For small streams, a choice of one every 50 miles would be the 
same aS Saying there should be one on every tenth stream averaging 5 miles 
in length. Most engineers quizzed by the writer intuitively chose from 50 to 
100 miles as a reasonable distance between stations, yet even the more con- 
servative of these two choices would require quadrupling the size of the present 
network, 
From spot sampling of stream density (map studies of perennial stream 
length divided by area) in various parts of the country, E, W. Coffay10 has 
jound the average to be slightly more than one mile length per square mile 
of area. Although this study has not been carried to completion (as of 1961), 
und the final answer may be slightly different, a rounded figure of 3 million 
miles of perennial streams in conterminous United States can be used safely 
jor estimating purposes. Deducting stations in Alaska and Hawaii but including 
ull other streamflow stations operated by the USGS and other Federal or 
state agencies (6680), the ratio comes out to one station now in existence for 
every 450 miles of streams. 

There is no large stream in the country with gaging stations as far as 450 
miles apart, on the average, so it must be concluded that stations are still 
larther apart on the small streams at present, In fact, stations now may aver- 
wee less thanone for every 600 miles of aggregate length on the small streams. 
Admittedly, a mileage basis is no way to design a stream-gaging network; its 
mly significance is to lend perspective as to the magnitude of the task of 
retting adequate data on all small streams, 

Other Segments of the Hydrologic Cycle.—An excellent and up-to-date 
ippraisal of the six principal types of hydrologic data was given in 19599 by 
he Task Group of Hydrologic Data of the Committee on Hydrology, Hydraulics 
Nivision, American Society of Civil Engineers (ASCE). It will serve here 
simply to quote excerpts of the Task Group’s conclusions which are pertinent 
0 this appraisal. 

The Task Group found that: 


“The existing federal networks of gages measuring hydrologic data 
are producing satisfactory records. Periodic reappraisals are desirable 
for bringing about improvements to keep pace with changing needs and 
conditions. A reluctance to change prevailing practices can lead to over- 
collection and failure to begin records needed for the future.” 


The report noted several new types of data of current and future importance. 
"hese included natural or background radioactivity of waters and changes due 
9 radioactive poilutants, data to evaluate effects of watershed treatment 
neasures, and artificial rain-making. 

‘In an earlier part of the Task Group report devoted specifically to runoff 
t was stated that: 


“The most obvious deficiency is in data from small drainage areas 
of 50 square miles or less. More data are needed on streamflow in the 


10 personal Communication with the author. 
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zone of tidal influence, and such data already collected should be as- 
sembled, analyzed and published.” 


INFERENCES AND TENTATIVE CONCLUSIONS 


Looking at the national stream-gaging program from several different 
perspectives leads to the following tentative conclusions: 


1. Present costs are not excessive; five cents per person per year or 
one-half of one percent of annual expenditures for water resources develop- 
ment, 

2. The United States is not unexpectedly high or low in gaging-station 
density in relation to other countries. Considering the high standard of living 
and extent of industrialization in the United States, it may be lower than ex- 
pected. On the other hand, the United States has considerably more stations 
than any other country of comparable size. 

3. The need for streamflow data apparently is related to population density 
and extent of industrialization, both of which are projected to grow very 
rapidly in the United States in future decades. 

4. Greatest present and future need is for data on thousands of small 
streams, Past methods of building permanent gaging stations wherever data 
are needed will not suffice even if funds and manpower could be doubled or 
tripled. 

5. To keep pace with the widespread needs, the program must be directed 
toward obtaining less complete data but on a broader coverage—particularly 
on smaller streams and tidal streams. Also some new types of data should 
be observed and published to meet recently emerging requirements. 


The phrase “areal stream gaging” is used to describe such a program, in 
which all available streamflow, flood, and drought records are evaluated along 
with topographic, geologic, and meteorologic informationto determine regions 
or areas of hydrologic similarity. Within such an area the shape of the hydro- 
graph is reasonably similar for streams of equal drainage area; peak runcff 
per square mile, average runoff per square mile, and base flow character- 
istics tend to be comparable, Within such an area only one or two complete 
gaging stations need be maintained continuously, with variations from place 
to place within the area being defined by reconnaissance methods. 


REVIEW OF THE NETWORK 


Development of Concepts.—The thinking and concepts previously expressed 
were actually the development of many years of careful consideration and 
study by many hydraulic engineers. The first person to document these ideas 
formally, as far as the writer is aware, was Walter Langbein in 1951. The 
pros and cons of these concepts were debated extensively within the USGS in 
the months that followed, and in 1952 a committee of six engineers was ap-. 
pointed to study this matter formally. This committee, of which Langbein was 
chairman, submitted its report in 1953. The report was reproduced and dis= 
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[tributed within the USGS for study, comment, and field trial of the ideas and 
procedures proposed. The concepts of the committee were stated as follows: 


“The principles that are proposed as guides inthe development of the 
gaging-station program in order to achieve maximum information rela- 
tive to costs are the following: (1) Apartfrom operational needs, stream 
gaging is a sampling process, and (2) within limits it is possible to cor- 
relate the flow of one stream with that of another, A gaging station pro- 
vides information on the rates and amounts of flow during a particular 
time in a particular channel and therefore samples a given time and 
given place. The variations that occur with respect to time at a gaging 
station are fairly well understood, but those that exist from piace to 
place less so. A well-balanced network, therefore, gives the two factors 
of time and place their proper consideration. The two factors can be 
brought into balanced relationship through use of correlation techniques. 


“These basic principles lead to the conclusion that a record can be 
viewed not only as a measure of the flow at one site but also as an in- 
dex to the flow of streams in the contigrous area—in effect magnifying 
its usefulness manyfold. These principles lead toward areal rather than 
spot stream gaging.” 


In November, 1955, Langbein andC.H. Hardisonpresented a paper!! which 
described the theory, procedures, and statistics involved inthe correlation of 
a short record with a nearby long-term gaging-station record. These tech- 
miques are very useful in applying the program concept of a hydrologic net- 
work of stations. The Presidential Advisory Committee on Water Resources 
Policy, in December, 1955, made recommendations* which added further 
impetus to this concept of a hydrologic network, made up of long-term and 
short-term stations. The USGS ran a pilot projectin the Colorado River basin 
to test the feasibility of the concepts and the problems involved in attempting 
to adapt an existing stream-gaging program to these network principles. The 
problems were many but none insurmountable. The results of the pilot project 
were deemed to warrant adopting the network principles as a guide for future 
programming. ; ; 

During the last half of 1956 and all of 1957, a formal and detailed review 
of the entire stream-gaging program in the United States and its territories 
was made by the USGS in accordance with these network principles. Since that 
time the USGS has been following these principles in its program design. 
Naturally, these concepts and the modification in program goals were dis- 
cussed with representatives of other Federal agencies in the Subcommittee on 
Hydrology of the Inter-Agency Committee on Water Resources. They were 
also discussed with local cooperating officials in state and federal agencies 
at different times during the consideration, testing, and adaptation of these 
ee ane during this period of more than 8 yr that these concepts have 
yeen developing have engineers of the USGS felt that they have known all the 
answers, Accordingly, there has been no great rush to get into print any final 
version of this type of network design. However, it is felt that the concepts 
lave been sufficiently tested and that enough answers are known to warrant 


. . . 
11 “Extending Streamflow Data,” by W. B, Langbein and C, H. Hardison, Proceedings, 
ASCE, Vol. 81, Proceedings Separate No, 826, November, 1955, 
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publication. Furthermore, it has become apparent that lack of wide dissemi- 
nation of the ideas and concepts has somewhat deterred their acceptance by 
the engineering profession and officials concerned with streamflow data pro- 
grams, 

Additional information has been presented by Langbein and W. G. Hoyt. 
This book, sponsored by the Conservation Foundation, contains a compre- 
hensive discussion of hydrologic data networks and programs. Therefore, 
only general principles, procedures, and results will be mentioned here. 


TABLE 1.—CLASSIFICATION OF SURFACE WATER RECORDS 


Streamflow Station Partial Record Station 


Hydrologic Network Hydrologic Network Hydrologic Network 
Primary (long term) Primary (long term) Primary (long term) 
Bil Areal (Ponds & lakes) None 


All Areal (including ne- 
cessary sup- 
porting records 


A12 Mainstream 


B12 


Mainstream (Rivers 
& tidal 
estuaries 


Secondary (short term) 
A 21 Areal 


A 22 Mainstream 


Secondary (short term Secondary (short term 
B 21 Areal (Ponds &lakes) | C 21 Flood crest 
B 22 Mainstream (Rivers & | C 22 Low flow 
tidal oa=-1 oy paid a Takes 
tuaries) bat e inventory 
C 24 Periodic streamflow 


A 23 Seasonal daily - 


hydrologic 
Water Management Water Management 
Long term_ Long term 
Compact 
A 32 Legal 
A33 ional 
A 34 Administrative Administrative 
A 35 Basin Accounting Basin Accounting ; 
36 Federal Power Comm.,|B 36 Federal Power Comm. ; 
Short term Short term Short term | 
41 Research & experi- B 41 Research & experi- C 41 Flood hydrographs nad 
mentation mentation and timing ; 
Detailed des B 42 Detailed design . 


B 43 Operational C 42 Seepage and low flow — 


> 


Classification of Stations.— The first step in the actual modification of the 
program was to review and classify all existing gaging stations as to their 
best use in a nationwide network. Table 1 shows the classification system 
adopted. Each main category was further divided into hydrologic network or 
water management stations, A water management station is defined as one 
needed to provide data for one or more specific operational, legal, or admini- 
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strative purposes. Its establishment, location, and length of operation depend 
almost entirely upon the needs of the project or projects for which it is oper- 
ated. Accordingly, it was decided at the outset that this type of station could 
hot be programmed in advance by consideration of general basin hydrology 
and general purpose needs. The programming problem was thereby narrowed 
own immediately to the design of a general purpose hydrologic network. 
However, in all cases for which a gaging station was found in the original 
review to serve hydrologic network purposes as well as certain water man- 
ugement purposes, the station was classified as a hydrologic network station 
‘(or economy reasons, 

_ The hydrologic network was broken down into primary stations and sec- 
mndary stations. The primary stations are the long-term stations to be oper- 
ited for an indefinite period of time. Secondary stations are intended to be 
of short-term operation, say from 5 yr to 10 yr. At the end of that time the 
station, and as much of the equipment as can be moved, would be relocated to 
some other point needing stream gaging. Both the primary and secondary 
slassifications were further broken down into areal and mainstream stations. 

Primary Stations.—Areal primary stations are index stations on streams 
hat must essentially be free from past regulation and, hopefully, free from 
xxtensive future regulation, diversion, or other development. In parts of the 
United States it is becoming increasingly difficult to find stations completely 
ree from regulation and diversion. This made the job more difficult but no- 
where prevented the use of the network principles. Thus an areal primary 
station is selected for representativeness and length of record; insofar as 
Ossible it is free from past and future regulation or diversion, and will be 
yperated for an indefinite period of time to obtain a long range time sample 
wf the hydrology of the section in which it is located. This station record would 
ye used as the independent variable in making correlative estimates of long- 
erm streamflow characteristics at other sites in the same hydrologic pro- 
vince, 

A mainstream primary station, on the other hand, may be considerably 
iffected by regulation or diversion. It serves only as as a record of actual 
low at that point and as an index of flow at other points upstream or down- 
stream on the same large river. An example of this would be the Mississippi 
iver at Vicksburg, at which point the river flow would correlate very closely 
vith simultaneous flow in the river at other points upstream or downstream 
or many miles, but would not bear any particular relationship to the flow in 
mall tributaries entering the Mississippi River in this vicinity. 

The areal primary stations selected were those which seemed to be most 
‘epresentative of the hydrologic characteristics of the area in which they 
yere located. Trial correlations of monthly mean flows for a 5-yr period, 
‘enerally 1951-55, were used for this purpose. It was found that stations with 
ither very large or very small drainage areas did not correlate well with 
urrounding stations. Therefore, most areal primary stations, as a result, 
re stations with medium size drainage area. Attention was given to this 
joint to insure that some primary stations were selected on smaller basins 
nd some on larger drainage areas in all parts of the country in order to 
void complete bias in the final network. 

Secondary Stations.—An areal secondary station would be operated where 
eneral streamflow information is wanted or will likely be needed in the 
uture, the length of the required record depending on the number of years 


@ 


: 


90 March, 1961 HY 2 


necessary to define a correlation with a nearby areal primary station ade- 
quately. A mainstream secondary station, likewise, would be operated for 
only aS many years as needed to define a correlation with a nearby main- 
stream primary station. Then either type of secondary station might be dis- 
continued or converted to a partial-record station. Because of their frequent 
relocation, the design of new secondary stations will feature prefabricated, 
demountable, and semi-portable type shelters that can be quickly knocked 
down and moved. 

Partial-Record Stations.— The principal emphasis on partial-records ison 
flood peak flows and base flows, the two extremes of river regimen, Flood 
peak partial-records are frequently obtained by use of the “crest-stage gage.” 
A small vertical pipe on a tree or bridge abutment has inlet holes near the 
bottom and is charged with a spoonful of granulated cork. After a flood the 
cork leaves a high-water mark on agraduated stick left in the pipe. The engi- 
neer visits the station periodically, preferably after each large rise, and notes 
the peak stage since the previous visit, removes the mark, and re-charges 
with cork. The cost of operating such a station is from 5% to 30% of a com- 
plete gaging station cost, depending on need for and difficulty of calibration 
of the site to develop a rating curve. 

Base flow measurements are made several days to a week or more after 
the most recent rain on an area. The flow includes no direct storm runoff and 
is essentially all from ground-water discharge. Base flow is systematically 
measured in the spring, when the ground-water table is ordinarily at its 
maximum elevation for the year, and then measured again in the late summer 
or fall, when ground-water levels and discharges to the streams are usually 
at or near the minimum for the year. Thus, two widely separated points on 
the correlation graph are obtained each year for several years to define the 
base flow at any particular site in comparison to a nearby secondary or pri- 
mary gaging station, The cost is about 5% to 10% of the cost of a complete 
gaging station. 

Streamflow Correlations .—The correlation technique used in the nationwide 
review of the stream-gaging network was the comparison of simultaneous 
monthly mean flows. A graphical solution was usually made, as shown in Fig. 
5. At least 60 months (5 yr) were plotted and a smooth mean curve drawn 
through these points, The standard error was estimated by drawing two curves 
parallel to the mean curve at suchdistancesfrom it as to exclude the top one- 
sixth and bottom one-sixth of the plotted points or months, thus including two- 
thirds of the points, or the part included within plus and minus one standard 
deviation, In this example, two-thirds of the months were between +22% and 
-18% of the mean curve. 

In general, the correlations of streamflow characteristics were found to 
be good in the humid east and southeast. Correlations were good to fair in 
the Midwest and northwestern parts of the country, and fair to poor in the 
arid southwest and states immediately east of the Rocky Mountains. 

Other possibilities for making correlations between primary and secondary 
stations or partial-record stations are available. Concurrent daily flow, 
during base flow conditions, is especially useful for low flow partial-record 
stations, Discharge for equal percent duration or equal recurrence interval 
could also be used instead of monthly mean flows to extent records at a 
short-term secondary station. The extremes of flow at the short term station 
may also be estimated by a relation developed with the long term station but 
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Discharge, in cubic feet per second; Sepulga River near McKenzie, Ala. Drainage area 470 square miles 


FIG. 5. CORRELATION OF SIMULTANEOUS MONTHLY MEAN FLOWS 
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it is usually preferable, especially where great extrapolation of the relation 
would be required, to make a study of flood frequency or drought frequency 
by other methods. There is the additional prospect that rainfall data can be 
used to improve the correlation between the flows of two streams, Multiple 
correlation using records for several nearby gaging stations may sometimes, 
but not always, improve the results, 12 

Scope of the Review.—Insofar as it was practicable, the USGS’ review of 
the network of gaging stations included stations operated by states or other 
agencies such as the Mississippi River Commission, the International 
Boundary and Water Commission, United States and Mexico, the Corps of 
Engineers, the Bureau of Reclamation, and agencies of the Department of 
Agriculture. The criteria for inclusion of stations operated by other agencies 
or states required that the station records be published regularly and that 
the station be suitable for inclusion in the hydrologic network. If a station 
record collected by another agency was considered accurate and readily 
available to potential users, it was felt that there would be no need for dupli- 
cation with an additional hydrologic network station operated by the USGS. 
Stations, operated by other agencies, that would be of use for water manage- 
ment purposes only, however, were not classified. These policies yielded 
coverage that is somewhat less than perfect for the nation as a whole. A con- 
siderable number of gaging stations are operated by other agencies in con- 
nection with research or operational projects, but the basic data have not 
been published though the evaluated results may have been included in special 
reports. Nevertheless, probably in excess of 95% of all existing stations in 
the country were included in the review. So the results are considered valid 
and useful from the standpoint of national coverage. 

Results of the Review.—Certain aspects of the review were quite direct 
and yielded definite results. For example, it was possible to assign a classi- 
fication of best use in the nationwide program to every existing gaging station. 
As might have been expected, even this was not timeless and there have been 
a number of reclassifications since the initial review was completed. However, 
the number of reclassifications has been very small compared to the number 
of stations, less than 1% in 3 yr, so the classification of existing stations 
may be considered fairly firm. 

Another firm finding was that some 700 specific basins in this country did 
not have a single gaging station within the area that met areal primary station 
requirements, Most of these basins are inthe Great Plains and Rocky Mountain 
areas where existing station density is low and variability of hydrologic 
characteristics greater than in the humid areas. The correction of this de- 
ficiency is being given high priority. 

The review showed that more than 1,200 existing secondary stations could 
be moved to other locations after a period of operation of 5 yr to 10 yr, and 
that a good share of these were ready for immediate relocation, It was further 
shown by brief correlation tests that about 1,150 gaging stations previously 
discontinued had accumulated sufficient length of good records to permit ade- 
quate correlation with primary stations, It was also found that some 1,250 
existing stations classified in the water management categories would qualify 
for the secondary hydrologic network (except for relocation after 10 yr). Thus, 
a total of about 3,600 points have been or are being gaged in the secondary 
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12 “Graphical Correlation of Gaging-Station Records,” by J. K. Searcy, U.S. Geol, 
Survey, Water-Supply Paper 1541-C, 1960, i q 
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hydrologic network even though only 1,200 stations are currently classified 
‘as secondary stations. 

A result, much less firm than the preceding ones, was the tentative con- 

‘¢lusion that about 2,500 additional short-term secondary stations would be 
meeded in the next ten years. This conclusion was reached by adding up the 
specific sites which, in the judgment of the 44 district engineers of the USGS 
‘across the country, would be needed in their districts. Thus, it is strictly a 
‘matter of judgment and cannot be proved right or wrong although it is the best 
‘guide available. 
__ The last, and perhaps most hypothetical, result of the review was the indi- 
cation that a total of about 10,000 active gaging stations in the United States 
is the maximum number that would ever be needed and justifiable at any one 
time. This figure was arrived at by adding the 700 additional primary and 
2,500 additional secondary stations needed, as previously described, to the 
existing network and assuming that for every station established thereafter 
another one would be dropped. The figure of 10,000 assumes a slow but steady 
increase in water management stations and a very large increase in partial- 
record stations for reconnaissance type data. No specific date was set as a 
‘target for obtaining 10,000 stations. A straight line projection on Fig. 1 indi- 
cates reaching this level about 1975-85, much later than the goals established 
in the three previous surveys. 

A side-effect, if not a result, of the review was the strengthening of the 
‘conviction of those participating in the study that a broadening of the types of 
data collected would be more valuable than simply adding more years of daily 
discharge record at all existing gaging stations. Data on flood plain inundation, 
water use, time of travel, flood profiles, bankful stage, water temperature, 
and chemical quality (to name only a few), are becoming more important and 
necessary every passing year. 

Trends Since the Review.—Fig. 6 shows the number of streamflow stations 
in each classification as of January, 1957, andJuly, 1959. Of the 700 additional 
areal primary stations the review found to be needed, about 80 had been es- 
tablished by July, 1959. Only a net of 36 additional areal secondary stations 
were in operation but nearly 350 new secondary stations had been established 
during the period. Similarly, for all classifications of streamflow stations, 
data were being collected at 660 new locations, but because of selective dis- 
continuation the net increase was only 115 stations, 

_ There were very few changes in number of stage stations during this 
period, 34 being established and 17 discontinued for a net increase of 17. 
However, the expansion in the number of partial-record stations during the 
period was particularly encouraging. The exact number in existence in 1957 
is not known but is believed to have been about 500 low-flow and 800 crest- 
stage stations. The count in September, 1959, shows growth to more than 
2,100 low-flow stations and about 1,950 crest-stage gages. 

At the risk of misleading the casual reader, the increases in the several 
categories and classifications can be totalled, though strictly speaking each 
is different from all others. It can be generalized that surface-water data 
are being collected at nearly 3,500 new points since January, 1957; the data 
are less detailed and require analysis for full use, but the cost in funds and 
manpower has risen only about the equivalent of 600 new complete gaging 
stations. 
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Remaining Problems.—Although the concept of a hydrologic network seems 
to be working out very well, not only in theory but in actual practice, there 
are quite afew remaining problems. The most immediate, if not the largest, 
of these problems is that of gaining acceptance of this type of program design 
by engineers and hydrologists throughout the country. So far there has been 
no outright disagreement with the concepts or opposition to the policies from 
anyone close enough to the problem to study it through. However, the ideas 
and the procedures have not been well advertised and thoroughly discussed 
in the literature as yet, to inform fully all who need to know about this ap- 
proach, Natural resistance to change until the reason for change is fully 
understood has arisen in a few isolated instances. These can probably be re- 
garded as the exceptions that prove the rule. Most of these instances arise 
when it is proposed to discontinue a secondary station which is no longer 
needed for hydrologic network purposes, Better understanding of concepts 
and program goals will gradually overcome this problem. 

Another fairly urgent problem appears to be the acquisition of funds to 
establish “bench mark” areal primary stations found to be needed. A “bench 
mark” station is one in a nationalora state park, a wilderness area, or other 
area where hydrologic change due to man’s activity does not appear likely. 
There are few existing stations where the hydrology of the basin will not 
change perceptibly in the next one hundred or more years. Wherever a basin 
can be found in which long-term, man-made change in the hydrologic regimen 
seems unlikely, an areal primary gaging station should be established to de- 
tect and document long-term trends as between wet and dry cycles. If an 
average of two or three such stations could be found in each state, they would 
represent a wonderful long-term investment in better understanding of 
regional hydrology and more reliable determination of the probable frequency 
of recurrence of given floods or droughts. 

Another problem arising as a result of this type of programming is getting 
all basic data assembled and published for wide use. Data from partial-record 
stations and other reconnaissance-type information must be either fully eval- 
uated and interpreted in published reports or the partial-record data and re- 
connaissance data themselves must be published. It is planned to publish data 
from crest-stage gages and base-flow stations along with continuous station 
records in water-supply papers of the USGS. 

The final problem which needs mention is that of getting underway more 
special studies and analyses in which all available hydrologic data are as- 
sembled into interpretive reports. These reports should contain maps, graphs, 
and tables to make the great mass of accumulated hydrologic data more readily 
useful to practicing engineers and other data users. In addition, shorter and 
less technical reports should be prepared directed toward school children and 
general public readers. These should contain the more significant facts to 
promote wise public use, control, and conservation of the nation’s water re- 
sources, The investment in data collection will be partly wasted if the knowl- 
edge gained is not fully disseminated. ‘ 


SUMMARY 
The stream-gaging program was reviewed from several different angles 
and found to be generally adequate. To overcome existing deficiencies on 
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a streams, and to orient the program better toward future needs, certain 
hanges were found to be desirable. Accordingly, a flexible system has been 
.dopted with a nearly permanent primary network of index stations and a 
-ontinuously changing secondary network of short-term stations. Both net- 
vorks are supplemented by a large number of partial-record stations de- 
higned to give specific data on flood peaks or low flow. 

The whole system is mutually supporting and correlations between primary, 
econdary, and partial-record stations will make data collected at each much 
more useful and valuable than if used alone. The network system provides a 
pgical guide for efficient expansion of the program and, if need be, for re- 
uctions with minimum damage. Above all, the program is designed to speed 

greatly reconnaissance-type coverage on small streams and to give maxi- 

um information for the manpower and funds expended. The growth rate, in 
rms of total stream-gaging stations being operated, has slowed down as a 
esult of network concepts. However the growth rate, in terms of total sig- 
ificant streamflow data being collected, has speeded up markedly. 
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STILLING BASIN DAMAGE AT CHIEF JOSEPH DAM4 
' 
; 


By R. H. Gedney,! M. ASCE 


| 
SYNOPSIS 


Chief Joseph Dam on the Columbia River, in the state of Washington, is 
he first dam downstream of Grand Coulee Dam. The dam, approximately 
“00 ft high, is a concrete gravity structure with a tainter-gate controlled 
| spillway, that discharges into a hydraulic jump-type, energy dissipator 
aving a single row of baffles and an end sill. Damage of the concrete in the 
sasin was discovered in 1957, two years after the project became operative. 
[The results of condition surveys made in 1957 and in 1960 are presented 
ogether with an evaluation of conditions that are considered to have contri- 
uted to the stilling basin damage. 


INTRODUCTION 


The occurrence of damage to spillway energy dissipators is not an un- 
Ommon event. However, there is seldom an opportunity to determine wheth- 
r the events causing this damage are attributable to construction or to 
perating conditions and to define the specific situations responsible. The 
epair of a spillway energy dissipating structure for a major dam after a 
roject is placed in service is often costly. Therefore, a full understanding 


Note.—Discussion open until August 1, 1961, To extend the closing date one month, 
written request must be filed with the Executive Secretary, ASCE. This paper is part 
f the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
iety of Civil Engineers, Vol. 87, No. HY 2, March, 1961. 

2 Presented at Seattle, Washington Meeting, Hydraulic Division Conference, August 
, 1960, ' 
1 Chf., Planning Sect., Corps of Engrs., Dept. of the Army, Seattle District, Seattle, 
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of circumstances leading to the occurrences of damage in Chief Joseph Dam 
stilling basin, that is a conventional structure constructed in accordance with 
-onventional procedures, should be helpful. It is interesting to note that dam- 
we to spillway energy dissipators in varying degrees has occurred at Grand 
SOulee, Chief Joseph, Rocky Reach, and Bonneville Dams, all of which are 
major structures on the Columbia River. 


PROJECT LOCATION AND SPILLWAY FEATURES 


Chief Joseph Dam is located 51 miles downstream from Grand Coulee 
Xam and backs water into the Coulee Dam tailrace. The dam primarily con- 
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FIG, 2,—CHIEF JOSEPH DAM, FLOW OVER COMPLETED SPILLWAY. 


ts of a concrete gravity gate-controlled spillway spanning the river chan- 
1, and an intake structure and powerhouse located just downstream, as 
own on Figs. 1 and 2. Use of the spillway is required annually, generally 
er a 5- to 7-month period during which time peak flows of about 350,000 
to 500,000 cfs can be expected. The historical flood on the Columbia 
ver at Chief Joseph site is about 730,000 cfs. The flood of record occurring 
1948 is 640,000 cfs, and the spillway design discharge is 1,250,000 cfs. 
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ormal head at the project during ordinary flood periods is in the order of 
agnitude of 150 ft from pool to tailwater and 200 ft from pool to stilling 
iSin floor. 

Flow over the spillway is controlled by nineteen tainter gates, 44.5 ft high 
nd 40 ft wide. During flood periods, the pool will ordinarily be maintained 
- Ele. 946 to 948 with a maximum rise to Ele. 957 for the spillway design 
ow. A paved concrete stilling basin that inducesa hydraulic jump for energy 
issipation is provided at the toe of the spillway. The stilling basin is con- 
scted to the ogee by a 100-ft-radius bucket section and is 167 ft long, from 
»int of tangency of the bucket to the end sill, as shown on Fig. 3. A single row 

baffles and a stepped end sill are located in the downstream portion of the 
illing basin. The stilling basin details and general plan are shown on Figs. 
and 5. The stilling basin floor is comprised of slabs designated in Fig. 5 
om Row D upstream to Row A downstream. 


TABLE 1.—HYDRAULIC JUMP DEPTHS 


Flow desi- 
gnation 
(2) 


villway? 1,250,000 0,90 
istorical aie 
maximum 730,000 6 


-Year 
‘requency 


400,000 


Pool at Elev, 957, all other pools at Elev. 946. 


, = Theoretical depth at toe of jump. 

1 = Theoretical velocity at toe of jump. 

2 = Theoretical required tailwater depth to form jump. 
9' = Actual tailwater depth available, 


The slabs have a minimum thickness of 5 ft and are tied to the rock 
undation with an extensive system of grouted anchor bars. Only the down- 
ream slab on which the baffles and end sills are located is reinforced; this 
ab is heavily reinforced to take the impact of high velocity flows on baffies 
d end sill. The floor of the basin is level at Ele. 743, about 12 ft below 
reambed, An extensively jointed granitic rock underlies the basin and is 
und at or near the surface of the streambed downstream of the stilling 
sin. The rock jointing varies greatly in strength and in many areas the 
ck is subject to erosion under high velocity flow. ee 
The stilling basin causes energy dissipation by hydraulic jump action. 
angle row of baffle piers, streamlined on the sides to reduce cavitation po- 
itial and a stepped end sill cause the jump to occur on the basin under max- 
um flow conditions and also aid materially in damping downstream mays 
tion. The jump conditions in the basin for spillway design, the ca : 
axximum flood and a 400,000-cfs flow (3-yr frequency), are given in Tal e : 
If allowance is made for the portion of the river flow that will be age e 
ough the powerhouse, the D2'/D2 ratio would be in excess of 1.0 for flows 
i 


rj 
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experienced as of 1960. The design features of the stilling basin were thorough- 
ly tested in model studies, using a 1:33 scale, 3-gate sectional model, and a 
1:80 scale general model that included the entire spillway. 

One of the outstanding features of the design is the coordinated effect of 
the baffle piers and end sill incausing formation of the jump. The baffle piers 
were located close enough to the end sillso that the end sill aids in creating a 
pressure rise in the vicinity of the baffles that minimizes the occurrence of 
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FIG, 4.—STILLING BASIN EROSION CONTOURS 


negative pressures. The basin was carefully tested in a 1:36 sectional mode 
to obtain an adequate length and depth and to avoid cavitation pressures or 
baffles and end sills. The model tests indicated an adequate margin of safe 
for satisfactory stilling basin performance for all flows and showed no in. 
dication of cavitation pressures for flows of less than 900 ,000 sec-ft. 
significant feature disclosed by the model tests was that a slight variation it 
the curve of the sides of the baffle piers from the theoretical shape specified 
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jesulted in high negative pressures at design flow not present when the exact 
hape was used, 


INITIAL REPORT OF DAMAGE TO SPILLWAY 
: Water was first diverted over a completed section of stilling basin in 
952. Thereafter, until the project was placed in service in 1955, a portion 
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FIG, 4, -CONTINUED 


f the stilling basin was utilized continuously each year under construction 

nditions. Continuous use of the spillway lasted until early in 1957 when 
ficient hydraulic capacity was installed in the powerhouse to permit com- 

e diversion of low water flows. Thus, initial inspection of the basin was 

layed until March 1957. : 

The first inspection of the basin, made over a two-day period, consisted 

imarily of four underwater traverses by a diver from one side of the basin 


104 March, 1961 HY 2 


to the other covering the end sill, the baffle piers and the paving from end 
sill to the spillway ogee section. This inspection revealed several large 
areas in which erosion had occurred to depths of zero to 12 in.; several ex- 
tensive areas in which erosion exceeding 12 in, had occurred and in which 
there was exposed reinforcing steel and marked damage to a number of 
baffles and several sections of end sill. The surface of the slab was eroded 
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FIG, 4.—CONTINUED 


in many areas to expose the 6-in. maximum aggregate used in spillway 
paving. Although no evidence of rock from construction operations was founc 
in the basin. There were numerous loose pieces of aggregate in the basir 
that had been rounded and reduced in size to 2-in. or 3-in.-diameter sphere: 
by washing action. Several loose pieces of broken 1-1/2-in. square reinforc- 
ing bars were also found in the basin, Eight or ten baffles were reported as 
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everely damaged as were localized sections of end sill. About half the baffles 
rere reported to have sustained negligible or minor damage, 


RESULTS OF SURVEY OF NOVEMBER 1957 


A detailed underwater survey of the basin was undertaken in November, 
957. A team of two divers using specially fabricated equipment measured 
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FIG, 4, CONTINUED 


evations over the floor of the basin on a grid spaced at approximately ae 
snters and also made detailed measurements of baffle piers and end sil 
-osion. anon 

‘The results of this survey are presented in detail in Fig. 4 and condensed 
| show overall major damage areas on Fig. 5. Fig. 4 indicates the erosion 
Htours with the amount of erosion indicated in inches. The symbols indi- 
€ 
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cated in Fig. 4 are defined as follows: SF indicated a sand finish, with no 
appreciable erosion of concrete; P is pea gravel in relief; AR defines aggre- 
gate in relief, fines eroded from between coarse aggregate, and coarse aggre- 
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FIG. 5.—STILLING BASIN" 
gate not eroded: T.S. denotes temporary sluice, that is presently (1960) filled 


with concrete. Scarcely a single slab within the basin failed to show some 
degree of erosion. Because the aggregate used in paving of the slab was” 
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sraded up to a 6 in. maximum size, any erosion of the slab concrete results 
in development of an extremely rough surface as indicated in Fig. 12, and 
us reported by divers. Major damage areas of sizeable proportions with 
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rosion exceeding 12 in. and more lightly eroded areas with erosion of less 
1an 12 in, are indicated on Fig. 5 by various types of hatching. 
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Fig. 5 shows that only in the “A” row of slabs did erosion exceed 12 in. 
The most severe damage appears concentrated in localizedareas immediately 


in front of baffle piers 


37 to 26, and 3 to 9, between these baffles and between 


the baffle piers and end sill. The areas of greatest slab damage appear to 


Total erosion, March I960 


LEGEND 


Increase in erosion, December 1957 to March 1960 


TAv Average erosion in slab,in feet and hundredths. IAv Average increase in erosion in slab , in inches 
TMi Minimum erosion in slab, in feet and hundredths. IMi Minimum increase in erosion in slab, in inche: 


TMa Maximum erosion in slab, in. 


feet and hundredths. IMa Maximum increase in erosion in slab, in inche 


2) No erosion measurement taken during 1957 su 
visual inspection only was made. 
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FIG. 6.—STILLING BASIN a 
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severe damage. Two small scour holes approximately 5 ft deep were found q 
downstream of baffle piers 22 and 23, and 30 and 31; other areas in which 
deep erosion occurred were generally less than 2 ft in depth. In the “A” — 
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row of slabs, that is the only row with reinforcement, the surface layer of 
steel is extensively exposed in locations at which exceeds 12 in. 

Detailed measurements were made and plotted of damage to each baffle. 
Of the forty-three baffles, it was found that twenty showed erosion of 6 in. or 


Note: The averoge tuilwoter elevation during the stilling bosin survey 
wos 770 feet, U.S.C.&GS, M.S.L. datum. 
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EROSION, 1960 SURVEY 


: 
less and could be considere 


d in satisfactory condition; seven showed erosion 


‘of 6 in. to 12 in., indicating near exposure of reinforcing steel; and sixteen 


showed erosion exceeding 12 in, that expose 


d reinforcing steel. The last two 


‘categories of baffles are indicated by shading on Fig. 5, Damages to baffles 
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TABLE 2,—PAVING EROSION DATA CHIEF JOSEPH DAM STILLING BASIN 


Total erosion March 19604 Increase in erosion Dec 57-Mar 602 


Average, | Minimum, |] Maximum,] Average, Minimum, | Maximum, 
feet feet feet inches inches 
(2) (3) (4) (5) (6) 


inches 
(a) Row of slabs containing baffle piers 


Slab No, 


(7) 


1-A 0,41 

2-A 0,82 8 

3-A 0,91 5-1/2 

4-A 0.30 b 

5-A 0.25 . Af; 2 

6-A 0.33 0.08 2-1/2 

7-A 0.97 0.40 10 

8-A 0,44 0.25 6 

9-A 0.71 0,14 3-1/2 
10-A 1.21 0,28 37 
11-A 0.64 0.23 4-1/2 
12-A 0,97 0,20 3 
13-A 0.69 0.07 3-1/2 
14-A 0.42 0.11 4-1/2 

1-B 0.32 4 

2-B 0,52 3 

3-B 0,48 B 

4-B 0,37 

5-B 0,20 1-1/2 

6-B 0.37 3 

7-B 0.33 4 

8-B 0.39 3 

9-B 0,37 3 
10-B 0.64 6 
11-B 0.52 4-1/2 
12-B 0.56 3-1/2 
13-B 0,28 3-1/2 
14-B 0,36 4 | 

Averages 0,50 : 


1-C 
2-C 
3-C 
4-C 
5-C 
6-C 
75=C 
8-C 
9=C 
10-C 
11-C 
12-C 
13-C 
14-C 
Averages 


soft 
i) 


rPrRoeoecoocooco 


a 
> 
bo bo 


4 Does not include erosion measurement 
considered representative of slab. 
b Visual inspection only, 1957 (except several isolated measurements), 
© Visual inspection made in 1957; no measurements were taken, 


s taken in construction joints which are 
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ceurred principally on the lower front face and sides with generally accom- 
anying damage to adjacent slab areas. However, of the solid shaded baffles, 
umbers 4, 26, 27, and 28 suffered damage principally to the tail of the baffle, 
md number 17 had a large piece knocked out of an upper corner, 


SURVEY OF MARCH, 1960 


To establish the rate of progression of erosion damage and to determine 
thether erosion might be attributable to operating conditions, a second under- 
rater survey was undertaken in March, 1960. The results of this survey 
ulated in terms of total average, minimum and maximum erosion in feet of 
ach slab in 1960 and the corresponding increase in erosion in inches since 
957 are shown in Fig. 6 and in Table 2. In the two flood seasons between 
urveys, peack flows of 377,000 in 1958, and395,000 in 1959, were experienced, 
oth of which are less than 4-yr frequency occurrences. It was not possible 
) make measurements in the “D” row of slabs in 1960 because of the hazard 
f ice falling off the tainter gate assemblies. During the survey, the divers 
rere able to distinguish many light colored fresh pockets throughout the 
asin in which aggregate had been recently eroded from the surface of the 
oncrete. There was also extensive occurrence of fresh impact marks on 
affles and end sill, probably caused by small rounded balls of aggregate 
bund in the basin. 

The average increase in erosion between the 1957 and 1960 survey in each 
f rows “A,” “B” and “C” was found to be approximately 1-1/2 in., although 
ne total erosion in both “A” and “B” rows of slabs appreciably exceeded that 
1 upstream row ‘C.” Correspondingly, the 2-yr increase in damage relative 
> the initial damage measured in 1957 is much greater in the upstream row 
f slabs. Additional surveys in the future willbe needed to determine whether 
tis trend will continue. 

In the 1960 survey, detailed measurements were made of only fourteen 
affle piers including ones representative of light, moderate and heavy dam- 
ge in 1957. In general, baffle piers with surfaces that were undamaged or 
1 sound condition in 1957 were also found to be in sound condition in 1960. 
1ereases in baffle pier erosion was limited to surfaces that were found to 
ave been markedly eroded in 1957. 

Changes in end sill erosion followed the same condition pattern as the 
affle piers with added erosion being generally limited to some areas that 
ad previously been damaged. The overall average increase in erosion of 
amaged end sill areas was only 3/8 in. 


CAUSE OF DAMAGE 


The damage observed in the last 2 yr of spillway service appears minor 
S compared to the pattern of major damage found in the 1957 survey. It 
ppears that most of the damage to the stilling basin has originated from 
onditions during project construction, The pattern of localized, deep and 
xtensive erosion between the end sill and certain baffles and the damage to 
ront faces and sides of these baffles indicates that much of the damage must 
ave resulted from scouring of foreign materials in the basin, particularly, 
ecause there has not been any development of new areas of severe damage 
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nce 1957, An analysis of conditions during project construction has verified 
at on several occasions, rock was introduced into the stilling basin during 
ver diversion, and this situation combined with diversion flow conditions 
eated a potential for major damage. 

Fig. 7 is a composite drawing showing the several stages of cofferdam 
nstruction for river diversion. In the first stage of diversion a cofferdam 
msisting of a series of steel sheet pile circular cells, numbered 1 through 
, were installed to divert the river to the left side of the river channel. 
iring the construction of these cells it became necessary to dump a sub- 
antial groin of 1-ton to 3-ton rock immediately upstream of cells 3 through 
in order to divert high velocity river currents away from the area of con- 
ruction of downstream cells. This rock groin, shown on Fig. 8, was only 
tially removed. Five low monoliths approximately at streambed level were 
mstructed in the first stage cofferdam together with four pairs of 8-ft-wide 

16 ft-high temporary sluices with invert elevations 14 ft above riverbed 
r use in passage of second stage diversion flows. 

First stage construction was completed in January, 1952; the cells were 
moved and second stage diversion, consisting of cells numbered 26 through 
, were installed and operative by April, 1952 so that the 1952 high-water 
is entirely diverted through the low monoliths and temporary sluices of the 
-st stage construction on the right side of the channel. Fig. 9 shows flow 
nditions in March, 1952, with the second stage cofferdam nearly completed 
st prior to high water. Note the white water over the remains of the rock 
oin upstream of the first stage construction and also the outside rock groin 
ing placed at the upper end of the second stage diversion cofferdam. In 
dition to rock wash from these sources, a near failure of cell 29 in the 
cond stage cofferdam during the 1952 high-water required emergency 
mping of rock around the outside of the cell. The high velocity flow around 
2 outside of the cofferdam undoubtedly washed a considerable amount of 
ck into the first stage stilling basin. There was also rock wash into the 
sin from placement of armor toe rock in the wraparound abutments, 
rticularly on the right bank side. With this set of circumstances, the 
ttern of concentrated damage downstream from the low block monoliths 
8, 10, 12, and 14 becomes understandable. 

In the third stage of diversion in 1953, Fig. 10, the cellular cofferdams 
d been removed and flow was diverted through the twenty-four temporary 
tices and the low monoliths that were still at streambed level, During 
od periods a very heavy concentration of flows resulted through the low 
molith diversion notch that, in turn, generated a strong eddy current in 
2 basin, This same pattern of flow was repeated in 1954, except that with 
> raising of the monoliths the diversion notch was now 100 ft above the 
lling basin floor and the corresponding maximum head was 136 ft from 
ol to stilling basin floor. Fig. 11 shows flow through the fourth stage di- 
rsion notch and the location of eddy currents. : 4 
The significance of eddy current action on erosion of the stilling basin 

illustrated by Fig. 13, that shows exposed reinforcing bars and erosion 

the streambed face of the end sill with the end sill concrete eroded to 
pose reinforcing bars bent upstream counter to the cee a of ae 
is apparent that the eddy current was of sufficient strength to pick up roc 
ym the streambed below the dam and carry it into the basin during the 


rd or fourth stage diversion. 


‘rd 
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FIG, 8,—FIRST STAGE COFFERDAM CONSTRUCTION SHOWING ROCK GROIN 
USED AS CURRENT DEFLECTOR 


FIG, 9,-SECOND STAGE COFFERDAM PRIOR TO OCCURRENCE OF HIGH WATER ~ 


AND SHOWING WHITE WATER OVER ROCK GROIN AND BERM CON- 
STRUCTION 


CHIEF JOSEPH DAM 


FIG, 10.—THIRD STAGE CONSTRUCTION DURING HIGH-WATER 
SHOWING EDDY CURRENTS 


FIG, 11.—FOURTH STAGE CONSTRUCTION DURING HIGH 
WATER SHOWING EDDY CURRENTS 
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The initial underwater inspection made in 1957, two years after the 
spillway was placed in full service, showed no rock in the basin other than 
a small amount of dislodged aggregate. This confirmed model tests that 
showed the completed basin in normal operation to be self-cleaning for 
discharges in excess of 300,000 cfs. However, flow conditions throughout the 
diversion period were so different from normal operation that it is probable 


FIG, 12,-UNDERWATER PHOTOGRAPH SHOWING ROUGHENED CONDITION 
OF SLAB AND EXPOSURE OF LARGE AGGREGATE 


that rock washed into the basin was retained as a tumbling agent to scour 
and abrade concrete surfaces. : 

Another damage factor during diversion was that the very heavy concen- 
tration of flows in the fourth stage diversion notch approached maximum 
spillway design conditions, but with an appreciably greater tailwater defi- 
ciency than that for which the basin was designed, as indicated in Table 3. 
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his lack of adequate tailwater depth created a temporary condition for 
otential cavitation damage during diversion. 
; CONCRETE QUALITY 


The possibility of stilling basin damage because of the quality of the con- 
>rete has been considered. Concrete for the stilling basin was designed 


DIRECTION 
OF FLOW 


i FIG, 13. UNDERWATER PHOTOGRAPH SHOWING REINFORCING 
i BARS IN END SILL BENT UPSTREAM. 


for a compressive strength of 3,000 psi. However tabulation of concrete 
properties from construction records (Table 4) indicates, this strength es 
yuirement was substantially exceeded in many areas. Slab and unforme 


surfaces of the top of baffles, end sill and ogee bucket section received 
steel trowel finish. Smooth-surfaced tongue fee rptiee ati 
fo: ? of the baffles and end sill. Concrete wo. 

or formed sections of the te in the north half of the stilling 


for 14 days. A large proportion of concre 
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basin was placed during freezing weather with protection from freezing pro- 
vided by tarpaulins salamanders. No special quality of concrete surface 
was specified. Three possible origins of damage to concrete other than rock 
abrasion are believed possible. 

Erosion of Concrete by Water-Jet Action.—The exposed aggregate condi- 
tion of the stilling basin floor and concave wearing at the base of and on the 


TABLE 3,—COMPARISON OF STILLING BASIN DESIGN AND FOURTH STAGE 
DIVERSION ENERGY DISSIPATION CHARACTERISTICS 


Total flow 


(1) 


Q per foot, cfs| Entrance velocity, 
feet per second 


(3) 
Spillway Design 


400,000 440, 113, 1,0 


Fourth Stage Diversion 


391,000 

(peak) 

350,000 
(exceeded 44 
days) 

300,000 
(exceeded 66 
days) 


23.0 


21.0 
18.0 


4 Dg is theoretical depth of tailwater for jump formation. 
D2! is actual depth of tailwater for jump formation, 


TABLE 4,—PROPERTIES OF CONCRETE USED 


Maximum} Cement 
Content 
in sachs 


Water 


Entrained| Slump | Compressive StrengthD 
content 


Air,? | at plant in psi 
percent 


Lo = 
cation 


(1) 


ee ae ee oy Le 


2 minus 1 1/2 in, 
Approximate overall average: cylinder tests 


Ta Penn > eer ee 


front face of many baffles could possibly be attributed to erosion caused by 
high velocity water-jet action, High velocity flows through the stilling basin 
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| The first erosion due to water-jet action was reported during the initial 
version of high water in 1952 when water flowed up to 15 ft deep over mono- 
ths 7, 9, 11, and 13. Inspection immediately after the flood-flow period re- 
paled light scour and an exposed aggregate condition of concrete on these 
jonoliths below Ele. 775. Extensive tests conducted by the United States 
ureau of Reclamation, Dept. of Interior indicates that concrete will stand 
.e flow of clear water at high velocities without appreciable damage, provided 
vat the flow is smooth and the direction and velocity of the water is not 
pruptly changed. However, water-jet tests at velocities from 100 fps to 
75 fps, showed pitting and surface roughening when the jet angle exceeded 
rg Damage from scour at moderate to high velocities is apparently a random 
| oh depending on the density and quality of the surface attained in 
| nstruction. It is expected that resistance to high velocity scour varied 
.roughout the basin as no special quality of finish was specified. 

- Exposure of Under-Strength Concrete to High Velocity Flows.—Some con- 
rete in the north half of the stilling basin was exposed to high velocity river 
ows ranging up to 40 fps within 60 days after placement. Concrete strength 
: this time may have been below normal because it was placed and cured in 
reezing and subzero weather during November 1951 and December 1951. 
pme concrete was exposed to 45 F temperatures and colder, high velocity 
ows less than one month after placement. 

Concrete Surfaces.—Paving of the slab using concrete with 6-in maximum 
ze aggregate resulted in exposure of large size aggregate with minor erosion 
* the slab surface. The exposure of this aggregate created an immediate 
avitation hazard that with cavitation damage could spread progressively to 
ndamaged adjacent areas. The model tests showed incipient cavitation dam- 
xe for minor deviations from the theoretical curvature on the sides of the 
affles even when using a smooth curve. In the prototype, side curves were 
rmed with narrow tongue and groove lumber, resulting ina marked devi- 
tion from the type of surface tested in the model. It is apparent that any 
‘regularities in the surface resulting from forming, finishing, or patching 
ould greatly increase the possibility of cavitation. In construction, the 
affles were formed only on the sides and front with the top and back slopes 
laced and finished without forms. The use of this open-type forming limited 
1e amount of vibration and compaction of the concrete because of sloughing 
f the back slope, thus creating conditions that caused rock pockets on formed 
reas, Although all such areas were patched, these patches are areas of 
eakness andcould havecreated irregularities or roughness to start cavitation. 


REMEDIAL MEASURES TO PREVENT OR MINIMIZE 
DAMAGE IN FUTURE CONSTRUCTION 


The primary cause of stilling basin damage at Chief Joseph Dam is believed 
) be scour caused by large size rock particles introduced into the basin by 
version flow conditions during construction, It is also possible that damage 
suld have been reduced or minimized by some higher quality control of con- 
rete surfaces. 

A contributing factor to the occurrence of this damage was undoubtedly 
lack of awareness of the damage potential during construction at both design 
ad field levels, Full consideration of the problem during the model testing 
age might have led to changes in design. For example, the tendency of the 


i 
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stepped-end sill to retain rock could have been reduced by sloping the face 
of the sill. It might have been possibleto raise the elevation of the low mono- 
liths 10 ft to 15 ft above streambed to prevent rock from entering the basin 
during diversion, but this would have necessitated higher cofferdams and the 
Single cell-type of cofferdam used was already at maximum height for the 
strength of steel sheet pile that was available. Thus, a double cell-type coffer- 
dam probably would have been required. 

A full consideration of the problem might have indicated a different plan 
of diversion during construction; however, the costs might well have exceeded 
costs of repairing damage to the basin. Greater restrictions placed on the 
operations of the contractor to minimize operations that introduced rock into 
the stilling basin would also have increased the costs, A higher quality con- 
crete surface would have been desirable and steps were taken several years 
ago by the Corps of Engineers to specify superior quality concrete of higher 
strength and better finish and employing small size aggregate for use in 
structures with high velocity flows. 

The steps taken to control damage during construction must necessarily 
be tempered by balancing ultimate costs of repair against the added costs 
of damage-preventative construction measures. 


CONC LUSIONS 


Damage to the Chief Joseph Dam stilling basin reported herein was caused 
primarily by rock entering the basin during the diversion phase of construc- 
tion. A greater awareness of this type of damage potential by both designers 
and field forces is necessary to avoid future recurrences. 
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FORECASTING RIVER RUNOFF BY COASTAL FLOW INDEX 


: By David M. Rockwood,! M. ASCE, and Carlton E, Jencks? 


SYNOPSIS 


An index method for forecasting April through September runoff for the 
‘olumbia River at The Dalles, Oreg., is presented, based on the relationship 
vetween winter runoff in coastal drainage areas and spring snowmelt runoff 
if the Columbia River. Reliability of the forecasts is comparable with those 
ased on precipitation and snow-course data. 


INTRODUCTION 


The Columbia River is one of the principal snowmelt runoff streams in the 
Inited States. Inasmuch as the flow of the Columbia River is primarily from 
nowmelt, an evaluation of existing snowpack water equivalent makes it possi- 
le to forecast the total runoff potential several months in advance. Such 
orecasts enhance the use of multiple-purpose reservoir storage for obtaining 
ptimum benefit of water resources for all water uses. The evaluation of the 
otal runoff potential of the Columbia River is necessary for coordinated 
stem operation of reservoir storage for flood control and other water use 
equirements, as well as for flood potential evaluation prior to completion of 
n adequate system of reservoirs for flood control. The Corps of Engineers, 
Inited States Army, is responsible for flood flow regulation in the Columbia 
fiver Basin. The prime basin forecast of the total runoff is for the gaging 
tation located at The Dalles, Oreg. (D.A. = 237,000 sq miles), at which point 


ee 


Note.—Discussion open until August 1, 1961. To extend the closing date one month, 


written request must be filed with the Executive Secretary, ASCE, This paper is part 
f the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
iety of Civil Engineers, Vol. 87, No. HY 2, March, 1961. 
1 yar, Engr., U. S. Army Engr. Div., North Pacific Div., Portland, Oreg. 

2 Meteorologist, U. S. Army Engr. Dist., Portland Div., Portland, Oreg. 
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the runoff from essentially the entire Columbia River drainage east of the 
Cascade Range is measured. 

Methods of forecasting seasonal runoff volume from snowmelt have been 
based primarily on use of precipitation or snow accumulation measurements, 
The most commonly used procedures employ the concept of indexes, whereby 
point measurements of hydrologic variables are correlated with runoff by 
various statistical techniques. Coefficients so derived from historical data 
provide the basis for forecasting future runoff, by applying them to observed 
index values. Forecasting methods of this type dependupon an adequate period 
of record for statistical treatment, and these index procedures may involve 
several variables affecting runoff. 

One of the objectives of the Corps of Engineers Weather Bureau, United 
States Dept. of Commerce (USWB), Cooperative Snow Investigations, was the 
development of practicable and reliable methods of forecasting seasonal 
streamflow resulting from snowmelt. The emphasis in the analytical work 
that was carried on by the Corps of Engineers’ Snow Investigations Unit in 
pursuit of this objective was in the rational evaluation of the water balance on 
basin areas involving snow, and the determination of the storage of water in 
the snowpack in relation to the other factors involved in the hydrologic cycle. 
The results of these investigations are available.3; 

During the course of the snow investigations, it was realized that index 
procedures are useful in deriving practical forecasting equations. General 
knowledge of the hydrologic cycle and causative meteorologic factors involved 
in the deposition of precipitation, as appliedtoa specified basin, are essential 
for the intelligent selection of suitable parameters for forecasting seasonal 
runoff by index methods. 

Due to the fact that the Columbia River Basin is located in the middle 
latitudes, in the zone of the prevailing westerly air flow, the sources of 
moisture supply are the maritime air masses of Pacific origin, particularly 
during the winter season. Precipitation occurs with the eastward movement 
of large-scale frontal systems from the Pacific. These systems, which are. 
usually occluded by the time they reach the coast, produce varying amounts 
of frontal precipitation. Additional precipitation results from orographic 
lifting of the moist air masses as they pass over the Coast, Cascade, and 
Rocky Mountain topographical barriers. Inspection of upper-air charts” shows 
that, during times of significant winter precipitation over the Columbia River 
Basin, the average air flow at the 500 and 700-mb level is generally from the 
west, but may vary from southwest to northwest. Further inspection reveals” 
that the basin is generally well centered in the westerly stream of air which, 
on the average, covers a width of about 15° to 20° of latitude. : 

The initial orographic uplift of the air stream over the Coast and Cascade 
Ranges results in deposition of a portion of the total precipitable water of the 
air mass before it enters the Columbia River Basin, east of the Cascades, 
The lower elevations and predominant maritime air-mass characteristics of 
the drainage areas in the Coast and Cascade Ranges result in higher compa 


tS 


3 “Summary Report of the Snow Investigations, Snow Hydrology,” U.S. Corps of Engrs, 
North Pacific Div., Portland, Ore., June 30, 1956, 

4 “ Application of Snow Hydrology to the Columbia Basin,” by Oliver Johnson and Pete 
B, Boyer, Proceedings, ASCE, Vol. 85, No. HY 1, January, 1959, 

5 U. S. Weather Bur., daily series synoptic weather maps, northern hemisphere sea 
level and 500 millibar charts, U. S, Dept. of Commerce, Washington, D, C, 
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res which in turn cause a large partof the precipitation to fall as rain, while 
the higher and more continental areas of the Columbia Basin nearly all of 
e precipitation falls as snow. 
The winter runoff from the Coastal streams is considered to be a reliable 
dex of the total winter precipitation falling on the Columbia River Basin, 
}cause it samples the total incoming moisture supply in the air masses dur- 
s times of strong zonal flow. Such an index applied to the Columbia Basin 
- a whole would be a measure of the total moisture supply to the basin and 
erefore may be used as a forecasting parameter. However, it would not 
jaluate the distribution of precipitation on the various tributary streams or 
ib-basins by the application presented herein for the Columbia River Basin, 
With this concept in mind, a project was carried on in 1954 as part of the 
rk of the Snow Investigations Unit, whereby forecasting relationships were 
‘rived for the Columbia River at The Dalles, using graphical correlation 
chniques. The period of record used in those studies was the 22-yr period 
1929 through 1950. A report on the study has been issued.® Subsequent to 
at time, nine additional years of record have become available for study 
951 through 1959). Furthermore, the North Pacific Division office acquired 
e of an electronic digital computer in 1956, which could be utilized for 
ultiple-regression analysis, Accordingly, the analysis presented herein has 
ken advantage of these conditions and is independent of that data previously 
esented,® although much of the basic data is utilized. 


OBJECT AND SCOPE 


A set of correlations of runoff with a coastal winter-flow index and other 
gnificant runoff parameters, by which runoff of the Columbia River at The 
ulles for the period April through September may be determined will be 
esented, Both winter and springtime effects on runoff are evaluated, and 
sults derived by arithmetic statistical methods are presented in graphical 
rm, for the evaluation of effectiveness of the parameters in estimating 
moff, and the range of values normally experienced. The accuracy of fore- 
sts is briefly compared with other methods of forecasting in common use 

the Columbia Basin. 

Considerable use is made of statistical techniques in deriving the relation- 
lips presented herein. It is assumed that the reader is acquainted with basic 
atistical terminology and analytical procedures, and for the sake of brevity, 
) general summary of statistical methods is included. For those who desire 
formation of this nature, textbooks on statistical analysis are available. 


DATA 


Columbia River Near The Dalles, Oreg.—The United States Geological Sur- 
y, Dept. of the Interior (USGS) has published continuous records of stream- 


6 “A Coastal Flow Index Method of Forecasting Seasonal Runoff for Columbia River 
r The Dalles, Oregon,” by David M. Rockwood, Snow Investigations Research Note 
te No, 23, U.S, Army Engr. Div., North Pacific Div., Portland, Oreg., September 30, 


Be eiceinai Methods,” by Herbert Arkin and Raymond R, Colton, (College Outline 
ries), Fourth Ed., Barnes and Noble, Inc., New York, N. Y., 1955. 
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flow for the Columbia River near The Dalles, Oreg., since June, 1878. The 
record for that station is the longest in the Pacific Northwest and is rated 
“excellent” by the Geological Survey. It is the furthest downstream of any 
rated gaging station on the main stemofthe Columbia River, and the drainage 
area is 237,000 sq miles. Annual and April-through-September runoff volumes 
were determined from this record. Appropriate corrections for major reser- 
voir storages, including Kootenay, Hungry Horse, Flathead, Pend Oreille, 
Franklin D. Roosevelt, Coeur d’Alene, Chelan, and Brownlee, were applied to 
observed flows. There are many irrigation diversions and small reservoir 
storages which affect the flows at The Dalles, but their net effect on annual 
and seasonal runoff is small in comparison with the total volume. This is 
particularly true with regard to year-to-year differences. Their effects were 
therefore ignored. 

Coastal Winter-Flow Index.—It has been pointed out that the use of a coastal 
winter-flow index is predicated on evaluation of runoff from areas which are 
at low elevations relative to the Columbia Basin as a whole, and also which 
are well located with respect tothe airstream during average winter air-flow 
conditions. In addition, records of sufficient length for analysis of data are 
required. The 31-yr period 1929 through 1959 was chosen for study, principally 
on the basis of continuity of streamflow records for coastal streams. The 


combination of streams selected for the coastal winter-flow index is as fol- 
lows: 


Stream-Gaging Station Drainage Area, Sq Miles 
Willamette River at Albany, Oreg. 4,840 
Wilson River near Tillamook, Oreg. 159 
Lewis River at Ariel, Wash. | 731 
Chehalis River near Grand Mound, Wash. 895 


Fig. 1 is a map of the Columbia River Basin and the coastal region, showing 
the location of the index streams. 

The Willamette River at Albany was selected because of its long period of 
record and because it samples the flow from the Cascade Mountains as well 
as the coast range. Due to its moderately high average winter temperatures, 
nearly 70% of the annual runoff occurs during the period October through 
March, Records for the station are termed “good” by the USGS. For the 
period through 1950, no correction was incorporated for reservoir storage 
because of the relatively small amount of reservoir capacity above the station 
(about 140,000 ac-ft of storage out of 10,000,000 ac-ft annual runoff). For the 
period subsequent to 1950, appropriate corrections for storagein Fern Ridge, 
Cottage Grove, Dorena, and Lookout Point Reservoirs were applied to observed 
streamflows. 4 

Wilson River near Tillamook samples the air flow north of the Willamette 
River at Albany, and it represents runoff from the west side of the coast 
range. Only 6% of its area is above 3,000 ft and about 80% of its annual f 
occurs from October through March. The records are termed “good,” and 
there are no upstream diversions or regulation. Lewis River at Ariel was 
used for sampling a west slope Cascade drainage area in Washington. I 
average elevation is about 2,800 ft, somewhat higher than the other 


March. The records are termed “good,” and appropriate corrections wer ° 
made for changes in reservoir storage. A lower elevation area in Washing 
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is represented by Chehalis River near Grand Mound. Here only 1% of its 
area lies above 3,000 ft, and 81% of the annual runoff occurs from October 
through March. 
Table 1 shows average monthly runoff for the period October through 
March, during the period 1929 through 1959, for each of the four coastal 
winter-flow index stations. Table 2 shows the total October-through-March 
flow for each station, expressed in percentage of the 22-yr average for the 
eriod 1929 through 1950. The four station average is indicated in column 6, 
Table 2. A minor recession correction, for the purpose of correcting the ob- 
served runoff for water temporarily in storage in the coastal flow index 
streams at the end of March, is shown in column 7. This correction was de- 
termined by recession analysis for Willamette River at Albany, in a manner 
similar to that shown elsewhere.” The coastal winter-flow index (CFI) values 
shown in column 8, Table 2, are those used in the multiple linear regression 
analysis to be described. 
Columbia Winter Runoff Correction Index.—Nearly all of the winter pre- 
cipitation from significant runoff producing areas in the interior of the 
| Seep River Basin falls as snow. Streamflow during the winter is largely 
ase flow recession, augmented by occasional low-level rain and snowmelt 
runoff. Only 28% of the annual runoff occurs in the period October through 
March, Nevertheless, significant variation in winter runoff does occur as the 
result of unusual weather conditions. Inasmuch as the basin moisture input 
indexes were correlated directly with April through September runoff for the 
Columbia River at The Dalles, it was necessary to account for the variability 
of winter runoff. The October through March runoff for the Columbia River 
at The Dalles, corrected for storage in Kootenay, Hungry Horse, Flathead, 
Pend Oreille, Coeur d’Alene, Franklin D. Roosevelt, Chelan, and Brownlee 
Reservoirs is tabulated in Table 3. The winter runoff correction index is the 
product of the October through March runoff and the costal winter-flow index 
(CFI); the values are tabulated in Table 3. 

Winter Temperature Correction Index.--In order to account for differences 
in evapotranspiration losses, a winter temperature index is included in the 
analysis. The index used in this procedure isa 5-station mean for the months 
of October, November, February, and March. Climatological records for 
Revelstoke, British Columbia; Missoula, Mont.; Spokane, Wash.; Yakima, 
Wash.; and Boise, Idaho, were used for this purpose, and values of mean 
monthly temperature are tabulated in Table 4. The temperature correction 
index is the product of the 5-station, 4-month mean temperature, times the 
coastal winter-flow index. The locations of the temperature index stations 
are shown on Fig. 1. i 2 

April Coastal Flow Correction Index.—April 1 of eachyear is a convenient 
point of separation between winter and spring effects of runoff in the Columbia 
River Basin. Springtime effects are considered in this study to be confined 
to the months of April, May, and June. Precipitation during April is normally 
associated with the zonal atmospheric flow pattern characteristic of the winter 
season, Accordingly, the index of moisture supply during April was the coastal- 
flow runoff. In this case it is only for the Chehalis River near Grand Mound, 
Wash. Only the one station was used because the index streams of the Willa- 
mette and Lewis Rivers may have significant amounts of snowmelt runoff 


_ 8 «Hydrology Handbook,’ ASCE, Manuals of Engrg. Practice, No. 28, New York, 1949, 
p. 72, 
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TABLE 1.—MONTHLY RUNOFF IN 1,000 ACRE-FT4 


Water Year October through 
December 


(2) 


(1) 


Average, 
1929-1950 


Accumulated 
Average for 
Period End- 
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TABLE 1,—CONTINUED 


Average, 
1929-1950 
Accumulated 
Average for 
Period Ending 
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TABLE 1,—CONTINUED 


@) (5) 


LEWIS RIVER AT ARIEL, WASH. D. A, = 7315 Miles?,> 


Average, 
1929-1950 


Accumulated 
Average for 
Period Ending 


Average, 
1929-1950 


Accumulated 
Average for 
Period Ending 


4 Data from U. S. Geol. Survey Water Supply Papers, 1929-1956 unpublished records 
by U.S. Geol. Survey, 1957-1959 9 
Stream flow corrected for reservoir storage 
c Estimated by correlation with Silete River at Silete, Oreg. 
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te the month of April, while the Chehalis River is virtually unaffected by 
Dwmelt runoff during April. Table 5 lists the April runoff for Chehalis 
ver near Grand Mound, together with the recession correction applied at 
: beginning and end of the month to account for changes in transitory storage. 
ese corrections were derived from the recession analysis. The “generated 
noff” (Table 5) represents the actual streamflow generated from precipi- 


: TABLE 2,—COASTAL WINTER-FLOW INDEX, OCTOBER THROUGH MARCH, 
i IN PERCENT OF BASE PERIOD AVERAGE, 1929-1950 


Coastal 


Index* 
(CFI), % 
(8) 


Grand 
Mound 
(5) 


eee ts 


es eed | 
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1929-1959 Mean Index - - - 107.6 


4 Column 6 + Column 7 — 


ion during the month of April. The April correction index is equal to the 
nerated runoff times the CFI. _ 
Spring-Precipitation Index.—During the spring and summer months, the 
rmal winter atmospheric flow pattern may be altered, so that the runoff 
m coastal streams may not adequately reflect the precipitation over the 
lumbia River Basin as a whole. The effect of spring-precipitation must 
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be evaluated as one of the factors affecting the total seasonal runoff. There- 
fore, a spring-precipitation index for the Columbia Basin as a whole wa: 
computed, using USWB and Canadian climatological stations having records 
for the period of study and presently reporting on the Columbia River Co- 


TABLE 3,—WINTER RUNOFF OF THE COLUMBIA RIVER NEAR 
THE DALLES, OREG.# 


Water Year October through March Winter Runoff 
Runoffb Correction Index® 


4 In 1,000,000 Acre Ft 


b Corrected for storage in Kootenay, Hungry Horse, Flathead, Pend Oreille, Frank 
lin D, Roosevelt, Coeur d’Alene, Chelan, and Brownlee Reservoirs 


© Product, Coastal Winter-Flow Index (Table 2) times October through March Run 
off (Divided by 100) 


Cee 
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+s 
operative Hydrometeorological Reporting Network. There are sixteen aa 
stations, the locations of which are shown on Fig. 1. Inspection of the geo- 
graphical locations of the stations reveals that they are fairly well distribut 
with respect to basin coverage. The arithmetic sum of the precipitation i 
inches at all stations for the months of May and June for each year was co 
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ted, and the average for each month was determined, These values are the 
ring-precipitation indexes for each of the months and are listed in Table 6. 


TABLE 4,—WINTER TEMPERATURE INDEX® 


October November 
(2) (3) 


ater Year 4-Month |Temp, Cor- 
i Average, rection 
Degrees F Index 


(1) 
1959 
1958 
1957 
1956 
1955 
1954 
1953 
1952 
1951 
1950 
1949 
1948 
1947 
1946 
1945 
1944 
1943 
1942 
1941 
1940 
1939 
1938 
1937 
1936 
1935 
1934 
1933 


(6) (7) 


46.14 
44,21 
38,30 
61,29 
31,33 
59,49 
43,76 


February March 
(4) (5) 


ako pp 
eae : ES St ell al ee ae Bl te Led 
SCONIWAOWWAHENHARROLMWNH 
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4 5-Station Average: Revelstoke, British Columbia; Missoula, Mont.; Spokane, Wash.; 
cima, Wash.; Boise, Idaho, : 

b Mean Monthly Temperature in Degrees Fahrenheit. o 

¢ Product, Coastal Winter-Flow Index (Table 2) times average temperature (divided 


100), ; 


e total index is the sum ofthe May and June precipitation for each year, and 
. May-June correction index is the product of this total times the CFI. 


ANALYSIS AND RESULTS 


It is indeed fortunate that this index is feasible for one of the principal 
ers of the United States, for which forecasting water supply is a prime 
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requisite for optimum use of multiple-purpose reservoir storage. Table ' 
presents comparative area-elevation data for the Columbia River east of the 
Cascade and the coastal winter-flow index drainages. On the basis of the 
average for the four coastal winter-flow index drainages, a smaller percentage 


TABLE 5,—APRIL COASTAL FLOW INDEX CHEHALIS RIVER 
NEAR GRAND MOUND, WASH.# 


April Monthly Recession Generated April Correction 
Water Year Runoff, 1,000 Correction, Runoff, Index> 1,000 
Acre Ft 1,000 Acre Ft 1,000 Acre Ft Acre Ft 
(1) (3) (4) (5) 


Mean, 1929-1959 


4D, A, = 895 Sq Miles 


aoe Product, Coastal Winter-Flow Index (Table 2) times generated runoff (divided t 


SSS eee 

a 
of area lies above 3,000 ft on the coastal index drainage areas than lies abov 
6,000 ft in the Columbia River Basin. In addition to the elevation differences 
there are also differences in air-mass characteristics which affect th 
seasonal distribution of runoff between the basins. Also, when considerir 


b 
2 RIVER RUNOFF 135 


unoff-producing areas, the lower elevation zones in the Columbia Basin 
rincipally comprise the arid Columbia and Snake River plains which produce 
.rtually no runoff. As a result nearly all of the runoff in the Columbia Basin 
‘generated from areas above 4,000 ft. 


ABLE 6.—COLUMBIA RIVER BASIN SPRING-PRECIPITATION INDEX (MAY - JUNE)? 


i Water Year Spring Precipita- May-June 
tion? Index, Correction 
i Inches Indexc 

| 


(1) (4) (5) 

3,960 
3,718 
4,514 
5.405 
2,537 
4,008 
4,033 
4,370 
3,804 
3,443 
2.434 
6.687 
3.859 
3,584 
2,634 
2,030 
3,609 
4,546 
2,580 
1,908 
2,364 
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4 16-Station Average, in inches 

b May and June 

© Product, Coastal Winter-Flow Index (Table 2) times May-June Precipitation Index 
vided by 100) 


On coastal index streams about 90% of the runoff is generated from areas 
low 4,000 ft in elevation. These differences are reflected in the average 
w characteristics, as illustrated by the fact that about 28% of the annual 
noff for Columbia River near The Dalles occurs in the period October 
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through March, while about 73% of the annual runoff is produced on the coastal 
tributaries for the same period. The difference between the two represents, 
principally, the differences in snow-water storage in the respective basins. 
An approximation of annual runoff amounts reveals that on the average, the 
total runoff west of the Cascade Range in Washington and Oregon (about 63,000 
sq miles) is about 150,000,000 acre-ft, as compared with 140,000,000 acre-ft 
of runoff for the Columbia Basin east of the Cascades (237,000 sq miles). 

Derivation of Forecasting Relationships.—The initial work on the coastal 
flow index method for forecasting Columbia River runoff’ was based on a 
graphical statistical solution through use of a coaxial diagram.9 One of the 
principal reasons for the choice of this method was that it could account for 
inter-relationships between the “primary” and “secondary” independent vari- 
ables to be expressed in the forecasting relationship. For example, the effect 
on runoff of spring precipitation is a functionof the winter snow accumulation 
as well as the magnitude of the spring precipitation itself. 


TABLE 7 


Location Percent of Area Above Given 


Elevation 


Columbia River near The Dalles 
Willamette River at Albany, Oreg. 


63 45 
26 13 5 


Lewis River at Ariel, Wash, 42 15 3 1 
Wilson River near Tillamook, Oreg: 6 0 0 0 
Chehalis River near Grand Mound, Wash, 1 0 0 0 


Average for four coastal winterflow 
index stations 


With the advent of the electronic digital computer, the solution of multiple 
linear regression equations became fast and relatively simple for statistical 
analysis. By combining certain independent variables as product functions, ii 
is possible to account for variation in the dependent variables resulting from 
inter-actions between the independent variables. This technique was utilizec 
by the Snow Investigation Unit in the multiple linear regression analysis o 
snowmelt rates, where there was known to be a physical law relating twe 
independent variables as a product. 

For the multiple linear regression analysis presented herein, the primar 
independent variable is the CFI, which indexes the total moisture input to thi 
basin. Inasmuch as precipitation is primarily in the form of snow, the CFI i 
an index of the total snowpack accumulation. The “secondary” independen 
variables of winter runoff, winter temperature, April runoff, and May-J: 
precipitation, are all dependent on the magnitude and areal extent of th 
snowpack in their effectiveness in producing runoff. Accordingly, the secondar 
independent variables are included as product functions with the CFI as” 

‘ 
SE ai ge Ie a tll a eS ee oe 

9 « applied Hydrology,” by R. K. Linsley, Jr., M, A. Kohler, and J, L, H. Paulhu: 

McGraw-Hill Book Co., Inc., New York, 1949, Appendix A, | 
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ealistic and rational combination so as to consider the variations in the 
unoff properly, After having arrived at this concept, it was found that cor- 
elations involving the “secondary” independent variables as the product 
inctions with the CFI actually provided stronger correlations than when they 
jere inserted as simple independent variables in the multiple correlations, 
Multiple correlations based on this principle were computed for several 
ombinations of independent variables for the periodof record. The parameters 
ere finally selected for use in the forecast equation on the basis of indexes 
hich rationally account for the various phases of the water balance for the 
asin, They do not necessarily represent the “best fit” of historic data of all 
arameters investigated, which might be selected without due regard to 
aysical explanation in terms of known hydrologic principles. 
The forecast equation was computed by use of the “stepwise” multiple 
mear regression program10 for the computer. All variables affecting runoff, 
icluding spring precipitation, were inserted in the regression analysis, in 
rder to obtain the most realistic weighting of each of the variables. The 
yuation so derived is equivalent to the forecast equation for conditions 
nown as of July 1 of each year. Forecasts as of April 1, May 1, and June 1 
ere computed by using mean values of each of the variables affecting runoff 
iter the date of the forecast. The basic forecast equation is 


= 0.7428 X, - 0.2836 Xp - 0.678 Xg + 0.0358 X4 + 6.50 X_ + 29.6 .... (1) 


Ihnen Y is the April through September runoff for the Columbia River at The 
alles, in millions of acre-ft; X, is the coastal winter-flow index (CFI) in 
srcent (Table 2); Xo is the winter runoff correction index, in millions of 
sre-ft (Table 3); X3 is the winter temperature correction index in degrees, 

(Table 4); X4 is the April generated Grand Mound runoff correction index, 

thousands of acre-ft (Table 5); and X5 is the May-June precipitation cor- 
ction index, in inches (Table 6). Table 8a lists the values of each of the 
irameters, together with the forecast and observed runoff, for each of the 
ars of study, based on the solution of Eq. 1. Results of this correlation 
“esent a measure of the effectiveness of all parameters used in estimating 
‘asonal runoff, and thereby represent forecasts that would have been made 
; of July 1 of each year. Tables 8b, 8c, and 8d show similar data for fore- 
Sts which would be made from data available as of June 1, May 1, and April 
of each year, respectively. 

The mathematically derived forecast equation can be shown in graphical 
rm to illustrate the relative magnitude of individual variables in accounting 
r runoff variance, the ranges of each of the variables for the period of 
cord, and the effect of the product function in weighting the “secondary” 
dependent variables. Fig. 2(a) is the graphical solution of Eq. 1, for the 
nter components of the forecast equation, for which the CFI is plotted as 
e abscissa and April through September runoff components are plotted as 
dinate values. The parameters of October through March runoff and the 
nter temperature index are shown above and below the central base line. 
‘dinate values above the central base line represent the evaluation of the 


10 “Operating Instruction for Use of the Stepwise Regression Procedure on the IBM 
), File 4260,” by M. A. Efroymson and O, Schriker, Esso Research and Engrg, Co., 
2cess Research Div., Linden, N. J. 
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first two terms in Eq. 1, as a function of CFI and October through March 
runoff. Ordinate values below the base line for specified conditions of CFI 
and winter temperature index provide the solution to the third term in Equals: 
Fig. 2(b) shows similar graphs for the effects of the indexes of April runoff 
for Chehalis River near Grand Mound and May-June precipitation, combined 
as product functions of CFI. 
: The means’ values and extremes of each of the index values for the 3l-yr 
period, 1929 through 1959, are indicated on the graphs. The forecast runoff 
for a specific set of conditions can be obtained by the algebraic sum of the 
brdinates of each of the four components shown in Figs. 2(a) and 2(b), plus 
the regression constant, 29.6 (all in unit of milliolns of acre-ft). In application, 
f course, it is simpler tosolve Eq. 1 arithmetically, but the primary purpose 
bf the graphical solution is to illustrate the effect of the product functions. 
If the secondary independent variables were entered into the regression 
analysis as simple variables, rather than as productfunctions, the parameter 
lines in Figs. 2(a) and 2(b) would be parallel. It is reasonable to assume, for 
these particular indexes, that their effect on runoff is increased with increasing 
amounts of winter snow accumulation, mainly asthe result of the greater snow 
covered area that would generally prevail. Accordingly, the diverging shape 
pf the parameters as shown in Figs. 2(a) and 2(b), resulting from the product 


TABLE 9 


Date of Forecast |Average Error,| Standard Error 


Correlation Coefficient 
106 Acre Feet | of Estimate Coefficient, of Determina- 
106 Acre Feet R tion D (= R2) 


(2) (3) (4) (5) 


‘unctions, is a more reasonable solution than using simple unrelated inde- 
»endent variables. 

Summary of Results.—The comparisons of forecast versus actual runoff 
isted on Table 8, are plotted on Figs. 3. Table 9 lists statistical comparisons 
f April through September runoff forecasts which would have been made, as 
# April 1 through July 1, based on solution of the forecast equation for the 
31-yr period, 1929 through 1959. Early season forecasts are based on assumed 
mean values of runoff index parameters for those variables which would enter 
he equation subsequent to the date of forecast. The mean April through Sep- 
ember runoff for the period is 96,420,000 acre-ft, and the standard deviation 
Mf the runoff is 19,030,000 acre-ft. The standard errors of estimate and the 
‘orrelation coefficients shown above were computed on the basis of following 
1umbers of degrees of freedom: 


April 1 27 (n - 4) 
May 1 26 (n - 5) 
June 1 254 (n - 54) 
July 1 25 (n - 6) 
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FORECAST RUNOFF, APRIL THROUGH SEPTEMBER IN 1,000,000 AC. FT, 


FORECAST RUNOFF, APRIL THROUGH SEPTEMBER IN 1,000,000 AC, FT. 


March, 1961 


0° acre-feet 
R = 0.960 


Z) 
WY 


(a) JULY 1 FORECAST 
0 


OBSERVED RUNOFF, APRIL THROUGH SEPTEMBER, IN 1,000,000 AC. FT, 
COLUMBIA RIVER AT THE DALLES, CORRECTED FOR MAJOR RESERVOIR STORAGE 


Average error = 4.69 x 106 acre-feet 
Standard error = 6.47 x 108 acre-feet 
Correlation coefficient, R = 0.947 


(See Table 1) 


OBSERVED RUNOFF, APRIL THROUGH SEPTEMBER, IN 1.000,000 AC. FT. 
COLUMBIA RIVER AT THE DALLES, CORRECTED FOR MAJOR RESERVOIR STORAGE 


FIG. 3(a),(b). FORECAST VS, ACTUAL RUNOFF 
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Average error = 6,48 x 108 acre-feet 
Standard error = 8.48 x 106 acre-feet 
Correlation coefficient, R = 0.906 


40) 
(See Table 12) 


FORECAST RUNOFF, APRIL THROUGH SEPTEMBER IN 1,000,000 AC, FT. 


20) 


(c) MAY 1 FORECAST 


OBSERVED RUNOFF, APRIL THROUGH SEPTEMBER, IN 1,000,000 AC. FT. 
COLUMBIA RIVER AT THE DALLES, CORRECTED FOR MAJOR RESERVOIR STORAGE 


Average error = 6.72x 106 acre-feet 
Standard error = 9.34 x 10© acre-feet 
Correlation coefficient, R = 0.880 


(See Table 13) 


Syx = 9-34 x 10 ac. ft. 


FORECAST RUNOFF, APRIL THROUGH SEPTEMBER IN 1,000,000 AC, FT, 


OBSERVED RUNOFF, APRIL THROUGH SEPTEMBER, IN 1,000,000 AC, FT. 
COLUMBIA RIVER AT THE DALLES, CORRECTED FOR MAJOR RESERVOIR STORAGE 


FIG. 3(c),(d). FORECAST VS, ACTUAL RUNOFF 


ier =" 25°ie 
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The standard errors for each of the regression coefficients for the variables 
in Eq. 1 are as follows: 


Variable Spb 
Xy 0.1345 
Xo 0.1538 
Xe 0.3100 
X4 0.0132 
Xs 1.150 


Comparison of Results With Other Index Procedures.—The prime forecast 
of seasonal runoff in the Columbia River Basin has been April 1, partly be- 
cause snow survey data as of that date has the most general coverage, but 
also because that date represents the time at which a fairly reliable forecast 
of streamflow may be determined and yet is early enough to take into account 
operational problems that may be incurred onthe basis of the forecast runoff. 
In 1954 the technical staff of the Water Management Subcommittee, Columbia 
Basin Inter-agency Committee, prepared a report for review and comparison 
of forecasting procedures for The Dalles, 11 in which comparative measures 
of accuracy of four procedures in general use at that time for forecasting 
runoff in the Columbia Basin, based on April 1 forecasts are shown, !1a 

The procedures were developed by four United States government agencies, 
including the Corps of Engineers, USWB, Soil Conservation Service, and USGS. 
They were based primarily on snow survey or winter precipitation data, but 
with minor adjustments for other runoff factors. Inasmuch as they were de- 
veloped for different forecast periods, a consistent measure of reliability 
was the standard error of estimate, expressedin percent of the average runoff 
for the forecast period. These values ranged from 11.0% to 13.2% for April 1 
forecasts,114 For the CFI method, the standard error of 9,340,000 acre-ft, 
in terms of the average April through September runoff of 96,420,000 acre-ft, 
is 9.7%. The correlation coefficients of the four procedures ranged from 0.71 
to 0.88, based on April 1 forecasts. The coastal winter flow index method has 
a correlation coefficient of 0.88, as of April 1. 


COMMENTS 


The derivation of this procedure of forecasting snowmelt runoff from the | 
Columbia River at The Dalles is based onthe geographic location of the basin 
with respect to the large-scale air-mass and atmospheric flow conditions © 
during the winter season, as well as the differences of the topographic char- — 
acter, principally with respect to elevation, of runoff-producing areas, between 
the Columbia Basin and the adjacent coastal regions of the Pacific Northwest. 
Because the method is based on an index procedure, some consideration | 
should be given to the various indexes presently used for forecasting runoff. F 
In general, precipitation indexes, either from rain gage or snow course water- | 
equivalent measurements, are based on samples which measure only an in- — 
finitesimal part of the total water supply for the Columbia Basin. ‘ 


11 “Review of Procedures for Forecast Seasonal Runoff of Columbia River near a 


Dallas, Oregon,” Columbia Basin Inter-Agency Committee, Water Mgt, Subcommittee, 
Portland, Oreg., 1954, (a) Table 9, p. 49. 
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In the mountains of the west, where thereis great areal variability of pre- 
Cipitation, the relation between gaged and actual precipitation amounts may 
vary from year to year as the result of meteorologic factors which, for the 
most part, are too complex for a rational analysis by existing techniques and 
intensity of the gage network. As an index, then, precipitation measurements 
cannot account for the random variability caused by inherent deficiencies of 
Measurements by gages or snow courses, as well as the variability caused 
‘by the character of individual storms. The coastal winter-flow index, though 
having limitations, possesses certain advantages over gaged precipitation 
measurements. It samples a much larger proportionof the total water supply, 
as illustrated by the fact that in the Columbia Basin, the ratio is ahout one 
part in 14 in terms of runoff. The accuracy in the measurement of streamflow 
‘is an inherent limitation, but such deficiencies also apply to other index pro- 
cedures for forecasting runoff. The principal limitation is the variability 
from year to year of the average air-mass characteristics, and the average 
atmospheric circulation pattern during times of significant precipitation. 

The winter-temperature index and Columbia River winter runoff account 
for part of the variability caused by the air-mass differences. They occur 
from year to year. As indicated by the comparison of forecast errors, the 
coastal winter-flow index accounts for the variability in runoff as well as or 
‘better than precipitation indexes as measured at presently operated rain gages 
‘or snow courses, for the Columbia River Basin as a whole. 

The effects of spring precipitation are accounted for in this procedure by 
the April generated runoff index for the Chehalis River near Grand Mound, 
and the May-June precipitation index, based on a 16-station mean, Use of 
statistical techniques for weighting individual indexes of runoff should be 
predicated on inclusion of all variables in the multiple regression analysis. 
By correlating runoff directly with just those values of the variables known 
as of April 1, itis possible to obtain an apparently better correlation for April 
1 forecasts than using the correlation for all variables, and then computing 
forecasts for April 1, using mean spring precipitation. Omitting the spring 
precipitation effects in the derivation of April 1 forecasts is unrealistic and 
simply forces a best fit of the historic data with incomplete information. It 
does not necessarily provide the best weighting of individual winter indexes. 

In the particular case of the winter indexes of runoff for the coastal 
winter-flow index method, correlating them directly with April through Sep- 
tember runoff provides a standard error of estimate of 9,190,000 acre-ft, 
compared with 9,340,000 acre-ft for April 1 forecasts, obtained from the 
general forecast equation and using mean values of spring precipitation in- 
dexes. The apparent reduction inthe standard error of estimates is fallacious. 
Similarly, direct correlation of indexes known on May 1 of each year would 
provide a standard error of estimate of 8,220,000 acre-ft, compared with 
8,480,000 acre-ft for the forecast equation. __ 

The use of product functions in multiple linear regression analysis of 
seasonal runoff variance may have application to forecast procedures based 
on snow course or precipitation gage indexes. Product functions should be 
used only where there is known to be a physical process relating two inde- 
pendent variables as a product. 

Application to Other Basins.—Since this procedure deals with largescale 
meteorologic phenomena, its application is intended primarily for large 
basin areas, rather than smaller tributary areas. However, where there are 
{ 
és 
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significant elevation differences between basin areas, it is believed that 
winter low-elevation runoff may be used to evaluate water stored in the snow- 
pack in the middle latitudes on basins smaller than the Columbia River pro- 
vided that the index runoff area is well centered in the air stream dominating 
the basin area. 


CONCLUSIONS 


The coastal winter flow index, in conjunction with other runoff parameters, 
provides a basis for forecasting seasonal (April through September) runoff 
volume for the Columbia River at The Dalles. The historical streamflow and 
climatological records used in the correlations have a sufficient period of 
record (31 yr), considering the number of variables used (six), to insure a 
reasonably sound statistical analysis. The use of the “stepwise” multiple 
linear regression program for the electronic digital computer expedited the 
correlation analysis used in developing the basic forecast equation. The basic 
approach of correlating runoff (the dependent variable) with indexes of all of 
the known significant factors affecting runoff (the independent variables), is 
considered to be hydrologically sound. Apparent improvements in fit of his- 
toric data by direct correlations with runoff of variables known as of the date 
of forecast (that is, April 1) are considered to be fallacious. 

The use of product functions of the primary and secondary independent 
variables improved the forecast relationship, in accounting for runoff variance. 
Product functions are believed to be valid in cases where the effectiveness of 
one independent variable is known to increase with increasing values of a 
second independent variable. Comparisons of forecasting methods for Columbia 
River at The Dalles, developed by government agencies for operational use 
in the Columbia River Basin, have been made by the technical staff of the 
Water Management Subcommittee, Columbia Basin Inter-Agency Committee. 
Use of these comparisons shows that the coastal winter flow index method 
is as accurate with regard to fit of historical data as other procedures based 
primarily upon winter precipitation or snow course data for forecasts made 
as of April 1. 
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EFFECT OF AQUIFER TURBULENCE ON WELL DRAWDOWN?2 


Closure by Joe L. Mogg 


JOE L. MOGG.!—mnr. Philip has called attention to a concept that was new 
» the author. He pointed out that inertial forces, rather than the inception of 
urbulence, initially caused the head loss to vary with an exponent of the ve- 
ocity that was greater than one. 

Mr. Philip’s excellent and numerous references, plus some others were 
eviewed by the author. G. Schneebeli2 describes experiments that also show 
nat the head loss begins to vary with an exponent of the velocity greater than 
ne after a critical Reynolds number is reached, However, dye introduced in 
ne flow, for the purpose of identifying streamlines, showed that the flow was 
aminar and not turbulent in this non-linear region. Turbulence was evident 
t higher Reynolds numbers when the dye became readily dispersed. Thus, 
nese experiments show that there is a broad region of flow between laminar- 
inear flow and turbulent-non-linear flow that is laminar, but not linear. The 
ause of the non-linearity could well be due to the inertial forces. 

The conclusion is, then, that turbulence seldom occurs in the vicinity of a 
umping well because most of the non-linear flow is actually laminar-non- 
inear flow. 


4 November, 1959, by Joe L. Mogg (Proc. Paper 2265). 

1 Field Engr., Edward E, Johnson, Inc., St. Paul, Minn. 

2 “Experiments on the range of validity of Darcy’s law and the appearance of turbu- 
ance in a filtering flow,” by G. Schneebeli, La Houille Blanche, March, April, 1955, 
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NEW APPROACH TO LOCAL FLOOD PROBLEMS® 
Closure by Herbert D. Vogel 


HERBERT D. VOGEL, ® F, ASCE.—The observation by Bernard L. Golding, 
. ASCE, to the effect that communities may make their own observations 
id determinations through local engineering forces and consulting engineers, 
sound in principle. The primary need isto obtain adequate data from what- 
rer sources may be available and make the best possible use of it. 
Legislation has been enacted (1960) by the Congress to permit the United 
ates Army, Corps of Engineers to carry out studies of this kind in water- 
eds and on rivers and tributaries under their jurisdiction. This provides a 
ol for the use of communities throughout the entire country and should en- 
uurage the approach to solution of flood problems that has been indicated. 
However the studies may be made, it is a certainty that results will be of 
nefit to individual communities during years of growth and that savings of 
‘eat magnitude will be made by the Federal Government in the future. 


é Paper 2336). 
January, 1960, by Herbert D. Vogel (Proc. 
5 Chmn., Tenn, Valley Authority, Knoxville, Tenn. 
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THE FOURTH ROOT n-f DIAGRAM? 


Closure by T. Blench 


T. BLENCH,39 F, ASCE.—It is gratifying to find that the comments on an 
tensibly practical paper have touched on so many of the fundamental points 
at occasioned it. However, if the most fundamental issue of all is not to be 
oided in reply, as it could have been had comments concerned only details 
the diagram and its use, the writer must confess to a third “object in mind” 
add to the two stated by Gerald Lacey, discuss its reason, and then answer 
ong lines that include detailed replies automatically. The third object in 
ind was to present practicing engineers with the means to free themselves 
5m the trammels of the pseudo-science of engineering fluid mechanics so 
at they can develop knowledge of pipe flow in terms of their own data and 
2 logic that they acquire in the exercise of their profession. The context 
d the references were intended to aidthis object, whereas the diagram pro- 
Jed the means for co-ordinating data. The reason for the object is that the 
-iter is one of those who deplores the long neglect of disciplined instruction 

basic science and mathematics in North American schools, continued into 
lege engineering curricula. In his opinion the neglect has been responsible 
r the displacement of the applied science of hydraulics by so-called fluid 
schanics necessarily based on evasions, part-truths, dogmatism, and so 
rth. From such a subject no advance in fundamental understanding is possi- 
2; in fact it misdirects research. For example, to get away from pipes, the 
gineering student is taught a Bernouilli equation that allows for friction by 
etending that pressure is the same in all directions about a point, omitting 
2 velocity head of fluctuating motion, often ignoring the kinetic energy cor- 
ction factor, and lumping all errors and omissions into a “loss of head.” 
e may imagine the wasted effort if he were asked to conduct fundamental 
search into the hydraulic jump ona 45° spillway chute at a place at which 
locity had become terminal. Here the shear stress on the floor equals the 
rect stress and one principal stress is zero, the shear stress along the 
yor in the jump is high, the kinetic energy and momentum correction factors 
d the turbulent velocity head at the control sections all merit attention. If 

had no fluid mechanics training at all, but good fundamental dynamics and 
ysics, his chances of success would be reasonable because he would know 
w to state his problem in terms of its relevant factors, translate it into the 
iguage of mathematics as necessary, then decide the extent to which it 
ould be tackled inductively or deductively with or without expert aid. He 
uld also research successfully if equipped to study hydraulics as presented 


@ January 1960, by T. Blench (Proc. Paper 2340). 
9 Prof., Dept. of Civ. Engrg., Univ. of Alberta, Edmonton, Alberta, Canada, 
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in texts40,41,42 that presume an adequate foundation of mathematics and 
physics, or draw attention to assumptions, idealizations and limitations and 
do not avoid an argument by “it can be proved that” or “according to so-and- 
so,” 

To avoid unintended offence to textbook writers the author states that he 
places blame squarely on neglected basic science education, has committed 
errors because of it, has experienced the difficulty of teaching inadequately 
prepared students, and cannot imagine a publisher accepting a non-conforming 
college book. He knows there is a changing attitude to education and that at 
least some colleges, including his own, are making desirable innovations. 

Because of the preceding opinions an attempt will now be made to answer 
specific points in an informal manner that will distinguish definite knowledge 
from various degrees of conjecture. 

The practical problem may be stated as “A Newtonian liquid flows under 
conditions of uniform dissipation of head (energy per unit weight) per unit 
length in a circular pipe of given material. Derive a formula to express the 
mean speed of flow, U, intermsofS, 9, y, and d and the boundary properties.” 

An omniscient being, OB, if asked to solve and explain, would probably 
start from the laws of intra andinter molecular action and combine them with 
a relativistic form of Newton’s second law. These would show that the motion 
would not be steady, in general, but the flow would “roll itself up” into eddies 
whose statistical mechanics could be worked out so that the probability of a 
specified velocity at any point during any prescribed time interval could be 
deduced. Thence the mean behavior at every point could be determined and 
finally the mean velocity of the section in a form . 


U =fn(g 8, dj) v} “roughness”)... s.. see. es se (220 


from which p is missing, for dimensional reasons, and in which OB would 
have to explain how to measure “roughness.” His answer would be the dy- 
namically correct one as it would have been deduced correctly from unas- 
Sailable premises in terms of undisputable laws. An ordinary hydraulics in- 
vestigator, HI, of the classic school would admit he did not possess OB’s 
omniscience so would proceed like his brother physicists who also had to un- 
ravel the mysteries of “fluids” (for example, electricity, heat) whose internal 
mechanism was not observable. That is, he would correlate measureable 
quantities covering the whole available range, and test formulas fitting the 
results to see if they suggested simple relations among dynamical entities - 
force, rate of work, and so forth. His knowledge of general physics would 
cause him to believe that if he were lucky enough to hit a formula identical 
with OB’s true answer it would be simple and would certainly express some 
simple relation among dynamical entities, Such a relation might be reducible 
to minimal type, analogous to “light travels so as to make the time of transit 
a minimum,” “free flow over a weir occurs so that the energy of flow is a 
minimum at a critical section... .” Although he would hope to discover the 
dynamically correct formula by his indirect approach he would expect the 
odds to be against him, particularly if the measureable quantities were not. 
well defined or did not cover an adequate range. As a scientist he would not 
convince himself that he had made a discovery until he had checked it against 


40 «Fluid Dynamics,” by L. Prandtl, Blackie & Son, 1952, pp. 125-130, 


41 «The Mechanics of Turbulent Flow,” by B. A, Bakhmeteff, Princeton Univ. Press, 
Princeton, N, J., 1936, 


42 “Fluid Mechanics and Heat Transfer” by J, M, Kay, Cambridge Uni 
5 .M. p v. Press, Ne 
York, N. Y,, 1957, a a 
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. available information and verified that it contained no physical contradic- 
ms. Then he might allow himself to state that he had found a relation that 
emed to have, or suggest, certain dynamical meanings within a specified 
nge of variable. The relative newcomer on the North American scene, the 
thusiastic fluid mechanics investigator, EFMI, does not have HI’s inhibitions 
t feels free, on occasions, to conjure formulas from the air, consider physi- 
1 absurdities in them to be unimportant so long as the answer looks practi- 
lly good, produce only enough data to justify his assumptions and extrapolate 
th happy abandon. 

HI has developed two main lines of attack. One is to regard the flow formula 
an equation between boundary resistance and driving force justifying the 
btting of non-dimensional shear stress, f, against the Reynold’s Number 
d/v) for various values of anunmeasureable relative roughness. This leads 
the standard friction factor diagram and the discovery of different phases 
turbulent flow. The practical and basic value of this diagram is that it is 
otted from readily measureable variables, except for “roughness.” So, if 
ta can be separated into groups of approximately constant roughnesses, it 
7eS a means of defining roughness in terms of whatever formula HI likes 
try as his current guess at OB’s and, of course, it shows him the range of 
nditions over which “roughness” is relevant to the flow. Then, if he finds 
ut facts justify describing roughness by a “roughness height” e in the rough 
undary range, he can define e according to Eq. 4 if he favors Manning, or by 
. 5 if he favors the fourth root type of formula, or by the logarithmic form- 
4 that EFMI likes, or by any reasonably fitting formula; whatever formula 
likes can be used, after algebraic conversionas with Eqs. 2 and 3, to draw 
es of equal relative roughness on the friction factor diagram in the proper 
ae and to label them with appropriate e values. For the purely practical 
in, any system of e’s is good enough provided the formula on which it is 
sed is a good practical fit, for fixed e, in the range of data of interest to 
n; but, of course, HI will never be satisfied till he thinks he has hit OB’s 
aamically correct formula that would fit the full range of its dynamical 
jlicability. The second line of attack follows logically from the first be- 
ise the phases of turbulent flow shown by the f diagram direct attention 
mediately to the question whether they are associated with some peculiarity 
the velocity distribution. Observations of velocity distribution, although 
ry difficult near the boundary at which conditions are important for ex- 
nations, showed the existence of a thin layer of non-turbulent flow for so- 
led “smooth boundary”® and indicated, as implied by D’Arcy long ago, 
t there was a universal velocity deficiency law within an inner somewhat 
lefinite zone®, 40,41 in the sense that, to suitably chosen nondimensional co- 
jinates all velocity data plotted from their turning points as origin would 
icate one curve there (See Fig. 3). A logarithmic curve is a good fit to the 
er zone and to the wall zone at which one curve is not strictly a fit. 
andt140 has indicated conditions (not those of a circular pipe) under which 
velocity distribution curve could be deduced to follow a logarithmic law 
2r a range. A parabola is a good fit, according to Stanton,’ over_most of a 
e but does not attempt to fit the non-conforming wall zone. A hyperbolic 
ction, according to the writer, fits a trifle better than a parabola, No doubt 
thematicians could find comparably good fits. Now, with a formula fitted 
a velocity curve a different way of looking at a flow formula becomes con- 
jient. The flow formula may be regarded as the consequence of integrating 
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the velocity distribution (not the deficiency one) over the cross section to ob- 
tain Q, and dividing by the cross sectional area to obtain U. If the formula 
used to express approximately, the true velocity curve is a good practical 
fit over most of the pipe then the flow formula deduced from it by integration 
must be practically useful over the same ranges of variables. However, as the 
velocity formula remains speculative, the chances of deducing OB’s dynami- 
cally correct flow formula from it seem remote. EFMI wrongly accepts the 
logarithmic fit to the velocity deficiency distribution over the whole pipe as 
absolute dynamical truth and uses it with various assumptions to deduce a 
logarithmic flow formula. This procedure he wrongly attributes to Prandtl 
despite the latter’s warnings.40 


—~— 
6.56 (r/R)* 
13 [cos h(r/R) -1] 


—2.5 log.(y/R) 
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FIG, 3 


With the preceding background we can now outline attacks along HI lines 
that should automatically answer specific parts of the comments. 4 
First let us consider velocity distribution in a circular pipe under the 
given conditions and use dimensional analysis. Then we may write: . 


ie u=fn(y, R, Umax & v) sec ri ee 


i) = dt Ry Urades Gps V)i Hebe een A eee (24) 


in which y is the distance of a point of interest from the boundary (assuming 
that the boundary can be replaced by a geometric surface), R is pipe radiu 
(with the same proviso), u is the time-mean velocity at y, To is the time- 
mean shear stress at the boundary. Given the mathematical ability we could 
always eliminate e from Eq, 23 by inserting its value from Eq. 24, and the 
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uswer would be of the form: 


ue fn(y, Be Upads Ton.P» v) EE eee Oe (25) 
nich can be reduced to: 
u 
ye aR al V pls Vac Vi Ui) oe Saiuy> dgek' cal (26a) 


DW, Suppose the smooth-boundary case shown by the f diagram, in which e 
in be dropped from Eqs. 23 and 24. Then we could eliminate Umax from 
em and Eq. 26a would reduce to: 


u/V, = fn(y/R, V,. y/v) CSYHOOUH) cots cccte ee terete (26b) 


does not seem possible to believe both this equation and Fig. 102 of reference 
i(a) quoted by Mr. Rajaratnam because that figure is fitted by a single line 
aplying that y/R can be neglected. Also, if that line is really the ultimate 
uth, then we must believe that du/dy has a definite value at the centerline 

a pipe, that is, there is a cusp in the velocity distribution. It is interesting 

note that the line has an equation that gives a mixing length distribution 

the inner 30% of the radius that is utterly different from that found by in- 
srting measured values in the formula that defines mixing length.43 The 
riter emphasizes that he is not seriously disputing practical goodness of fit 
the line over the range of its data, but is disputing its dynamical exactness 
id, thence, the EFMI implication that the dynamically correct flow formula 
is been discovered beyond doubt, unlimited extrapolationis permissible, and 
e subject is forever closed. 

Reverting now to Eq. 26a we note that it is consistent with the general ex- 
srimental finding that, for any boundary, there is a universal velocity de- 
siency curve provided y is not too small. For it does not contradict the 
\ySical possibility of: 


i Sie he = fn(y/R, V,y/v) Fait. ioe aero A cy nk OHS 


.d f-diagram experience shows that viscosity can cease to affect certain 
‘oblems appreciably when Reynold’s number exceeds a certain amount. 
irther, arguing as we did to derive Eq. 26b from Eq, 26a, but eliminating 
instead of Umax, we would get 


u/V,, = fn (y/R, Uae Me) oversee eee @ G5 Oy 8 78 € (28) 
at would not contradict the possibility of nature arranging: 
) (Umax - ¥)/Vy = fn (y/R)- +++ sre ecapene oD) 


r this smooth boundary case. However, the function would be discontinuous 
ee the existence of a laminar zone has been demonstrated experimental- 


This seems to sum up all that is definitely known about velocity distribu- 
ms. For practical reasons, or for analysis, various deductions may be made 
om various assumptions, but the results cannot be dynamically correct gen- 
ally except by fortunate cancellation of errors. For example, the “thickness 


43 «The Mechanics of Turbulent Flow,” by B, A. Bakhmeteff, Princeton Univ. Press, 
inceton, N, J., 1936, p. 73, Fig. 49. 
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of the laminar sublayer .. . given by the equation... (11)” that Mr. Rajarat- 
nam quotes on the authority of reference 15 is one of many quantitative defi- 
nitions possible in terms one of many possible assumptions. It is derived by 
assuming that the laminar zone (reference 15, Fig. 102) is fitted by a straight 
line To = udu/dy, instead of a parabola as it should be for dynamical and 
functional correctness, and that the turbulent zone is fitted by the logarithmic 
formula that we have just demonstrated cannot be dynamically correct. The 
intersection of these two is found by trial and error solution of V+6/y =N 
= 2.5 log N+ 5.5. Some proofs seem to make a kind of physical argument out 
of the constancy of V,6/v, and thence claim special validity. Actually any 
single curve that fits the closely bunched points of Fig. 102 reference 15 will 
have a single intersection with any single curve that is supposed to fit the 
laminar film so will give a fixed V,6 /v for the intersection. It should not be 
forgotten that the laminar film is of indefinite and fluctuating size that cannot 
be measured, so its “thickness” is merely a conventional measure. However, 
the convention may be based on a dynamically correct relation or not, so the 
definition might lead back to something useful. The writer chose to define 6 
in terms of an argument of Prandtl’s adapted from a slightly different case;® 
allowing for the actually fluctuating motion in the laminar film the argument 
can be used almost word for word. The different definitions will not agree nor 
will results deduced from them. 

Another class of not quite perfect deduction concerns flow formulas and 
suggests that there is not necessarily a mathematical objection to channels 
and pipes having the same functional type of flow formula—a point that inter-— 
ested Mr. Lacey. One starts from the deficiency law: 


(Umax ~ ¥)/Ve = fn (F/R) oe eee eee eee eee (30) 


and averages the terms over the cross section (by multiplying by 2 mr x dr 
integrating over the whole cross section, and dividing by m7 R2). No matter 
what one assumes the average value of Uma, is itself, and of uis U. The 
average value of the right side depends on the function and, of course, the 
experimentally determined values of the parameters in the function. If the 
function is taken as A (r/R)2, in which A = 6.56 agrees with Stanton over 80% 
of R, the average is A/2, from the property of a paraboloid of revolution that 
the mean ordinate is half the maximum. If the function is the standard loga- 
rithmic, minus 2.50 loge y/R, that attempts to fit the non-conforming wall 
zone, the average is 3.75 and is practically more accurate than 6.56/2 = 3.28 
that deliberately ignores the wallzone. The author has found that 13 [cosh (r/ R) 
- 1] fits Nikuradse data very closely over 85% of R and gives an average of 
3.435 (See Fig. 3). From the viewpoint of aiming at finding OB’s dynamically 
perfect flow formula it is important to note that any function, fn(r/R), will 
lead to: : 


R 
Cnn We B= yf) ne nee 


in which B is a constant, and even if the whole process is repeated fora 
broad channel the only difference in Eq. 31 will be in the value of B (see Eq. 
34). However, as was shown via Eq. 26a, the right hand side of Eq. 27 is not 
quite a function of r/R alone, so B should be slightly variable, and some doubt 
arises whether we can ever obtain a dynamically correct flow formula from 
B = constant. To preceed from Eq. 31 towards deducing flow formulas we use 
the definition f = 2 g dS/V2 and the dynamical fact that To = 4 g dS, so that 
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y. 31 converts to: 


tee eer OFT Harn Ue OTE Jp es (32a) 
Nectope cm A egies ie: a 5 leat ar ate eek » » » (32b) 
ping to ate Vor eB mG/R).. Se eee (32c) 


q. 32a shows that no further progress towards deducing the form of f is 
ssSible without knowledge of the relation of Umax/V, to the factors on 
nich it depends; this knowledge is not contained in the velocity deficiency 
q. 27. The HI method of obtaining the knowledge would start by obtaining 
ta and studying them, The writer would not like to assert, without extensive 
‘search into literature, that such data as were obtained preceded a decision 
believe the logarithmic type of answer of engineering textbooks. For “rough 
uundary” one can quickly obtain Umax/V,, and e/r data of the Nikuradse’s 
periments on artificially roughened pipe from the diagrams in reference 41. 
the reaction between the two variables is to be proved logarithmic the data 
‘e plotted on single log paper and will give an excellent straight line. If they 
*e to be proved exponential the data are plotted on double log paper and will 
ve an equally excellent straight line of slope 1/8. If the equation of the for- 
er line is inserted in Eq. 32a the standard logarithmic expression for f in 
rms of e/r will result. If the equation of the latter line is inserted it will 
it give the writers fourth-root Eq. 5 but would not be expected to give the 
tter’s exact functional form even were it lucky enough to be the ultimate 
mamically correct one. (An important lesson of these alternative plots is 
at a logarithmic function can agree with a suitable power one over large 
nge of variable). 

The relation of Ujax/Vx to what causes it for smooth boundary does not 
em to have been attempted directly. Instead, the sequence of indirect de- 
ction seems to have been that explained neatly in reference 12, Sections 25 
d 28. Briefly, Prandtl deduced velocity distribution along a plane wall in an 
finite fluid with constant shear stress throughout and far enough from the 
ll for viscosity to be ineffective (reference 40, pp. 125-129). Having deduced 
as logarithmic he proceeded, by an ingenious but not quite convincing argu- 
ent, to extend it to fit the zone where viscosity mattered and came out with 
€ expression: 
k u/V, = log V, y/v +C; E Sai Ries tol, ea (83) 


» then pointed out that it would fit Nikuradse’s experiments on smooth pipes 
S an approximation” and that “for flow inpipes” there was actually a “devi- 
ion from theory.” The writer does not know of his having recommended or 
pported, as dynamically exact, integration of this flow formula to give a 
2an value of U, and thence an expression for f. However, it is this integra- 
yn with-the constants fitting Nikuradse’s smooth pipe f data, that gives the 
arithmic smooth boundary line on the Moody diagram, and its formula is 
ind associated in textbooks and articles with the names of Prandtl and von 
rman in a way that seems to have caused engineering readers to believe 
ut these hydraulicians deduced and recommended it. So, when Mr. Rajarat- 
m says “it is widely known that for values of the Reynolds number higher 
in 105 the flow resistance is better approximated by the Prandtl- Karman 
uation. . .” the writer feels that he has been misled by ESMI and the state- 
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ment should read “itis widely mis-stated by textbooks that . . . by the wrongly- 
named Prandtl-Karman equation.” All that the Nikuradse smooth-pipef values 
show is that even a so-called “smooth” pipe starts to become rough at high 
enough Reynold’s numbers, If his “smooth” pipes had been replaced by even 
“smoother” ones his smooth line should have been expected to start turning 
up at higher Reynold’s numbers still; if they had been a little less “smooth” 
then presumably he would have used them to obtain the constants in the loga- 
rithmic formula although they now pass off as extensions of the Colebrook- 
White ones and not “smooth.” Problem 5 was framed deliberately for the in- 
terested reader so that he could decide for himself whether there was some- 
thing dubious about the logarithmic smooth boundary formula at Reynolds 
Numbers higher than those used in the lab. Problem 4 was to convince him 
that, if the Moody Diagram is wrong, it may be alarmingly wrong for very 
large conduits. 

A definite quantitative answer to the queries about canal distributions is 
prevented by inadequate data. As the writer’s interest in flow formulas arose 
from the possible dynamical implications of Mr. Lacey’s regime formulas, 
he has never been impelled to arrange observations of his own to determine 
the constant in pipe velocity deficiency formulas that would suit channels. His 
inquiries from laboratory workers on channels and his long field experience 
with canals support, very strongly, the view that the functional form of the 
velocity distribution down the centerline of a broad canal is that along the 
radius of a pipe. Narrowness probably does not change the functional form 
except by introducing an extra non-dimensional ratio such as b/d and altering 
the value of a parameter. Exact determination is difficult in laboratory flumes- 
because they are narrow and the absolute depth is small. In the field current 
meters are useless near the free surface and the bed, and many beds are 
mobile. Moreover, even in lined canals the water is often turbid and laboratory 
observations show that turbidity can alter, radically, the coefficient 2.5 in the 
logarithmic deficiency formula—one result given privately by a well-known 
authority showed a change to 4,8 with a suspended load concentration of 15.8 
grams per litre, but a very good agreement with 2.5 over 85% of the depth for 
clean water. Finally wind is a major disturbing element in the field. A mathe= 
matical approach to the problem seems to show a contradiction between the 
two common practices of assuming that the velocity deficiency expressions 
are the same, complete with constants, over the radius of a pipe and the depth 
of a broad canal and that, inuse of the Moody diagram or the Manning formula, 
a channel can be replaced by a circle of the same hydraulic radius. If the 
velocity expressions are the same then the right hand of Eq. 31 becomes; 


BoE fh nS) ar eee cece ees (0M) 


that will not give the same value of B for the canal and the pipe. (Here ris 
depth from surface and R is overall depth). But if this is so then Eq. 32a 
shows that the velocity distributions are different for the two cases, in whi h 
it is supposed that Vx and f are the same. This raises the question whether 
the parameter in the velocity formula might differ between pipe and channel 
in such a way as to make B+ fn (r/R) the same for both. None of the functions 
of this article seems adequate for the purpose, but a mathematician mig 
discover one. However it should not be forgotten that B is not quite constant 
in fact. Physically the writer has not been troubled by the difference between 
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irface and center flow of canal and pipe. He has visualized mean surface 
ow in terms of the mathematical identity of the two cases (1) a vortex ap- 
roaches a rigid plane wall and (2) equal and opposite vortices approach each 
her. But another way he has visualized the mean motion at the surface of a 
nannel as if the mirror image were superimposed to turn it into a pipe. 
erhaps the whole preceding discussion will explain why the writer regards 
ie integration of speculative velocity distributions to obtain a flow formula 
5 analogous to using the presumed behaviour of the constituents of the atom 
» Obtain the simple laws of electromagnetism. The cautions procedure seems 
» be to use Eq. 32c to fix the constants of a practical velocity distribution in 
‘rms of f. 

The commentators’ interest in the justification for using channels to aid 
1 testing for a universal flow formula obviously calls for a little explanation 
1 amplification of the references. The writer’s interest in dynamics and 
nysics, his years of practical canal engineering, andhis study of the success 
ith which Mr, Lacey developed his regime theory along the logical lines of 
‘assical physics44 have all helped him to think of the flow formula in terms 
‘ either the equating of resistance to driving force or the selection, by na- 
ire, Of one particular flow for any one problem from among the infinite num- 
=r of flows that one can imagine as kinematically possible. It is difficult, 
ith an interest in physics, to believe that a formula giving the relationship 
stween dynamical entities in simple observed circumstances will be com- 
lex, or that its functional form will be altered seriously without a radical 
qange in the dynamical set-up. Complications may be added by geometric 
Inditions, or by extraneous phenomena being measured with the ones of 
iterest, but behind them is the simple relation. Fundamentally, in pipes of 
ither kind of boundary, or in channels of mobile boundary, all carrying New- 
nian fluids, the boundary “puts on the brake” and the flow responds by 
dling-up into an eddying form that does not seem to be altered essentially 
y alteration in the nature of the brake. Why, then, should the functional form 
‘ the flow equation be different for the three cases or, for that matter, for 
1annels with rigid boundaries? The answer seems to be that there should be 
) difference, but that the term that expresses the “braking action” will have 
fferent detailed expressions for the different cases. At one time the author 
‘ied to prove the point from somewhat advanced dynamical principles but 
eked the necessary skill. However, the normal HI procedure that he likes 
. think of as “plotting formulas” instead of points, had to be followed. That 
, a search had tobe made for a thread running through the various competing 
rmulas for smooth and rough rigid pipes and for mobile boundary canals 
he last made available in physically suggestive form by Mr. Lacey’s monu- 
ental work.44 The thread discovered is the equation given in the paper un- 
2r the heading “Universal flow formula,” with its values of x. The “theory 
' flow,” such as it is, amounts to believing, till something more definite 
rns up, that the odds against this functional form having no physical signifi- 
aynce are very high. Accepting this belief the form merits a try in practical 
3e and, far more important physically, ought to initiate a search for a dy- 
mical theory that would prove or disprove it. It is not believed that the ex- 
sessions for x are beyond improvement; in fact, for the regime case, it 
2ems very likely that small b/d ratios will require an adjustment, Obviously 


eS ee ee 
be : : d- 

44 Flow in Alluvial Channels with Sandy Mobile Beds,” by Gerald Lacey, Procee 
vs. Instit, Civ. Engrs., Vol. 9, February, 1958 pp. 145-164, plus discussion, Vol. 11, 
), 219-251, October, 1958, 
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many more data are required from the field to answer Problem 5 finally and 
thereby test the theory decisively over the practical range for smooth bound- 
ary. The data will have to be numerous, and from different observers, be- 
cause accuracy of measurement will not be too good and the danger of acci- 
dental bias is high even in laboratory work. 

Specific replies to remaining points are: 

(1) The Prandtl proof is not valid for a canal. The outstanding reason is 
that it rests on the assumption of constant shear stress, but canal shear is 
proportional to depth from surface, Flow of wind along the earth’s surface 
probably provides nearly ideal conditions not too close to the ground and not 
too far up. 

(2) Referring to the Webster and Metcalf data, the variation of n does not 
seem sufficient, in terms of reasonable errors of observation and the uncer- 
’ tainty of where to measure d in a corrugated pipe, to pass judgement on the 
matter. Note that, for a known discharge, n varies as the 2.67 power of di- 
ameter. 

(3) The reference to Morris’ work is appreciated. The author should have 
quoted it when he made the qualification, inthe “Introduction,” that the Moody 
diagram applied for “. .. rigid boundary roughnesses that can be expressed 
in terms of a single roughness height.” The oversight is regretted. 

The recent references by Mr. Van Sickle merit study in original. The 
writer believes they emphasize his own view that deviation between the loga- 
rithmic and the exponential formulas is too small over most practical ranges 
to permit a decision to be made between them. For the rough boundary loga- 
rithmic formula the best present test, in terms of real data, seems to be as 
in “Inconsistence of Manning Formula and Moody Diagram,” and for the 
smooth formula the best test is to study the few very large conduits, with 
Reynold’s Number greater than about 5 x 108, in which conditions are a bit 
off smooth and there are enough data to show that the transition curve is con-. 
sistent with Blasius (inthat it could bend into the Blasius line) and inconsistent 
with the logarithmic curve in that it heads right across it. Thus Fig. 1 of 
reference 29 is for Reynold’s Number well below 106 and indicates only that. 
the aluminum pipe appears to be just starting to go into transition on its way 
towards roughness; some of the points are below the log line a little and the 
Figure has the logarithmic line incorrectly drawn in its upper portion. In 
reference 31, Fig. 3 shows four sets of observations with marked inconsis- 
tencies that puzzled the commentator, Mr. Campbell, and caused the author, 
Mr. Burke, to admit that he believed conditions were not of uniform flow. Mr. 
Campbell noted that the one set of the four that gave the approximate trend 
of the logarithmic smooth boundary line fell below it. A very interesting fea- 
ture of the reference is the plot of observed velocity distribution in its Fig. 7. 
This plot deviates from the so-called Prandtl one (that was not for pipes) 


Finally, references 35 and 36 are very informative but do not seem to add 
significantly, for purposes of the present discussion, to the discussion of the 
paper of Hickox, Peturka, and Elder. They include Raam’s method of est ~ 


and correlates it with f. This procedure Ought to make the f-n diagram, or 
any other, useable outside the range where roughness height becomes large 
compared with diameter and is on a par, for practical utility, with using ex- 
perience, rather than formal curves, in the transition zone, 
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JNIFORM WATER CONVEYANCE CHANNELS IN ALLUVIAL MATERIAL2 


Discussion by Marcel Bitoun 


MARCEL BITOUN,!9 M. ASCE.—The writer was interested in this paper 
om the standpoint of the engineer who is confronted with the design of large 
‘rigation canals. The scope of the material covered is the widest yet pub- 
shed, and their review of the data collected by others and by themselves is 
most valuable contribution to the study of the behavior of man-made alluvial 
1annels. 

Careful reading of the paper and examination of the basic data reveal that 
ost of the canals carrying flows smaller than 1000 cfs follow the various 
orrelations studied relatively well. However, as was pointed out, the larger 


TABLE 8,—COMPUTED MANNING COEFFICIENTS 


Canal Number Discharge, cfs 
(1) (2) 


nals, conveying flows in excess of 1000 cfs, seem to depart markedly from 
e general trends indicated by some of the correlations. In particular, Fig. 8 
ows a limitation of the mean velocity to approximately 3 fps in the Punjab 
1d Sind canals when R“S becomes larger than 0.005. 

This can be expressed differently: if the Manning roughness factors n are 
ymputed by application of the Manning formula to the Punjab canals, one 
‘tains values substantially higher for the larger canals than for the smaller 
nals. For comparison purposes, an average value of n for the smaller 
nals is 0.0215. See Table 8. 

The Blench-King formula, Eqs. 15 and 29, is also inaccurate. Table 9 shows 
e results of the application of the Blench-King slope formula and the Lacy 
rmula to the large canals. 


4 May 1960, by D. B. Simons and M, L. Albertson (Proc. Paper 2484), 
10 Warza Engrg. Co., Chicago, Ill. 
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It is apparent that the application of these formulas lead to slopes sub- 
stantially flatter than the actual slopes. 

Although the authors indicate that a reasonable estimate of the design 
slope can be obtained by using the correlation of V and R2S or by using the 
Blench-King regime slope formula, this accuracy does not appear to be suf- 
ficient for design purposes, as the determination of the slope of irrigation 
canals is one of the most critical items. 

A reliable formula for flow on movable alluvial materials should give con- 
sideration to the grain size of the bed material and to the bed configuration. 
Such a formula was proposed by H. K. Liu, M. ASCE and S. Y. Hwang, A.M. 
ASCE (69). The validity of its form was verified on the Punjab canals, and 
computed slopes for the large canals are shown in Table 10. 

The coefficients Cp listed in Table 10 are slightly different from the co- 
efficients Cz recommended by Liu and Hwang. They were derived by the 
writer from an analysis of the Punjab data. The necessity of making minor 


TABLE 9,—COMPARISON OF REGIME SLOPE FORMULAE 


Measured Slope 
(per thousand) 
(2) 


Computed Slope 


Blench-King 
(3) 


adjustments to Cy for application of the formula to different areas is attributec 
by the writer to the fact that the mean grain size is not sufficient to com. 
pletely define a bed material. 9 
The accuracy of the formula 


AE re so Urabe ar re hR a bol ect BC (70 


is definitely better than that of the regime slope formulas recommended D1 
the authors, It is within acceptable limits for design purposes and it is sug- 
gested that its use be considered. q 

Another anomaly is the fact that the Punjab canals, reputed to be stable, 
are characterized by widely different tractive forces on the bed and largely 
different width-depths ratios. According to Lindley, Lacey, and Blench, fo: 


canals constructed in similar soils, the relationship between the width-dep ' 
ratio and the velocity is | 


W 
Fy Oe Wi igae 2 ghee casemate (71 
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Analysis of the Punjab data, whose soils are described as relatively uni- 
rm, does not confirm this relationship. 

The correlation of W/D and V is shown in Fig. 17(a), The width used in 
is correlation is as defined by Blench: W = A/D. 

It is apparent from Fig. 17(a) that the various canals are characterized by 
dely different values of K. For example, in comparing large canals No. 2, 

7, 8, and 9, the bed material mean diameter varies only from 0.41 mm to 
43 mm, and V varies from 2.54 to 3.25 fps, whereas K varies from 4,75 to 
75. For the ensemble of the Punjab canals, K is comprised within an even 
der range: 0.70 to 7.75. 

However a comparison of groups of canals operating under similar con- 
tions of velocity reveals the following: 


1. For the same mean diameter of bed material, the tractive force on the 
d increases up to a certain value as the width-depth ratio decreases. 

2. The maximum value of the tractive force exhibited by the canals and 
e mean diameter of the bed material both increase as the width-depth ratio 
creases, 


TABLE 10,—SLOPES COMPUTED BY THE LIU-HWANG FORMULA® 


Computed] Measured |Blench- 
anal No,| Material Slope (per}|Slope (per King» 
thous and)|thousand) 


av = Cy) RX SY 
b Slopes computed by the Blench-King formula as shown on Table 9, 
© Slopes computed by the Lacey formula as shown on Table 9, 


) illustrate this, Table 11 shows the values of T, d509 and W/D for the group 
canals operating at velocities of approximately 1.4 fps. 

It is suggested that some of the canals could have been constructed with 
lower value of W/D without danger of bed scouring. For example, canal No. 
a in which the tractive force is 0.026 psf could have been operated with a 
uctive force as high as that of canal No. 14 (0.044 psf) without being eroded. 
Thus, instead of searching for relationships common to all the canals under 
yestigation, envelope curves corresponding to the limit of stability should 
defined. Such curves have been drawn on Fig. 17(b). They represent the 
nimum width-depth ratios that characterize the stable Punjab canals and 
duld correspond to tractive forces on the bed close to the limit of stability 


or. 
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mpatible with a mean grain size of bed material. These curves are pre- 
nted here only in order to illustrate the concept of limiting conditions. 
irther investigations are necessary to determine accurate numerical values. 
The writer grouped the forty-two Punjab canals according to their eleven 
sdian grain sizes, then selected the highest tractive force value for each 
Oup and tabulated these values opposite their respective grain sizes as 
lows: dso 


max T, 
(mm) in psf 
0.43 0.110 
0.42 0.111 
0.41 0.106 
0.34 0.084 
0.30 0.088 
0.29 0.052 
0.28 0.047 
0.26 0.054 
0.19 0.048 
0.17 0.049 
0.15 0.050 


TABLE 112 


Canal Number 


(1) 


. 


DMDNUOWHOMNEA Wo 
PrmNwWwwwnnnHPworHHEE 


4 Group 1 corresponds to d50 
Group 2 corresponds to d50 
Group 3 corresponds to d59 


GN eS ee ee 
ee 


These values of T canbe compared to the shel nee eee by E. W. Lane 
als. Fair agreement is shown by Fig. 16. 

it ively awatire that, for similar mean velocities, the bed ae 
ial in the stable Punjab canals becomes coarser as the width-depth ratio 
‘reases, This is put in evidence by the family of envelope curves ae on 
‘, 17(b). Considerations of fluvial dynamics indicate that when eR tee 
ne is established, the material that forms the bed of an alluvial c mii me 
wracterized by a coarser mean diameter than the original soil in whic e 
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channel was excavated. During the transition period that precedes establish- 
ment of regime conditions, the material in place is subjected to a process of 
sorting or segregation. The fines contained in the bed material are put in 
suspension by turbulence and carried away as wash load. If the material ir 
place contains a sufficient proportion of grains witha size such that they can- 
not be removed, an armor of a certainthickness develops on the bed and pre- 
vents further segregation. The median size of the armor increases with the 


Ib / sq. ft 


Lane -Recommended values for canals: = 
With high content of fine sediment 


in the water. 


TRACTIVE FORCE , 


in the water. 
With clear water. 
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FIG, 18.—LIMITING TRACTIVE FORCES 


forces acting on the bed. However, although this bed is not being eroded 
stable alluvial canal, it is moving as bed load. Stability requires that the 
load movement of the canal bed be replaced by a supply of solids into 
canal at its headworks, equal in size and in quantity. 
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In the case of the Punjab canals, it could be that if some of the wider canals 
ith relation to depth) had been constructed somewhat narrower, their bed 
aterial would have been segregated to a larger size without causing further 
stability. They did not silt probably because the solid charge was not ex- 
issive. 

It would thus be possible to predict the extent to which segregation will 
cur in a new canal onthe basis of the tractive forces imposed, the gradation 
the material in place, and the solid charge permitted to enter the canal at 
2 headworks (or imposed by the parent river when sand control devices are 
it present). The design can be conducted soas to allow segregation affecting 
2 bed material up to a large grain size, if such grains are available in a 
fficient amount in the material originally in place. This would lead to lower 
dth-depth ratios. However, the width is limited in practice by the fact that 
osion of the banks could then occur. Stability of the banks could be deter- 
ined by relating the tractive force on the banks to the friction angle of the 
aterial in place and the slope of the cut (or fill). 

However, stability cannot be insured unless the bed load movement is de- 
rmined for the chosen size of canal and segregated bed material size, and 
less the supply of sediments into the canal at the headworks is maintained 
ual in size and in quantity to the material migrating downstream as bed load. 

When discharge and sediment concentration in the parent river is variable 
roughout the year, stability of the canal must take this into account. The 
iriation of sediment concentrations should be correlated to the variation of 
scharge in the canal imposed by irrigation requirements, An integration 

the aggradation and degradation should provide a yearly cycle resulting in 
erall stability. 

In that respect, it would have been interesting to find among the data pre- 
nted detailed information on seasonal variations of flow in the canals in- 
stigated and on the sediment allowed to enter the canals. 

An approach to the regime theory should not be based on a statistical 
alysis producing averages. The recognition of limiting conditions leads to 
sults equivalent to those obtained from an approach based on the tractive 
rce theory. These approaches, however, should be supplemented by due 
nsideration given to solid transportation as one of the determining factors 

stability. Addition of a bed load formula to the basic equations of stability 

necessary and should enable a complete determination of the regime of 
uilibrium. Yearly average stability implies an integration of allowable ag- 
adation and degradation with respect to fluctuating sediment supply and 
iter requirements. 


ADDITIONAL BIBLIOGRAPHY 


“Discharge Formula for Straight Alluvial Channels,” by H. K, Liu and 
S. Y. Hwang, Proceedings, ASCE, November, 1959. 


MONIR M. KANSOH.!1_This study of existing stable canals in the United 
ites is actually a continuation of the regime studies that were initiated in 


11 Prof, of Hydr, and Water Power Engrg., Alexandria Univ., U.A.R. 
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the Punjab by Kennedy and that yielded the empirical equation: 
FALE 9 LA eda pe A em . (72a) 


using foot units, or 
0.64 
Vip 0-546 Dn one reie einem iain a) 94h sind (72b) 
using metric units. 

The diameter of the suspended material in the River Nile in Egypt © 
much less than 0.20 mm with a concentration varying from about 3500 ppm 
. down to about 200 ppm with an average of about 1600 ppm, nearly similar to 
the conditions in the Sind in which the “mean size of bed material is within 
the limits 0.0346 mm to 0.1642 mm” and “the suspended sediment load ranges 
from 3590 ppm to 156 ppm.” This isto be compared with the conditions in the 
Punjab in which “the average diameter of the bed material is approximately 
0.43 mm” and the mean silt intensity is, like that in the canals investigated 
by Simons and Bender, about 328 ppm, According to the authors, the tractive 
force theory may not be applied in Egypt where water is not clear and where 
the bed and sides of the canals are, in most cases, composed of fine cohesive 
materials. Design of stable canals should, then, depend on regime formulas. 

In Egyptian canals, because the silt is finer than that in the Punjab, the 


critical velocity is sometimes taken as a of the value given by Kennedy: 


3 
0.64 
V_ = 0.56 D REE, a Pr Ware ee es (67/5 
using foot units, or 0 (13a§ 
Vi = 0.96.D tots «+1520 .0 aav 


using metric units. 
Observations on thirteen canals in different parts of Egypt by K. Ghaleb 
produced the empirical equation 


NgoriDBA Diet anpecove- rime (v4) 


using metric units. 


These three values are rather excessive and it was recommended, espe- 
cially for canals in Lower Egypt, that the formula 


4 


Vo separ Dh. cal cherie’ ti brcuete) cakes ae a el (75) 


be used in which m = 2 whereas 5 varies between 0.26 (for fine stilt) to 0. 30 


(for coarse silt), using nietinié units, 3 


The form of the canal cross section is supposed to be predetermined. If 
it is, as is usually the case, a trapezoidal cross section with a bed width Ba ‘ 
a depth D and with sides sloping at an angle @ with the horizontal, then 

A=(B+Dcot 6)D=Q/V=Q/(pD")......... . (76a) 
= (A/D) - D cot 6 = Q/(p D™*+) - Deot 6 ....... (76b) 


P= Q/(p D™*+!) + D(a cosec 6 - cot @) oa: cue ae 


Sei gees + D (2 cosec 6 - cot a ry: 


To fix the dimensions of a stable canal in order to convey a discharge Q, 
the critical velocity as given by the empirical regime formula is substitute 


and 


for) 
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io the Chezy equation with the aid of one of the formulas given by Bazin, 
unning, Kutter, and others. The Manning formula is often recommended 
cause, in addition to its simplicity, it gives satisfactory results. 

Then 
oS, p™ - ui8) R2/3 


ot units) or 


etric units). 

‘The depth D can be computed after putting R in terms of D as previously 
own, 

The bed width can be directly computed because 


Ba asa 
(p D ) 
‘The computation is much simplified and even avoided by using diagrams, 

mograms, or tables. 

‘The relation between B and D, or between W and D might as well be ob- 
ned from empirical regime equations like Eq. 4, or Eqs. 13 and 14. 
After a careful examination of a large number of recognized canals in 
ypt, Molseworth and Yennidunia recommended the following relations be- 
sen B and D: 


1/2 


D=0.10(> +4)B Shad dep oat Se tage 


DB 5.0.00154 Bi. 8)" Bi. .1,ca lene ccbsorinns « 2AT9D) 


D above 1.62 m and below 1.62 m respectively, S being the water slope in 
canal in centimeters per kilometer. 
They also recommended the ratios 


D=Cy Boece ec ecee ee cce see» (80a) 


=n B is less than 2.00 m and 
p-c, BY? ta 5 EAM GB yok Seep (1) 8) 


2n B is greater than 2.00 m, in which Cy = 1.00 and Cg = 1.75 for deep 
ins while Cj = 0.90 and Cg = 1.45 for shallow drains. 

These relations are represented in Fig. 19, that gives the desirable pro- 
tions of stable trapezoidal canals and drains with 1:1 side slopes. The 


gram consists of two groups of curves: 
1. The dotted curves that correspond to the value 


K-ar”/9-9 (2 s/?) See eae AOD 
ng metric units. p 
2. The ratio-curves corresponding to canals of different slopes, deep 
ins, and shallow drains. 
The point of intersection of the appropriate K-curve and the proper ratio- 
ve determines the desirable depth and bed width of the canal. 
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Interpolation is sometimes necessary when the computed value of Q/| gl/2 


ies between two K-curves, or when the slope lies between two of the slopes 
mdicated by the ratio-curves, 

The same relations of B and D, for canals, are the basis of a standard 
liagram used in the Egyptian Ministry of Public Works for the design and 
maintenance of stable trapezoidal canals with 1-on-1 side slopes when the 
solids in suspension do not exceed 2000 ppm. This diagram, Fig. 20, is suit- 
ible for one class of canals whose Manning roughness coefficient n = 0.025, 
while Fig. 19, suits any value of n. 

To use the standard diagram, Fig: 20, a vertical straight line indicating 
he discharge is extended to intersect the right hand inclined straight line 
mdicating the chosen or selected slope. From the point of intersection, a 
iorizontal straight line is extended to the left. This horizontal straight line 
ind the left hand inclined straight line corresponding to the same selected 
slope intersect at a point that directly indicates the desirable depth and bed 
vidth of the canal. 

The procedure may sometimes be altered by fixing the depth and trying to 
ind slope and bed width corresponding to the given or required discharge. It 
nay be the bed width that is first known while the depth and the slope are to 
»e determined. 

Hesitation usually accompanies the selection of the suitable slope indi- 
‘ated by one pair out of the two groups of inclined straight lines. To reduce 
his hesitation, H. El Defrawy suggested a particular selected slope for each 
‘ange of discharge, depth and bed width. These slopes are shown in thick lines 
n Fig. 20. 

A. M. El Banna proposed a diagram, Fig. 21, that is a combination of the 
‘tandard diagram, Fig. 20, and an old nomogram known as Wallace diagram. 

For using Fig. 21, it is sufficient to extend the straight line in both direc- 
ions joining the value of the discharge Q with the equal number on the scale 
ine. The straight line will directly give the desirable depth, bed width, and 
lope. In choosing the number on the scale-line, a little adjustment is some- 
imes necessary so that the bed width may be a round number, 

Adopting the dimensions given by this diagram, permits a velocity of flow 
laimed to be a critical one. This velocity is said to have the value 


a m 
Vo =PytPyD cee ceerercreceee (82) 


n which 

Py =P. = 0-15 « OR See PAR ca pubetet re BS) 
a = 1,00 and D is the depth of flow when the canal conveys the maximum de- 
aand, 

The authors state that Q could be correlated with D but, for design pur- 
oses, R is a more meaningful measure. This may explain why Q is given in 
2rms of P and R in Figs. 3 and 5, In fact, neither P, R nor the top width is 
f a practical value in fixing the dimensions or proportions of a canal cross 
ection as is the ratio B/D or W/D, and so on. 

The many diagrams given in the paper are, in the writer’s opinion, not 
1ore useful than any one of Figs. 19, 20 or 21. The former diagrams could 
referably have been amalgamated into a single diagram, nomogram or abac, 
) give a direct solution for the problem of fixing the slope, dimensions and/or 
roportions of stable canals under different conditions. 
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Nevertheless, the information given offers a direct relation or proportion 


fi the average width and the depth of stable canals, The width depth ratio can 
ve Obtained as follows: 


From Fig. 2 the perimeter P nearly equals 1.12 W. Noting Fig. 3 and Eqs. 
and 19 


_ 2.668 0.50 
War quae of Dameey Si Lasse eet nae (84) 


ipproximately. 
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FIG, 21.—PROPOSED DIAGRAM FOR THE DESIGN OF TRAPEZOIDAL NON 
SILTING CANALS IN EGYPT WITH 1 ON 1 SIDE SLOPES 


Similarly, taking the values of the hydraulic radius R from Fig. 1 (roughly 
.81 D) and from Eq. 12, from Fig. 5, and Eq. 21, for sand bed and natural 
erm, and from Eq. 22, for coarse non-cohesive materials, that is 


0.43 _ 0.247) 0.361 ure tere (83) 
D= (Gai oF far) ae 


lence, the required ratio may be taken as 


PE Si | ROME oes avin so e(86a) 
or sand bed and natural berm, or 


Asuna. ieee ths no + 02 4 (BBD) 
yr coarse non-cohesive material. : 
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The last kind of material was mentioned, with this description, many times 
in the paper. The writer wonders if there exists a coarse cohesive material. 

One of the important points is the limitation of the validity and applicability 
of regime equations that, as was stated, are valid only within the range of the 
observed data or are only valid for the limited range of conditions on which 
they are based. Two of these conditions were mentioned, in the same page, 
whereas one, a very important one with regard to the range of applicable 
slopes, was ignored. It is true that the range of slopes was mentioned as 
varying from 0.000058 to 0.000388, but this informationis most confusing and 
most misleading as explained by the following example. 

Considering a uniform flow in a stable trapezoidal canal of a bed width ZB; 
depth D and sides sloping at anangle @ to the horizontal, and combining Chény 
and Manning formulas with the Kennedy equation that, in order to facilitate 
computation, may be put in the form 

M= DP p?/s siete a: ssa Nal pat ew a (87) 


instead of the original form of Eqs. 72. Then 
2/3 _ 1.486 2/3 ,1/2 


oe oe BS as age ve. (88) 

D= oN rs/4- aoe 33/4 ( (B+Dcot ¢)D (89) 
pn pn (B +2 Dcosec 6) * 
B+2Dcosec 6 = = s3/4 (B 2D cot 6h ..42% (90) : 
of: : Sat oly - ga25 s°/* p cot 6 -(2cosec @)D.. (91) _ 
Hence ase 3/2 


53/4 (cot 6) - 2(cosec 6) 


3/2 
eu 22888) 33/4 
pn 


If n is given an average value 0.025 and if 6 is considered to be 45° and, 
Supposing that the canal conveys water carrying fine silt such that p in Kennedy — 


equation can be assumed less than o. = or, Say, 0.396 then 
7 
1.486 /4 
“ ( 1.486 ‘- 2 3/4 
-(5 a 


Bes 
D 
-396 x 0.025 


1500 s*/* 2.828 _ s*/*- 0.00188 (93b) 


or 
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Real values of B and D necessitate that the ratio z should be positive that 


means that in the present particular case of water carrying fine silt flowing 
uniformly in a canal of n = 0,025 and @ = 45°, the slope S should vary from 
0.000058 to 0.000232 only, whereas the slope is stated, in page 66, to vary 
from 0.000058 to 0.000389 without limitation or reference to any conditions. 

Even in the “Seaagphaal narrow range of slope between 0.000058 and 


0. ial the ratio 2 = D drops quickly from infinity to zero. Reasonable values 
of = may then be eee only within the range of slopes from about 0.00010 


to a 0.00020, for example, under the described circumstances, 

In addition to the care necessary to design stable canals, measures must 
ibe taken to reduce silt troubles, for example, by feeding the canals and 
branches from the outer side of curved reaches of rivers and main canals, 
‘as well as providing the canals with proper intakes comprising sills, sand 
traps and similar precautions. 


a 4 8s 
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DRAG FORCES IN VELOCITY GRADIENT FLOW2 


Discussion by Donald Van Sickle 


DONALD VAN SICKLE,® A.M, ASCE.—The writer was associated with the 
authors at the planning stage of the investigation and designed and tested the 
force-measuring device that is described briefly in the paper under Mr. 
Moore’s direction, for a previous investigation. ? One of the primary aims in 
ithe development of this device was to make it adaptable for use in subsequent 
investigations, and it is gratifying to the writer that apparently this aim was 
achieved. 

The authors very carefully point out that this was only an exploratory in- 
vestigation, presumably because they felt that the limited amount of data did 
‘not permit more than general conclusions to be reached, There are, however, 
a number of points that could perhaps be given additional attention. 

Eq. 4, in which, presumably, p is pressure intensity in pounds per square 
inch and r the cylinder radius in inches, gives the local force, f, in pounds 
per inch of height. For use in Eq. 3 however, f must be in pounds to give a 
dimensionless cp. It would seem likely that the authors actually multiplied the 
“{” obtained from Eq. 4 by some unit of height, possibly the orifice spacing, to 
determine the local force to be used in Eq. 3. Could the authors clarify this 
point? 

To the writer’s knowledge, the variable shutter gate used for achieving the 
desired velocity gradients has not been reported previously in the literature. 
Could the authors provide more details on the construction and operation of 
the gate and the range of effects achieved? How does it affect transverse 
velocity profiles? Is it possible to achieve a negative velocity gradient, that 
is with the maximum velocity at the channel bottom? Where were the velocity 
profiles of Fig. 4 taken with respect to the gate and the measuring cylinder? 
For how great a distance downstream did the gradient persist before boundary 
layer formation produced a more normal velocity distribution? 

The statement is made that “the measured total drag included, of course, 
the wave drag at the free surface. An approximate evaluation of the wave drag 
was made, based on the profile of the wave surface around the cylinder.” 
Could the authors give more details concerning the procedure for evaluating 
wave drag and, in addition, some quantitative indication of how “satisfactory” 
the agreement was between measured and computed total drag? The total Cq 
values for the various runs would be of interest. 

In the description of equipment the authors mention 24 piezometers in the 
measuring cylinder. Fig. 3 shows curves for only 10 of these, with numbers 


4 July, 1960, by Frank D. Masch and Walter L, Moore (Proc. Paper 2546). 
Hydr, Engr., Turner and Collie Cons, Engrs., Inc., Houston, Tex, , 
7 «an Experimental Investigation of the Drag Characteristics of a Model H-Section 
Bridge Pier,” by Donald VanSickle, Thesis presented te the Univ. of Texas at Austin, 
Tex., in 1958, in partial fulfillment for degree of Master of Science, 
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1, 4, 7, and 14 through 24 being omitted. Similarly, on Fig. 5, the local drag 
coefficient is shown for eight points in one graph, ten points in another, and 
thirteen points in the others. Because the depth was the same for all cases, 
presumably the same number of points would have been available for all. 
Could the authors explain the missing piezometers? 

In Fig. 4, the upper two curves indicate that velocity measurements were 
taken at the water surface, because the depth was stated to be 0.8 ft. Were 
readings actually taken at the surface, or was the depth for these curves 
greater than 0.8 ft? 

The label for ordinates in Figs. 4 and 5 is more correctly given as “Dis- 
tance from the Bottom, Feet” as described in the text. 

Based on the stated flow depth of 0.8 ft and the average velocity of approxi- 
mately 2 fps, the Froude Number for the authors’ tests would be on the order 
of 0.4. In the writer’s experience with the equipment used, Froude numbers 
of this magnitude and higher caused some oscillation of the cylinder, even 
with the dash pot arrangement, Did the authors find that these oscillations of 
the cylinder caused similar pressure oscillations on the manometer bank? If 
so, was any difficulty encountered in taking readings on 13 (or 24) manometers 
simultaneously ? 

In Fig. 1, the arrows indicate that on the upstream side of the pier, at the 
channel bottom, the induced flow isinthe upstream direction. Was any attempt 
made to inject dye from the bottom piezometer to check this? 

Had the authors considered the possibility that the surface wave could be 
at least partly responsible for the induced flows? With a standing wave at the 
upstream stagnation point, and a “standing depression” at the downstream 
stagnation point, the differential head between the stagnation points and the 
surrounding fluid could induce flows similar to those observed. 

In Fig. 3, those portions of the curves lying between 90° and 270° are shown 
as being positive. Because the cosine inthis region is negative, the pressures 
must also be negative to produce a positive product, Are these pressures 
negative with respect to surrounding fluid pressure or some other datum and 
how was the “datum pressure” measured? Because the cosine between 35° 
and 90°, and between 270° and 325° is positive, negative pressures must ob- 
tain between these points to give the negative product shown. Do the authors ~ 
consider the 35° and 325° points (approximately) the points of separation? Do 
the authors consider it unusual that the absolute magnitude of the pressure at 
the downstream stagnation points for curves 10 through 13 is at least twice — 
that at the upstream stagnation points? 

The authors use wave forces on piers and wind forces on stacks as ex- 
amples of engineering problems in which velocity gradients are of significance. _ 
It should be pointed out that in both types of problems pressure fluctuations — 
and flow accelerations are alsopresent. While the authors are evidently aware 
of this, they seem to give the impression that the velocity gradient is more — 
important than acceleration of the flow. Do they feel that their results could, _ 
in general, be said to apply to flow with acceleration even though the experi 4 
ments were performed in a steady, non-uniform flow? “aa 

The writer certainly does not wish to detract in any way from the value } 
of the authors’ work, While he feels that knowledge of the variation in Cq with 
velocity gradient is important, he also feels that knowledge of the variation. 
in cg with acceleration of flow is of at least equal importance. In all probability — 
both velocity gradient and flow acceleration are of significance, and must be 


— ee. 
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msidered together. While several papers®, 9,10, Ilhave been published con- 
‘rning the effect of accelerations on the resistance of cylinders, this is the 
cst paper, to the writer’s knowledge, that is concerned with the effect of ve- 
city gradients on resistance. For this reason it is a most welcome contri- 
ition. 

Although neither the aforementioned papers nor the authors’ paper give 
formation that can be used directly by designers, they do illustrate the im- 
tance of velocity gradients and flow acceleration on drag forces. The writer 
shes to congratulate the authorsfor making the results of their investigation 
ailable to the profession and to express the hope that they have thus stimu- 
‘ed efforts to determine reliable design criteria, 


a Fluid aun b A aca ia ee Motion,” by S, Russell Keim, Pro- 
es wate Pats on Cireulat Cylindrical eee Foiieeni K, E, Beebe, 
ee ee Dinasher by 4. a K, Laird, C, A, Johnson, and 
ater eee aoe Moting Cylinders,” by Alan Dy K. Leind, Charles A. 
ison, and Robert W. Walker, Proceedings ASCE, Vol. 86, November, 1960 p, 43, 
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REDICTING STORM RUNOFF ON SMALL EXPERIMENTAL WATERSHED? 
Discussion by Jaime Amorocho 


JAIME AMOROCHO, 4 F, ASCE.—The author has presented interesting data 
1t show the inherent limitations of the unit hydrograph approach as a method 
predicting the outflow response of a watershed to individual inflow episodes. 
wever, some of the interpretations given in the paper to explain the appar- 
ly irregular behavior of small watersheds in this regard deserve further 
rutiny. 
In accordance with the unit hydrograph principles, the ordinates of any 
rm hydrograph are directly proportional to the ordinates of a unit hydro- 
uph of the same basic duration. The factor of proportionality is the effective 
ensity of the storm expressed as a multiple of the unit rainfall excess. 
is may be said more formally by stating that the unit response to any net 
low episode is an invariant function of the physical characteristics of the 
chment., The term “net inflow” is meant to signify the portion of the gross 
low that becomes surface runoff. The principle of superposition implicit 
this general assumption is extended to any time sequence of inflow rates 
i precludes the possibility of using different unit response functions (unit 
lrographs) as the inflow rate changes. 
In very few if any of the field studies performed to date can it be claimed 
t the unit hydrograph theory has been tested strictly. This is because in 
ictically no instance has it been possible in nature to ascertain with any 
‘ree of assurance that a storm pattern that caused a flood hydrograph 
ually measured, was proportional in all its elements to the storm used to 
‘ive the unit hydrograph under analysis. Other factors, such as the uncer- 
ity in the determination of the infiltration rates used for the construction 
the unit hydrograph and for the predicted flood hydrograph contribute very 
ortantly to the weakness of such tests. Hence, the fact that the approximate 
lines of recorded flood hydrographs have been reproduced in a number of 
tances on the basis of a unit graph previously derived may appear to be 
zely fortuitous. For every instance of approximate agreement there have 
n very many instances of gross discrepancy. It is a well-known fact that 
unit graphs derived from large floods usually differ from those derived 
m minor floods. 
The watersheds investigated by the author are, by virtue of their small 
al extent, particularly well suited for a critical check of the application 
the unit hydrograph concepts. As a direct result of the reduced areas that 
ce the occurrence of very complex geometric patterns of rainfall with 
3p, non-stationary intensity gradients less likely within the catchment 
ndaries, the conditions required for close overall proportionality of in- 
7 can be expected to become more possible, Larger watersheds normally 
USSR Sil ht LE a? SATIS GIF PII a 


| August, 1960 by Neal E. Minshall (Proc, Paper 2577). 
L Acting Asst. Prof., Univ. of California, Berkeley, Calif. 
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produce an effect of attenuation of the rainfall variability that occurs at dif 
ferent locations, provided no major discontinuities exist in the storm patter 
over the catchment. In these cases the outflow hydrograph becomes rathe 
insensitive to inflow intensity changes recorded at individual precipitatio 
stations, and its shape tends to reflect the catchment characteristics only 
However, because major rainfall discontinuities are more likely in extensiv 
watersheds, the conditions conducive to proportionality of inflow are, in ef 
fect, rare. 

The author has derived unit hydrographs on the basis of floods produce 
by short storms of different, “reasonably uniform” intensities, The hydro 
graphs shown in Fig. 6 were presumably reduced to the same unit rainfa. 
excess duration in order to render them consistent with each other, This i 
not indicated in the paper, but because the observed rainfall excess durations 
as reported in Table 3 are approximately the same, no large error is to b 
expected if the conversion was omitted. The striking fact is that the hydro 
graphs represent distinctly different outflow responses. 

A basic premise of the unit hydrograph procedure is that the relationship 
between storage and outflow are linear.5 This results from the requirement 
for a linear solution with respect tothe net inflow rate of the storage equatio 
for a hydrologic unit. To illustrate, consider abasin receiving a constant rat 
of inflow. 

The storage equation can be stated as 


in which I is the inflow rate, Q represents the outflow rate and dS/dt denote 
the rate of change in storage, This equation has an immediate linear solutio 
in I provided S and Q are related linearly. For if we can write 


Si K: Qis. ores Gland ett 3% Hist ee 
K being a constant, then Eq. 2 becomes 


FeQ UE oe susiest ot eae 
that integrates to 3 


Q= 1(1 -e rad 
In this equation, the coefficient a = 1/K in the exponent is a characteristi 
parameter of the basin. The Sonstant K has the dimension of time and rep 
sents the lag due to a “linear reservoir” as defined by Eq. 3. The equa 
represents the rising limb of the instantaneous unit hydrograph. The equatic 
of the recession limb follows from integration of Eq. 4 for I=0att=t 


Q=ae 2 1S) can entalt od cee eee 
The subindex refers to the values at the moment of cessation of the inflo’ 
Now, because by Eq. 5 the discharge at peak is 


a es q 
nial werent Q,=10 e 2 eee 
Bet (eats) 25° 4) errant a 


It is seen that Eqs. 5 and 8 are linear jioaie 3 of I, that is the bas 
premise of the unit hydrograph method. J. C. I. Dooge,5 M. ASCE, has mas 


5 Journal of Geophysical Research, J. C, I. Dooge, Febnuary 1959, 
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exhaustive analysis of the unit hydrograph theory in which variability of 
inflow rate is accounted for. 
In natural watersheds, the assumptions of linear reservoirs and linear 
nslation processes can not be taken for granted. The terms S and Q are 
‘mally related by non-linear functions which contain additional variables 
+h as friction and velocity gradients. Introducing also the space variability 
she inflow rate, the storage equation has the general form 
I(t, x, y) =Q+ in £(Q, Ss, X, J, z) penta ee (9) 

It can be easily surmised that this equation will, in general, not have a 
2ar solution in I, 
Practically all cases in nature are of this complexity; therefore, it may be 
ted quite generally that the classical unit hydrograph concept can almost 
rer be applied strictly. It shouldbe emphasized that the preceding conclusion 
rerns not only the proportionality between the ordinates of the runoff hydro- 
yph and the unit hydrograph but the additive superposition of the effects of 
“iable inflow rates as well. 
Some rather interesting conslusions may be derived from Figs. 5 and 6, 
t illustrate the previous statements. The derivation of unit hydrographs 
m storms of varying intensities give rise to inaccuracies that depend on 
yjective factors that inevitably creep into the process of separation of the 
ferent component runoff episodes. Inaccuracies also result, as mentioned 
viously, from the determination of the infiltration rates, that, in the pres- 

case, was made on the basis of a completely empirical criterion. Never- 
less, some measure of verification of the latter was afforded by the ob- 
‘ved actual runoff, and the clear trends apparent in the alignment of points 
tted in Fig. 5 appear reliable. From these trends it may be concluded (1) 
t the shape of the unit graph is a non-linear function of the inflow rate; and 
that depending on whether the unit graph corresponds to an initial period 
inflow or to inflows that occur after some runoff created by previous rain 
taking place, the response is different although the intensity may be the 
ne, This may be interpreted as a “field proof” of the inapplicability of the 
nciple of superposition in this case. 
It is obvious that the exact relationships suggested by the data of Fig. 5 
ly only to the particular watershed described inthe paper. Other relation- 
ps will exist in different catchments, but as. long as their empirical deter- 
1ation is possible, the type of analysis that follows should have general ap- 
sation. The only requirement is that the overall geometry of the inflow pat- 
n should not vary greatly within the basin boundaries for different rain 
snsities. 
It is noted that the lines of the upper and lower graphs of Fig. 5 can be 
ss-plotted to define a direct relationship between the hydrograph peak and 
time to peak. Scaling values from the figure, it is found that the following 


ressions result: he 
For unit response early in storm: Qa = 15.25 ty : a wt(l0) 


- 0.793 
For unit response late in storm: dmax = 28-55 t) ea) 


which dmax is the maximum ordinate of the unit response and t, is the 
e to peak. These equations represent the loci of the peaks of hydrographs 
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such as those shown in Fig. 6. These hydrographs may be considered as 
families of curves that must satisfy the following conditions: 


1. To have maxima passing through the lines described by equations such 
as Eqs. 10 or 11. 

2. To encompass an area of unity between each curve and the time axis 
(by the definition of the unit hydrograph). 

3. To have 0 value at time equal zero and to be asymptotic to the time 
axis. 

4. To be continuous and have a positive second derivative near the origin 
and a first derivative equal to zero at the origin. 


These conditions are met for example by some asymmetric frequency dis- 
tributions. It may be noted by reference to Fig. 6 that the curves tend to be- 
come symmetrical as the time topeakincreases. We might, therefore, choose 
a position in which this symmetry is closely approximated and fit therein a 
normal frequency distribution to represent the unit response. The ordinate 
at time zero would have a finite value, but the error thus introduced in meet- 
ing conditions 2 and 4 would probably be small. This basic function could 
then be skewed and peaked simultaneously by introducing appropriate para- 
meters in such a way that Eqs. 10 or 11 would be satisfied for each value of 
tp. A similar approach could conceivably be used with other functions, in 
order to arrive at some empirical expression for the unit response as a func- 
tion of the inflow intensity. 

It is desirable, however, to use a function that is meaningful in terms of a 
theoretical analysis of watersheds, J. E. Nash® and Dooge® derived inde- 
pendently the general expression of the outflow function for a catchment with 
n linear reservoirs and a corresponding number of linear channels. This ex- 
pression has the form of a Poisson distribution and represents the unit hy- 
drograph of a linear catchment of this type: : 


q (t) = 


in which a = x has the same meaning as in Eq. 5. 


aT SPE 


Toa wt A 


The a anat of a and n can be computed from observed outflow hydrographs 
by the method of moments as described by Nash. A simple procedure, tha 
assures that the maximum value of q satisfies Eqs. 10 or 11 exactly, although 
the fit may be less precise elsewhere, consists of finding the solutions for a 


and n from the system formed by Eq. 12 and the relationships 3 
r 
Gag ee Ps. ee Sa a Sy 
and rd 
ty Pbepth on alsl ite he eifsk ser due Sue aee tees Ce 


obtained from the data of Fig. 5. The effect of the imperfect fit of the equatio or 


to the unit response derived from the data is of little practical consequence 
as will be discussed later. 


For q maximum, the first derivative of Eq. 12 must be zero: 


n a 
; ee n-2 -at n-1 -a 3 
A Sec Fer |o-0) neailie a picasa ties os »| “ohh 


6 Proceedin ings, J. E, Nash, Genl. Assembly of Toronto, Internatl, Assn. of Scientifi 
Hydrology, September, 1957, 
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q’ te qh et th ( 1) t 2-2 n-1 
i Tn) Th i= p = ty =a eee cee (15b) 
tence 
(n= 1) =a 4 és SEC dea rorel e ae (16) 
combining Eqs. 16, 12, 13 and 14 yields 
4 n (n- 1)" ein) e r 
max Pra (7) em iene = 6 I eoeeeeeee (17) 
“igh ee; 
(n-1)2e RE 25 6 A- AY ; 
OE a“) = 1 | Oe ee eee Pees eee ee (18) 


explicit expression for n in terms of the experimental parameters 5, 6,, 
id Ay can be obtained from Stirling’s approximate formula for Jn): 4 


Pn) ~ (2x)? (a) YM? exp, | -(n-1) + EG] +--+ 9) 


, when substituted in Eq. 18, yields, 


1 
poe fo, 0852 ee A ce (20) 
2<n< 9 it is found that 
ai 
rpm). (] Clee Ree ey oa: Be’ ei (21) 
refore, Eq. 20 becomes 
n=27 62 6,2 pe e917... (22) 


ewise, from Eqs. 13 and 14 we get 


a=27 0% 6, Pr-M 4, OEE st eee (23) 
1 


‘rom the data presented in the paper we find that 6 = 1.170; 54 = 0.483; 
0.710; and Aq = 0.547 as a result of which 


FeO Toso ack da LT. seanccscessnanceea ees ee 2S) 


a = 4.1569 2-873 + 0.3518 tie a-Si se 125) 


‘ig. 11 shows a comparison between the hydrographs presented by the 
ior in Fig. 6 and the curves computed by means of Eqs. 12, 24 and 25. 
t is noted that the fit becomes increasingly better for higher values of the 
ow intensity. Similar fits computed by the method of moments were found 
ye slightly better for the curves corresponding to low intensities but be- 
.e increasingly difficult at higher intensities. The precision of the fit is 
lin the order of magnitude of the error to be expected from the field data, 
nay be surmised from the scatter of the points shown in Fig. 7. In general, 
htly higher runoff volumes would correspond to the areas below the fitted 
yes during the first hour than to the areas below the curves shown in Fig. 


é 
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Eqs. 20 and 23 suggest that the hydrologic behavior of a watershed may be 
characterized by functions that define the number of equal theoretical linear 
reservoirs and the theoretical lag time of one linear reservoir in terms of 
the inflow rate, The substitution of these expressions in Eq. 12 gives an equa- 
tion that can be considered as the approximate solution of Eq. 9 when I does 
not vary in time, An interesting implication of these equations is that a water- 
shed may be regarded as equivalent to a linear catchment only when the in- 
flow is constant. As soon as the inflow rate changes, the number of theoretical 
reservoirs becomes larger or smaller and the response varies accordingly. 
Therefore, neither n nor a may be considered as invariant characteristic 
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parameters of a catchment. This conclusion explains the inapplicability of the 
principle of superposition in this type of hydrologic problem, and may no 
only be extended to the critical evaluation of the unit hydrograph method bu 
to aaa related flood routing procedures, such as the Muskingum method, as 
well. ; 

The preceding analysis suggests a practical procedure for obtaining the 
response function for small watersheds in which uniform precipitation rate: 


> 
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are recorded. This procedure would consist of first selecting a number of 
individual runoff hydrographs that can be separated and definitely ascribed 
to periods of reasonably uniform rainfall excess. The values of peak discharge 
and time to peak should then be plotted against rainfall intensity and from 
these relationships, by an analysis similar to that described in this discussion, 
the expressions for n and a in terms of I could be derived. It should be kept 
in mind that whatever form these expressions may acquire, their nature is 
strictly statistical and empirical. Hence, from a practical standpoint there is 
no object in scrutinizing the physical meaning or the dimensions of the para- 
meters that enter their configuration. 

An investigation is in progress (in 1961) at the University of California, 
Berkeley, Calif., in which these concepts are being studied in model catch- 
ments. A full theoretical and experimental treatment of this subject, initiated 
in 1959, is being carried out under a program of the University’s Water Re- 
sources Center. 


‘“EDICTING STORM RUNOFF ON SMALL EXPERIMENTAL WATERSHEDS2 
Discussion by Merwin D. Dougal 


MERWIN D. DOUGAL, ! A. M. ASCE.—Research and experience have well 
licated the many variables that must be considered in developing a usable 
nfall-runoff relation. Many that have been developed present unique prob- 
us in their use in estimating runoff on ungaged watersheds, The author’s 
rk extends known techniques both in methods of approach and in providing 
ormation on a particular area of the nation. 
The writer has studied8 the problem of a modified infiltration index and the 
-e of infiltration rather than either total retention or runoff estimates, but 
is apparent that either of the latter could be obtained indirectly. The ante- 
jent precipitation index approach was used, and season (or week-of-year) 
ration of excess rainfall and average intensity of rainfall were the other 
lependent variables given primary consideration. Certainly the Ralston 
eek watershed near Iowa City, Iowa is classical in hydrologic research. 
The use of the antecedent precipitation index, computed by Eq. 1, appears 
be widely applicable. As recommended by Kohler and Linsley,9 the index 
nuld be used in conjunction with season of the year or temperature. The 
hor states, for the various methods and lengths of period tested in com- 
‘ing the index, that the best correlation with retention was a 30-day index 
is 20% of the departure from normal of the 30 to 60 day precipitation. It is 
er stated that an antecedent retention index was used in an attempt for re- 
ement, but the results showed no noticeable improvement over the 30-day 
lex finally selected for use. The recession factor “k” in Eq. 1 was evaluated 
0.95 for the Edwardsville, Illinois watersheds. 
The selection of 0.95 for the recession factor, a relatively high value, in- 
ates that previous precipitation affects soil moisture in the watershed to a 
ystantial and somewhat lasting degree, because the API values recess quite 
wly. This means that if appreciable precipitation has previously occurred, 
effects on soil moisture linger on. Perhaps the author might comment on 
effect of the impervious clay-pan layer onthe relatively low amount of soil 
isture capacity indicated. 
The writer initially investigated three values of the recession factor “k” 
the Ralston Creek watershed: namely, 0.84, 0.90, and 0.94. In the result- 
graphical correlation, the values of 0.84 and 0.90 were about equally ca- 
le in enabling first approximation month-of-year curves to be drawn, but 


7 Staff Engr., Iowa Natural Resources Council, Des Moines, Iowa. 
8 “A Study of Factors Which Affect Infiltration Rates on the Ralston Creek Water- 


.d,” by Merwin D, Dougal. Thesis presented to the Graduate College of the State 
iversity of Iowa, at Iowa City, Iowa in February, 1958, in partial fulfillment of the re- 
rements for the degree of Master of Science. 

9 «predicting the Runoff from Storm Rainfall,” by M. A. 
search Paper No. 34, Dept. of Commerce, Weather Bureau, 
, 1951, 


Kohler and R. K, Linsley, 
Washington D.C,, Septem- 
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little correlation was evident for the value of 0.94. The mean of the first two 
values, 0.87 was subsequently selected as most applicable for the Ralston 
Creek watershed. This watershed is about 3.0 square miles in area, relatively 
rolling to hilly. The soil is of loessial origin over glacial till with a Clinton 
silt loam soil surface. The Ralston Creek watershed experiences high infil- 
tration rates. 

In regard to the methods and lengths of period available in computing the 
antecedent precipitation index, the writer would like to add several others. 
The problem is selecting an initial API value and a proper time interval be- 
fore a storm in which to recess the initial index. The author does not state 
what initial value was used in his study. The writer investigated four proced- 
ures in the Ralston Creek study: 


1. Beginning with a zero value of API 60 days prior to the storm. 

2. Beginning with the normal 10-day precipitation for the period prior to 
30 days before the storm (first storm on May 1, use as the initial value the 
10-day normal precipitation prior to April 1, and begin API computations on 
April 1). 

3. Beginning with the observed average 10-day precipitationfor the 30 day 
period prior to 30 days before the storm (first storm on May 1, use as initial 
index the 10-day average precipitation of March and begin API computations 
on April 1). 

4, Same procedure as (3), except the 10-day average precipitation is prior 
to a 3-week period before the storm, 


The four procedures, when applied in the Ralston Creek study to periods 
prior to storms indicated below, yielded the values indicated in Table 5. The 
results indicate that the initial value need only be reasonably assumed, To 
illustrate, for the storm of August 4, 1924, the initial value for procedure (2) 
was 1.58 in., the 10-day normal for June, For procedure (3) the observed - 
10-day average precipitation for June was 2.47in. As indicated previously, the — 
respective API values after a 30-day recession period were almost identical. 
It seems possible to the writer that the use by the author of a different “k” 
factor in a 30-day index might have givenidentical results to the 30-day index 
plus the 20% departure used in the study. 

A comment by the author on his method used in computing API values for © 
successive storms within a day would be appreciated. For instance, in Table : 
1, the same API valueis listedfor three storms occurring on August 14, 1946, 
Presumably the API value is higher for successive storms on the same day 
either by direct addition, or by recessing on an hourly basis. : 

The writer concluded, in the analysis of API in the Ralston Creek study, 
that it is simpler to carry the API computations forward continuously from the 
first storm in a particular year to the last storm, Less computations are 
involved than if separate API computations are made for each storm, In addi- 
tion, this continuous method would be used in flood forecasting oa 

The writer agrees that a large volume of data is needed to develop graphi- 
cal correlations in a hydrologic analysis. Over 258 storms in a 1927 to 1956 
record were used by the writer and all were utilized to advantage. However, 

about 100 of the storms were selected as having “ideal” storm characteristics 
and were used in initial correlation efforts of API and season, month-of-year, 
or week-of-year. 

The author indicates the variation inseason is season that can be expected. 
Because the final grouping refined some periods to 30-day periods, it is of 
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interest to inquire of the author’s method for allowing for the seasonal vari- 
ation in plotting the data. The writer attemptedto use temperature as a para- 
meter (departure from the mean weekly temperature for the week prior to the 
storms under consideration). It was concluded that if all other variables were 
about constant, temperature effects were of importance, but normally were 
too obscured by other variables to be of value. 

The writer believes that the abscissa in Fig. 3 could be more aptly termed 
“Duration of excess rainfali, in hours” because the plotted points indicate a 
single total storm event and not the variationin retention during a storm, The 
author also assumed that infiltration capacity is independent of rainfall in- 
tensity. This is a perplexing problem to analyze because in cultivated areas 
the intensity may effect surface cover conditions, especially in the spring 
months. The comparison of estimated and observed runoff indicated in Fig. 4 
shows best results for runoffs above one inch, and it is above this that the 


TABLE 52 


Date of Storm Procedure 


August 4, 1924 


March 14, 1934 0.16 

May 3, 1946 0.81 

May 28, 1947 0.80 
*k = 0,87 


@ API (k = 0.90) 


engineer is normally concerned in estimating storm runoff, For use in esti- 
mating past yields, individual estimates might be in error, but an acceptable 
average should be obtained over a period of days or weeks. 

The writer, in the Ralston Creek study, was unable to correlate certain 
storms that appeared on further analysis to possess a common characteristic, 
that is, a fair amount of intense rainfall had occurred within a few hours (24 
hr period maximum) before the beginning of the storm being scrutinized. This 
was finally attributed to an additional variable, namely, the variation or dif- 
ference between a rising and falling API at the time a storm occurs, It is 
easily reasoned that there is considerable difference in soil surface conditions 
whether the API is falling or rising to the index existing immediately before 
occurrence of the storm rainfall. Toillustrate, assume an API value of 2.0 in. 
exists immediately prior to each of two hypothetical storms. For one storm 
the API of 2.0 in. has resulted from the log arithmic recession of a 5 in. to 
6 in. rainfall a week or more before the storm in question. For the second 
hypothetical storm the API of 2.0 in. resulted primarily from a 2 in. intense 
rainfall a few hours before the storm in question. For these two hypothetical 
storms, all other conditions being the same, the obvious difference in soil 
moisture at the surface and in soil surface conditions will cause infiltration 
values for the two storms to be much different. The writer introduced shift- 
curves to account for this phenomena that enable, in effect, a shift to be made 
to a higher API value when arising API situation was encountered, thus giving 


a lower infiltration rate. 
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The author has accomplished an admirable evaluation of hydrographs on 
small watersheds. The variation of the peak rate of runoff with both rainfall 
intensity and pattern of rainfall intensity is of interest. In regard to the vari- 
ation with intensity, is it possible that the high intensity causes surface depths 
to reach a point at which not only is the flow velocity increased appreciably, 
but surface vegetation is altered; that is, flattened, disrupted, and so on? The 
author’s results, as indicated in Table 4, reveala good correlation, with most 
values for large amounts of storm rainfall being within 10%. 

The writer would like to inquire if the author compared his synthetic hydro- 
graphs with those that would be computed by present techniques used by the 
Soil Conservation Service and developed by the Agricultural Research Service, 
U. S. Department of Agriculture. An examination of Fig, 9 and Fig. 5 indicates 
that the former is based on the occurrence of high intensity rainfall late ina 
storm. Perhaps the author can state which of the two effects, early or late, 
is more commonly experienced in nature. 

The results of the author’s study are stimulating. Additional research ef- 
forts in other parts of the nation should be encouraged to aid the hydrologist 
and the engineer in his design capacity. Arecompilation and analysis of Ralston 
Creek data has been completed at the University of Iowa, Iowa City, Iowa, and 
is to be published (1961) in cooperation with the Iowa Highway Research 
Board and other agencies. This report should also add to the field of knowledge 
concerning small watersheds. 


ROBERT L. McFALL!° and BENA. JONES, JR.11_ The validity of the pro- 
posed method was tested with data from a central Illinois watershed main- 
tained by the Agricultural Experiment Station of the University of Illinois. The 
watershed is located near Monticello, approximately 120 miles northeast of 
Edwardsville, Illinois. From the limited data available, two storms over the 
82-acre watershed were selected for comparison. 


TABLE 6,—STUDY OF CENTRAL ILLINOIS WATERSHED 


Date of Peak Runoff Rates, in | Total Runoff, in inches Excess 
' Storm inches per hour precipita- 
observed | computed | observed computed tion, in 
inches 
(1) (3) (4) (5) (6) (7) 
June 27, 1951 0,51 0,58 0,56 0,50 
July 9, 1951 1,48 0.68 0.83 0,82 0.91 


The watershed is predominantly silt loam soil and has slopes that are 
generally less than 4%, typical of the dark-colored, moderately permeable 
gently sloping prairie soils of east central Illinois. : 

The time retention curves proposed by Minshall were used in computing 
runoff rates. The 30-day antecedent precipitation index (API) was computed. 
using the same formula used by Mr. Minshall except the K value was changed 
to 0.90. (There is some question in the minds of the writers as to whether a 
30-day index should be used, but data are lacking to test this parameter.) 


10 Instr., Univ. of Ilinois, Urbana, Ill, 
Assoc, Prof., Univ. of Illinois, Urbana, Ill, 
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The results of the study are shown in Table 6 and one of the hydrographs 
is presented in Fig. 12. Table 6 shows that the observed peak runoff from 
both storms was less than the computed peak runoff. But the observed total 
runoff was higher than the computed value for one storm and the reverse was 
true for the other storm. A study of the hydrograph shows that a change in 
the shape of the portion of the hydrograph that falls below the 15% runoff rate 
could change the volume of runoff considerably. Thus, the total runoff is a 
weaker comparison tool than the peak discharge rate and the time-rate change 
of the hydrograph. 

Because the time retention curves presented by the authors were computed 
from data from clay-pan prairie soils, the use of a soil retention factor, in 
addition to the seasonal parameters, is suggested. Although the time retention 
curves provide a means for determining the runoff, their use requires six 
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families of curves for each variation in soil retention characteristics. The 
data and retention curves presented were studied, and it appears feasible to 
construct a coaxial correlation that would contain, in addition to the seasonal 
parameters, information about the soil characteristics. This proposed method 
will greatly facilitate the determination of precipitation retention when working 
with widely varied soil conditions. More data and exploration are needed to 
confirm this opinion. 

The tests conducted indicate that the method proposed by Minshall for syn- 
thesizing hydrographs will lead to peak runoff values that are higher than ob- 
served flows. It is impossible with the limited data available to predict whether 
total runoff will be higher or lower than observed values, In these cases the 
difference in values does not exceed 15%. These are good results considering 
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the limited knowledge about some of the factors. The writers believe Minshall’s 
method shows promise and its applicability to a number of soil, cover, an 
climatic conditions should be investigated by other workers. 3 
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